Method for measuring the concentration of trace gases by scar spectroscopy

ABSTRACT

The present invention is relative to a method of ring-down spectroscopy in saturated-absorption condition, for measuring a first concentration of a gas through a measurement of the spectrum of a molecular transition of said gas.

TECHNICAL FIELD

The present invention relates to a method to analyse saturated-absorption cavity ring-down (which hereafter will be referred to as SCAR for brevity) signals, detected with a SCAR spectroscopy apparatus, for measuring with high accuracy the concentration of trace gases, in particular for trace gases present in part per trillion (10⁻¹² or ppt). Radiocarbon detection at concentration levels much lower than those present as natural abundance is an example of the applicability of this method.

TECHNOLOGICAL BACKGROUND

The so-called method of ¹⁴C (carbon-14 or radiocarbon), is a radiometric dating method based on the measurement of relative abundances of carbon isotopes. The method of ¹⁴C allows dating organic materials (bone, wood, textile fibers, seeds, wood, coals . . . ), thus containing carbon atoms.

Carbon is a chemical element essential to life and present in all organic substances. It is present on earth in three isotopes: two stable ones (¹²C and ¹³C) and a radioactive one (¹⁴C). The latter turns by beta decay into nitrogen (¹⁴N), with an average half-life of 5730 years, thus this isotope would disappear in the long run, if not continually reinstated. The production of new ¹⁴C regularly occurs in nature in the high layers of the troposphere and in the stratosphere, by the capture of thermal neutrons, secondary components of cosmic rays, by the nitrogen atoms present in the atmosphere. The dynamic balance between production and radioactive decay then keeps the concentration of ¹⁴C in the atmosphere constant, where it is mainly bound to oxygen in the form of carbon dioxide.

All living organisms that are part of the carbon cycle are continuously exchanging carbon with the atmosphere through breathing (animals) or photosynthesis (plants) processes or they assimilate it by feeding on other living beings or organic substances. Consequently, until a body is alive, the isotopic ratio of ¹⁴C and that of the other two carbon isotopes remains constant and equal to what is found in the atmosphere. In particular, the current natural isotopic ratio (abundance) in the atmosphere is:

$r_{0} = {\frac{14_{C}}{C} \approx {1.2 \times {10^{- 12}.}}}$

After death, these processes end and the organism does not exchange carbon with the outside anymore. Then, as a result of decay, the isotopic ratio decreases on a regular basis according to the formula:

r=r ₀ e ^(−Δt/τ),

where Δt is the time elapsed from the death of the organism and τ is the average lifespan of ¹⁴C.

By measuring the amount of ¹⁴C present in organic remains, the age thereof is obtained by applying the following formula:

Δt=−τ ln(c/c ₀).

The measurement of ¹⁴C is possible, given the low concentrations present, with the method of mass spectrometry (AMS, Accelerator Mass Spectrometry): using a mass spectrometer, the concentration of ¹⁴C present in the sample is measured. This method is able to obtain reliable measurements for concentrations in the order of

$\frac{14_{C}}{C} \approx {10^{- 15}.}$

However, the costs related to AMS equipment are relevant, in the order of millions of euros, and the overall dimensions of the same equipment is substantial, with high operating voltages. What above stated as an example stands for Carbon as well as any other substance which is present in traces and need to be measured.

Additionally, the detection of trace gases is generally very relevant in various technology fields. Apart from dating, the measurement of the amount of radiocarbon is important in biomedicine or in environmental and earth science.

A class of techniques for measuring the concentration of a gas is spectroscopy. Spectroscopy is a scientific technique that analyzes the spectrum of electromagnetic radiation emitted/absorbed by a source split in its wavelengths and hence it analyzes the properties of atoms or molecules that are the source/sink of such radiation. In these spectra, the lines of absorption or emission can be studied.

The origin of a given spectral line can be an electronic, vibrational or rotational transition of the molecule of interest. For example, in the infrared, the main origins of a spectral line are not transitions between energy levels of electrons, generally dominating in the visible spectrum, but transitions between molecular vibrational energy levels.

The conventional cavity ring-down (CRD) spectroscopic technique was devised over 25 years ago, using pulsed laser first and then continuous emission laser. The advantages of CRD are mainly two, as detailed below: the signal is immune to amplitude fluctuations of the radiation source used; the linear absorption coefficient is measured directly and then, knowing the total pressure of the gas and the line-strength of the spectroscopic transition, the concentration of the molecular species to be measured.

Typical CRD spectroscopy includes an apparatus comprising a laser that is used to send a highly coherent radiation beam to a high-finesse cavity consisting for example of two highly reflective mirrors (for example with a reflectivity R>99.9%). When the radiation emitted by the laser at frequency ν is in resonance with a cavity mode, the radiation intensity increases in the cavity due to constructive interference phenomena. The laser is then quickly switched off, or moved away from the resonance cavity, in order to measure the exponentially decreasing intensity of light that escapes from the cavity. During this decay, the light is reflected thousands of times by the mirrors, following a round trip, yielding an effective path a few kilometers long.

If a gas or a mixture of gases that absorb light is placed inside the cavity, the intensity of photons trapped decreases by a fixed percentage along each path inside the cavity due to the absorption from the medium in the cavity and due to reflectivity losses. The light intensity inside the cavity is then given by an exponential function of time:

I(t)=I ₀exp(−t/τ).

CRD spectroscopy measures how long light employs for its intensity to decay to 1/e of its initial intensity value, and this value of ring-down time τ is used to calculate the concentration of the absorbing substance in the gas inside the cavity.

The operating principle is thus based on the extent of a decay rate rather than an absolute absorbance. The decay constant τ is called “ring-down time” and is dependent on loss mechanisms inside the cavity. For an empty two-mirror cavity, i.e. without an absorbing medium inside, the decay constant τ₀ is dependent on mirror losses (transmission, absorption and scattering) and various optical phenomena such as diffraction:

$\begin{matrix} {\gamma_{c} = {\frac{1}{\tau_{o}} = {\frac{c}{l}\left( {1 - R + X} \right)}}} & \left( {{Eq}.\mspace{14mu} 1} \right) \end{matrix}$

Where γ_(c) is the decay rate without absorbing gas, in ms⁻¹, c the speed of light in vacuum in cm/ms, l is the length of the cavity in cm, R the reflectivity of the cavity mirrors, and X takes into account various other optical losses different from gas absorption. The equation uses the ln(1+y)≈y approximation for y close to zero, which is the case in the working conditions of the CRD.

A gas inside the cavity absorbs energy by increasing the losses according to the Beer-Lambert law and then the intensity of the light decays more quickly, after the laser has been switched off. Therefore, assuming that the gas fills the entire cavity, the decay time becomes:

$\begin{matrix} {{{\gamma \left( {v - v_{o}} \right)} = {\frac{1}{\tau} = {{\frac{c}{l}\left\lbrack {\left( {1 - R + X} \right) + {{\alpha_{g}\left( {v - v_{o}} \right)}l}} \right\rbrack} = {\gamma_{c} + {\gamma_{g}\left( {v - v_{o}} \right)}}}}},} & \left( {{Eq}.\mspace{14mu} 2} \right) \end{matrix}$

where τ is the cavity ring-down time due to all intracavity losses, in ms, γ(ν−ν₀) is the decay rate with absorbing gas in the cavity, in ms⁻¹, γ_(g)(ν−ν₀) is the absorbent gas contribution to the total cavity decay rate in ms⁻¹,

${\alpha_{g}\left( {v - v_{o}} \right)} = {\frac{1}{c}{\gamma_{g}\left( {v - v_{o}} \right)}}$

is the absorption coefficient of the transitions of the tested specific gas; ν₀ is the resonance frequency of the absorbent molecular transition, in Hz; ν is the absorbed light frequency in Hz. This is called a linear regime as α is considered independent of the intensity of the radiation. γ_(c) is considered independent on ν for frequencies close to the resonance frequency ν₀.

In other words, the cavity ring-down event occurs by abruptly stopping the radiation from the laser that impinges on the cavity and is characterized by a power transmitted which decays according to the exponential function exp(−γt), where γ=1/τ and t is the time measured from the moment of interruption of the incident wave.

If the cavity has an internal linearly absorbing medium, the constant γ is simply the sum of two terms: γ_(c)=1/τ₀ representing the empty cavity contribution and γ_(g)=1/τ−1/τ₀ which represents that from the medium absorption. Therefore, with two separate measurements, one with empty cavity and one with absorbing medium, the value of γ_(g) can be ideally determined.

Recently, a new technique was presented related to laser spectroscopy, called “saturated-absorption cavity ring-down spectroscopy” (hereinafter briefly SCAR), described in “Saturated-absorption cavity ring-down spectroscopy” written by G. Giusfredi et al., Phys. Rev. Lett. 104, 110801 (2010), which has proven that high sensitivity can be achieved. This technique is called below ring-down spectroscopy under saturation of absorption. This means that the intensity of radiation in the cavity that is set is much greater than the saturation intensity of the molecular transition to detect.

SCAR spectroscopy uses a non-linear effect bringing the absorbing medium to saturation, i.e. the intensity of the laser beam is such as to lead the molecular transition of the gas of interest—resonant with the laser—to saturation. In other words, the wavelength of the electromagnetic radiation emitted by the laser is adjusted so that it is in resonance with the transition of interest and the intensity of the radiation itself is increased or adjusted so that this transition is brought to a saturation condition. From the studies reported by the Applicants, they have figured out how to take advantage of the fact that a cavity containing a gas in high saturation conditions behaves almost as an empty cavity in relation to radiation, i.e. when the laser is switched off, the emission of photons follows a curve similar to that of the empty cavity at least for a first time interval. This scheme is called “effective empty cavity scheme”. Therefore, in an experiment where saturation is reached inside the cavity using a beam having a sufficient intensity and then turning off the same, measuring the radiation emitted with a photodetector, a curve is obtained that for a first part follows the decay pattern in an empty cavity. After a certain period of time, however, the behavior of the radiation emitted is no longer that of an empty cavity, since many photons have already left the cavity that contains a gas that is no longer under saturation conditions, thus for a second time interval the decay curve is the curve that one would get if the cavity was filled by a non-saturated gas, i.e., one gets back to the linear absorption regime. Thus, by measuring the ring-down radiation transmitted by the cavity, the two decays γ_(c) and γ_(g) are measured independently, using both the saturation condition and the linear condition from the same decay event. In other words, the value of γ_(g) is encoded in the “deformation” of the exponential decay.

Since in principle all information about the gas absorption to be obtained is contained in the “non-linearities” of the decay event (very small in the case of trace gas sensing at the ppt level or less), the SCAR spectroscopy minimizes the following errors that are introduced in the detected decay signals and cause not to be in an ideal condition, ultimately preventing two measurements as in conventional CRD spectroscopy:

-   -   the non-monochrome condition of the wave emitted by the laser         which is incident onto the cavity;     -   the imperfect immediacy of the interruption of the wave to be         “turned off” to measure the ring-down time;     -   the fluctuations of the resonant frequency of the cavity;     -   the imperfect matching of the spatial mode of the incident wave         emitted by the laser to the cavity mode, which can also vary         overtime;     -   the reflectivity inhomogeneity of the mirrors forming the         cavity, which combines with the alignment fluctuations of the         incident wave;     -   the dependence of γ_(c) on frequency, when going from “full         cavity” to “empty cavity” different longitudinal resonant modes         of the cavity are used, one coincident with a region of         absorption of the medium (i.e. the gas introduced into the         cavity the concentration of which is to be measured) and one in         a region of transparency;     -   the fluctuation of this dependency over time; and     -   The non-linearities of the detector and detection electronics in         a dynamic range of several decades of the detected decay.

In order to give a quantitative idea of the resolution obtainable through the SCAR technique, let us consider the special case of radiocarbon dioxide detection in natural abundance. In optimum conditions of temperature and pressure, the deformation from a pure exponential produced by radiocarbon along decay signals is of the order of 1 μV out of 3 V. This set very stringent limits on the residual non-linearity which can be born.

The measurement is possible thanks to the noise present in the signal to be captured, which, in good approximation, has zero, or at least constant, average. By averaging many events (each decay signal stored in memory is the average of 1280 events), it is possible to increase the resolution of digitization by approximately 35 times.

SUMMARY OF THE INVENTION

The invention relates to a method for the measurement of the concentration of a gas on the basis of the SCAR method outlined above, i.e. by measuring the temporal evolution of the ring-down in a condition of saturation of the transition of interest.

As indicated above, the SCAR spectroscopy methodology is based on a fit of a curve of the ring-down signal detected in conditions of resonance between the electromagnetic radiation emitted by the laser and the cavity and between the electromagnetic radiation emitted by the laser and the gas transition that has distortions with respect to a perfectly decreasing exponential, exponential which represents the decay curve in linear absorption condition. Just because a saturation condition is used, all the measurement is based on non-linear absorption of trace gas located within a high-finesse cavity containing the gas itself. The exponential curve is thus distorted. The extent and shape of the distortion contains information on the absorption that the gas would have in linear condition and thus on its concentration.

However, retrieving the linear absorption coefficient of the sample gas from the non-exponential behavior of the SCAR experimental decay curves is a challenging task.

The goal of the present invention is therefore to obtain from the SCAR data samples an accurate measurement of the gas concentration. This method includes a step to obtain a fit of the experimental data and from the fit parameters to obtain the desired concentration. Applicants have understood that in order to obtain an accurate estimate of the gas concentration from the experimental data obtained via a SCAR method, a physical understanding of the phenomena is needed to obtain meaningful parameters from the fit performed.

The Applicants have developed two different embodiments of the method of the invention to obtain accurate concentration determinations using the data obtained during a SCAR decay, or multiple SCAR decays. The first embodiment may give less accurate results, having under some conditions a slightly higher error in the final concentration determination, but can be performed “on-line”, that is, while the SCAR measurements are performed. One advantage of this embodiment is that the method of the invention can be performed during the measurements to determine whether there are evident errors in the procedure and whether a stop and a correction are needed. This embodiment can be performed on-line because it gives the result of the SCAR measurement one by one, that is, it fits each decay curve as soon as it is recorded for a given frequency. This method, of course, can be performed also off-line, saving all the measurements data for the different frequencies and evaluating the concentration of the gas using the method of the invention afterwards according to this embodiment. In case of relative values, that is when a ratio between a sample of gas having known concentration and a sample of gas having unknown concentration is desired, the two embodiments of the method of the invention are substantially equivalent in accuracy. The second method is implemented using a plurality of SCAR data relative to a number N of scar decays at different frequencies. Therefore, it cannot provide real time values, because as an input it needs all the data relative to all decays at different frequencies from which the concentration is to be determined. However, this second method has the highest accuracy, in particular it generally obtains the best value of the spectral area.

The apparatus used in any of the two main embodiments of the method of the invention to obtain SCAR decay data to be fitted includes a laser source emitting coherent electromagnetic radiation at a fixed wavelength/frequency. The laser source emits radiation with an intensity such as to work under conditions of saturation, as described below. The laser preferably emits infrared radiation. The infrared represents the optimal spectral range, since it has the most intense absorption bands of gaseous molecular species to be detected. Preferably, the laser emits radiation at a single wavelength, and it is a tunable laser (that can be tuned), i.e. whose wavelength/frequency is selectable within a certain range. This selection occurs in a known manner by means of modification of the wavelength which are known in the field. As is known, no laser is perfectly monochromatic, the laser of the invention is preferably a laser with a very narrow line-width, for example a CW laser wherein the line-width is less than the width of a cavity resonance mode.

In addition, the apparatus includes a resonant cavity with high finesse, e.g. comprising two reflecting mirrors arranged in such a way as to form a closed optical path for the electromagnetic radiation emitted by the laser source. Preferably, the round-trip length is finely adjustable (within a range of a few microns with sub-nm resolution) so as to put the cavity in resonance with the laser radiation, so as to couple the electromagnetic radiation to the cavity. The movement of the mirrors occurs through appropriate moving means known in the field. Preferably, the reflectivity of the mirrors of the cavity must be greater than 99.9%.

The resonant cavity is intended to contain the gas whose concentration is to be detected.

The arrangement of the laser source and of the cavity can be adjusted in such a way that the single wavelength radiation emitted by the source is coupled to the cavity, so that it is trapped therein, forming a resonant radiation.

If the laser is not monochrome, it is necessary that only one wavelength is substantially coupled to the cavity.

The type of laser source used in the present invention is for example described in the article written by Galli et al., Opt. Lett. 35, 3616 (2010). Other types of laser sources can be used, the intensity of the radiation emitted must however be such that the intensity of radiation I inside the cavity is much greater than the saturation intensity I_(S) of the molecular transition to be detected, i.e. I>>I_(S). Preferably, the saturation parameter, i.e. the ratio between the intensity of the radiation within the cavity and the saturation intensity of the molecular transition to be detected, or I/I_(S), must be >10.

The radiation/cavity combination enables an increase in the intensity of the radiation inside the cavity due to the resonance phenomenon. For a linear cavity, being λ the wavelength of the radiation coming from the laser and D the distance between the two mirrors, the cavity must satisfy the equation

$D = \frac{m\; \lambda}{2}$

with m integer.

The gas whose concentration is to be measured is located inside the cavity. Therefore, the cavity preferably includes an inlet and an outlet for the gas contained in dedicated containers. The cavity also includes, preferably, means for the thermal adjustment of the cavity itself so as to adjust the gas temperature, and pressure gauges to adjust the gas pressure inside the cavity.

The measurements are carried out by bringing the cavity in resonance with the electromagnetic radiation emitted by the laser, injecting it with the laser radiation to beyond a predetermined filling level. This coupling between the radiation and the cavity may be initially performed by the movement of one or more mirrors so as to vary the total length of the cavity and reach a cavity-laser resonance condition.

Preferably, as mentioned above, the laser emits radiation at different frequencies (or wavelengths) and radiation at various frequencies is introduced into the cavity, by scanning a certain range of frequencies. Such a scan is carried out across a molecular transition frequency of the gas to be studied and whose traces are to be detected, in order to measure the spectral area of the molecular transition, which is proportional to the gas concentration. Typically, the scan should be approximately from −5 to +5 HWHM Lorentzian widths or from −3 to +3 HWHM Gaussian widths, the largest of the two. For detection of radiocarbon is also preferable to extend the range beyond the interfering resonance of N₂O, then preferably the scanning should extend at least from −250 MHz to +400 MHz around the target radiocarbon resonance P(20).

Furthermore, the apparatus comprises a detector, in particular a photo detector, suited to detect the radiation leaking out of the cavity.

The electrical signal from the photodetector is sent to a processor. The processor processes such an input signal by a fit of the signal with a function that takes into account the effects of saturation, to calculate the spectral area of the linear absorption coefficient and thus the gas concentration present within the resonant cavity.

The detection is performed as follows: after the radiation has been introduced inside the cavity and a saturation condition has been obtained (i.e. in which the intensity of the radiation is much greater than the saturation intensity of the molecular transition of interest), the coupling between the cavity and the radiation emitted by the laser source is interrupted by any means known in the art and the decay of the radiation within the cavity itself is measured. For example, the coupling between the cavity and the laser source can be interrupted by switching off the laser source and/or by rapidly changing the emission frequency.

In order to take a ring-down measurement, the laser and the cavity are brought into a resonance condition, i.e. a condition in which the wavelength of the laser is a multiple of the optical path in the cavity. Therefore, all the measurements described below, are carried out in a condition where there is an initial resonance between the electromagnetic radiation emitted by the laser and the cavity. This resonance, hereinafter briefly called cavity-laser resonance, is “switched off” to perform the actual ring-down measurement.

Furthermore, to perform a measurement of the gas concentration, a further resonance condition is imposed, which is the resonance condition between the electromagnetic radiation emitted by the laser and the molecular transition of interest. This additional resonance will be briefly called laser-gas resonance hereinafter. The measurements are made around the laser-gas resonance, i.e. the frequency of the electromagnetic radiation emitted by the laser is “in a neighborhood” of the frequency of the transition of interest. Thus a range around the frequency of the molecular transition of interest ν₀ is determined, that is a ν_(min)≤ν₀≤ν_(max) range is selected and, starting preferably from one end of the range, the frequency of the laser is varied at pre-determined step Δν. For each frequency inside the selected range, a SCAR measurement is performed, that is, the cavity is brought to a resonance condition to this specific laser frequency selected, the intensity of the radiation inside the cavity is brought into saturation condition, the laser is switched off, and the decay of the intensity of the radiation is recorded.

Preferably, the operation above is repeated for all the frequencies in the selected range at each step Δν. The data so obtained are then fitted according to the first or the second embodiment of the invention.

The molecular concentration N_(g), in cm⁻³, of the target gas to be measured and positioned inside the cavity is directly related with the integrated area of the absorbance profile of the target absorption by:

$\begin{matrix} {N_{g} = {\frac{1}{{cL}_{s}(T)}{\int{{\alpha_{g}\left( {v - v_{o}} \right)}{dv}}}}} & \left( {{Eq}.\mspace{14mu} 3} \right) \end{matrix}$

with:

-   -   L_(s)(T) is the line-strength of the absorbent molecular         transition at temperature T, in cm/molecule     -   c is the vacuum velocity of light, in cm/s     -   ν₀ is the resonance frequency of the absorbent molecular         transition, in Hz     -   ν is the absorbed light frequency, in Hz     -   α_(g)(ν−ν₀) is the absorption coefficient, in cm⁻¹

The absorption coefficient α_(g)(ν−ν₀) is determined in SCAR spectroscopy by measuring the contribution to the cavity ring-down decay rate due to the absorption of the gas, γ_(g)(ν) as described in Eq. (2). Therefore, in order to determine the gas concentration, it is sufficient to calculate the cavity ring-down decay rate due to the absorption of the gas for a plurality of spectral frequencies around the resonant one, i.e. by measuring the spectral profile around the frequency of the target line.

Some results obtained with SCAR spectroscopy is reported in: G. Giusfredi, S. Bartalini, S. Borri, P. Cancio, I. Galli, D. Mazzotti and P. De Natale, “Saturated-absorption cavity ring-down spectroscopy”, Phys. Rev. Lett. 104, 110801 (2010) and in P. Cancio, I. Galli, S. Bartalini, G. Giusfredi, D. Mazzotti, and P. De Natale, “Saturated-Absorption Cavity Ring-Down (SCAR) for High-Sensitivity and High-Resolution Molecular Spectroscopy in the Mid IR”, in Cavity-Enhanced Spectroscopy and Sensing-Springer Series in Optical Sciences Volume 179, 2014, pp 143162, where however the method of the invention has not been disclosed.

In the following, the spectral line shape of emission of the transition examined in the target gas is considered. Three main factors generally defines an absorption line: central position of the line (e.g., the central frequency ν₀), the strength of the line, and shape factor (or line profile, g). Each line has a finite width (referred to as natural broadening of a spectral line). Several processes may result in an additional broadening of a spectral line of the molecules: such as collisions between molecules (referred to as the pressure broadening); due to the differences in the molecule thermal velocities (referred to as the Doppler broadening); and the combination of the above processes.

The experimental conditions, always with I>>I_(s), can be of two different types. They will be called in the following “homogenous condition” and “inhomogeneous condition”. In the inhomogeneous case, the line shape is dominated by a thermal Gaussian distribution and therefore the Lorentzian width due to collision and natural broadening is small, in other words when w_(G)>>w_(L), where w_(L) is the HWHM Lorentzian width and w_(G) is the HWHM Gaussian width. In a homogeneous case, the thermal Gaussian distribution dominates. Further, there is a third “intermediate” condition in which the two processes are both relevant.

The area normalized profile of the spectral line is called in the following:

g(ν−ν₀ ,w _(R))  (Eq. 4)

with

-   -   ν₀ defined as above     -   w_(R) is the half-width-half-maximum of the profile, in Hz

This can be written, in the three different experimental conditions above identified:

-   -   a) Inhomogeneous broadening (IB)>>homogeneous broadening (HB)         (i.e. w_(G)>>w_(L))

Gaussian profile of Gaussian width w_(G) (HWHM):

${g\left( {{v - v_{o}},w_{R}} \right)} = {{G\left( {{v - v_{o}},w_{G}} \right)} = {\frac{\sqrt{\ln \; 2}}{\sqrt{\pi}w_{G}}{\exp \left\lbrack {{- {\ln (2)}}\left( \frac{v - v_{o}}{w_{G}} \right)^{2}} \right\rbrack}}}$

-   -   b) IB<<HB (i.e. w_(G)<<w_(L))

Lorentzian profile of Lorentzian width w_(L) (HWHM):

${g\left( {{v - v_{o}},w_{R}} \right)} = {{L\left( {{v - v_{o}},w_{L}} \right)} = {\frac{1}{\pi \; w_{L}}\frac{1}{1 + \left( \frac{v - v_{o}}{w_{L}} \right)^{2}}}}$

-   -   c) IB≈HB (i.e. w_(G)≈w_(L)) (condition for maximum absorption,         normally)

Voigt profile of Voigt width w_(R) (HWHM), w_(R)={w_(L), w_(G)}: g(ν−ν_(o),w_(R))=V(ν−ν_(o),w_(G))

${g\left( {{v - v_{o}},w_{R}} \right)} = {{V\left( {{v - v_{o}},w_{R}} \right)} = {\frac{\sqrt{\ln \; 2}}{w_{G}w_{L}\sqrt{\pi^{3}}}{\int_{- \infty}^{\infty}{\frac{\exp \left\lbrack {{- {\ln (2)}}\left( {v^{\prime}/w_{G}} \right)^{2}} \right\rbrack}{1 + \left\lbrack {\left( {v - v^{\prime}} \right)/w_{L}} \right\rbrack^{2}}\ {dv}^{\prime}}}}}$

The peak normalized line profile is defined as:

$\begin{matrix} {{\overset{\_}{g}\left( {{v - v_{o}},w_{R}} \right)} = {\frac{1}{g_{o}}{g\left( {{v - v_{o}},w_{R}} \right)}}} & \left( {{Eq}.\mspace{14mu} 5} \right) \end{matrix}$

with

-   -   ν₀, w_(R) and g(ν−ν_(o),w_(R)), defined as above     -   g_(o)=g(0,w_(R)) is the peak normalization factor of the         profile, in Hz⁻¹.         -   a) IB>>HB (i.e. w_(G)>>w_(L)), g_(o)=G_(o),

${\overset{\_}{G}\left( {{v - v_{o}},w_{G}} \right)} = {\frac{1}{G_{o}}{G\left( {{v - v_{o}},w_{G}} \right)}}$

-   -   -   is the peak normalized Gaussian profile.         -   b) IB<<HB (i.e. w_(G)<<w_(L)), g_(o)=L_(o),

${\overset{\_}{L}\left( {{v - v_{o}},w_{L}} \right)} = {\frac{1}{L_{o}}{L\left( {{v - v_{o}},w_{L}} \right)}}$

-   -   -   is the peak normalized Lorentzian profile.         -   c) IB≈HB (i.e. w_(G)≈w_(L)), g_(o)=V_(o),

${\overset{\_}{V}\left( {{v - v_{o}},w_{R}} \right)} = {\frac{1}{V_{o}}{V\left( {{v - v_{o}},w_{R}} \right)}}$

-   -   -   is the peak normalized Voigt profile.

In order to satisfy the condition that the intensity of radiation emitted by the laser and then present in the cavity should be much higher that the saturation condition, the saturation condition may be calculated, if desired. The saturation intensity is given by, in Wm⁻²:

-   -   a) Homogeneous case (i.e. w_(G)<<w_(L)):

$\begin{matrix} {I_{s} = {\frac{c\; ɛ_{o}\hslash^{2}\gamma_{//}}{2\; \mu_{ba}^{2}}\frac{\gamma_{//}}{g_{o}}\frac{{2\; \gamma_{R}} + \gamma_{//}}{\gamma_{R} + \gamma_{//}}}} & \left( {{Eq}.\mspace{14mu} 6} \right) \end{matrix}$

-   -   -   with:             -   ε_(o) is the vacuum electric permittivity, in C²m⁻¹J⁻¹             -   ℏ is the Plank constant/(2π), in Js             -   c is the vacuum velocity of light, in m/s             -   g_(o) defined as above (g_(o)=L_(o), when the lineshape                 is essentially Lorentzian)             -   μ_(ba) is the electrical dipole of the absorbent                 transition, in Cm             -   γ_(R) is the relaxation rate due to collision that                 redistribute the population between the different                 rotational sublevels inside the same upper and lower                 vibrational states, in s⁻¹             -   γ_(//) is the population decay rate from the upper to                 the lower state, in s⁻¹                 -   In the case of γ_(//)<<γ_(R) and for upper                     rotational quantum number J′ equal to J or J−1, with                     J rotational number off the lower level (as trace                     gas detection of ¹⁴CO₂ by using the (00⁰1-00⁰0)                     P(20) transition)

$\begin{matrix} {{I_{s} \cong {\frac{c\; ɛ_{o}\hslash^{2}}{\mu_{ba}^{2}}\frac{\gamma_{//}}{g_{o}}}} = {\frac{c\; \hslash \; k_{o}^{3}}{\pi \; A_{ba}}\frac{\gamma_{//}}{g_{o}}}} & \left( {{Eq}.\mspace{14mu} 7} \right) \end{matrix}$

-   -   -   -   -   with:                 -   ε_(o), ℏ, c, μ_(ba), g_(o), γ_(//) defined as above                 -   k_(o) is the wave vector of the light at ν_(o), in                     m⁻¹                 -   A_(ba) is the Einstein A coefficient of the                     absorbent transition, in s⁻¹

    -   b) Inhomogeneous case (i.e. w_(G)>>w_(L)):

$\begin{matrix} {I_{s} = {\frac{c\; ɛ_{o}\hslash^{2}}{2\; \mu_{ba}^{2}}\frac{{2\; \gamma_{R}} + \gamma_{//}}{\gamma_{R} + \gamma_{//}}\gamma_{\bot}\gamma_{//}}} & \left( {{Eq}.\mspace{14mu} 8} \right) \end{matrix}$

-   -   -   with:             -   ε_(o), ℏ, c, μ_(ba), γ_(R) and γ_(//) defined as above                 -   In the case of γ_(//)<<γ_(R) and for J′ equal to J                     or J−1

$\begin{matrix} {I_{s} \cong {\frac{\hslash \; k_{o}^{3}}{\pi \; A_{ba}}\gamma_{R}\gamma_{//}}} & \left( {{Eq}.\mspace{14mu} 9} \right) \end{matrix}$

-   -   -   -   -   with:                 -   ℏ, c, k_(o), A_(ba), γ_(R) and γ_(//) defined as                     above.

    -   c) Intermediate case (i.e. w_(G)≈w_(L)):         -   The coefficient I_(s) is the saturation intensity, given by

$\begin{matrix} {I_{s} = {{\frac{c\; ɛ_{o}\hslash^{2}}{2\; \mu_{ba}^{2}}\frac{\gamma_{//}}{g_{o}}\frac{{2\; \gamma_{R}} + \gamma_{//}}{\gamma_{R} + \gamma_{//}}} \cong {\frac{c\; ɛ_{o}\hslash^{2}}{\mu_{ba}^{2}}\frac{\gamma_{//}}{g_{o}}}}} & \left( {{Eq}.\mspace{14mu} 10} \right) \end{matrix}$

-   -   -   where γ_(R) is the relaxation rate due to collision that             redistribute the population between the different rotational             sublevels inside the same upper and lower vibrational             states, while γ_(//) is the population decay rate from the             upper to the lower state. The last quasi-equality holds when             γ_(//)<<γ_(R).

Next, we have to consider that in a high-finesse cavity the gas interact with a laser beam that in the best condition has a TEM₀₀ profile with an intensity given by

I(ρ,t)=I ₀(t)e ^(−2(ρ/w)) ²   (Eq. 11)

Where ρ≡√{square root over (x²+y²)} is the radial coordinate and I₀(t)=I(ρ=0, t) is the peak intensity taking as axis z the axis of cavity, while w is the spot size radius of the laser beam, i.e. the radius for which the amplitude of the field is 1/e times that of the axis. Moreover, w is considered to be constant along axis z in the cavity. The power of the beam is

$\begin{matrix} {{P(t)} = {\frac{\pi \; w^{2}}{2}{I_{o}(t)}}} & \left( {{Eq}.\mspace{14mu} 12} \right) \end{matrix}$

The dependence on time t remembers that in a CRD event the intensity decays as consequence of the interruption of the laser beam that is sent into the cavity.

A saturation power is defined, in watt:

$\begin{matrix} {P_{s} = {\frac{\pi \; w^{2}}{2}I_{s}}} & \left( {{Eq}.\mspace{14mu} 13} \right) \end{matrix}$

where w is the spot size radius of the laser beam, i.e. the radius for which the amplitude of the field is 1/e times that of the axis, in m as well as a saturation parameter

$\begin{matrix} {Z = {\frac{I}{I_{s}} = \frac{P}{P_{s}}}} & \left( {{Eq}.\mspace{14mu} 14} \right) \end{matrix}$

where I and P intensity and power of absorbed light, respectively as defined in the above equations (11) and (12).

For developing the following formulas, two approximation may be done, according to the case:

-   -   1) adiabatic approximation: the relaxation rate of the molecular         population and coherence γ_(//), γ_(R) and γ_(⊥) are greater         than the power relaxation rate γ_(c) of the radiation inside de         cavity;     -   2) local approximation: due to molecular collisions, the         diffusion time across the laser beam is slower than the         relaxation time of the molecular population in the exited state.         In other words, the molecules have a negligible displacement         with respect to the beam size during the time interval necessary         to approximate their equilibrium condition with the radiation.

For the above adiabatic approximation the dependence on t is negligible in a round trip time of the radiation inside the cavity, therefore such dependence is omitted here.

Eq. (12) of the intracavity power decay in case of non-linear losses can be rewritten as:

P(t)=P(0)exp(−γ_(c) t)ƒ(t)  (Eq. 15)

with:

-   -   P(0) intracavity power at the beginning of the decay event, in W     -   γ_(c) defined in eq. (1)     -   ƒ(t) with ƒ(0)=1 is the non-linear faction which takes into         account the non-linear effects as saturation.

According to theory, still in a saturation condition, that is in an intracavity saturated absorption where P(0)>>P_(S) defined in Eq. (13) and I>>I_(s), the absorption coefficient in cm⁻¹ of Eq. (3) and the intracavity power of equation (15) follows the following equations in the three different experimental conditions named above a), b) and c).

-   -   a) homogeneous regime (w_(L)>>w_(G))         -   In this regime, both adiabatic and local approximation are             applied.

$\begin{matrix} {{\alpha_{g}\left( {{v - v_{o}},I} \right)} = {\alpha_{o}L_{o}\frac{\overset{\_}{L}\left( {{v - v_{o}},w_{R}} \right)}{1 + {\frac{I}{I_{s}}{\overset{\_}{L}\left( {{v - v_{o}},w_{R}} \right)}}}}} & \left( {{Eq}.\mspace{14mu} 16} \right) \end{matrix}$

-   -   -   with:             -   ν, ν_(o), w_(R), g_(o), g(ν−ν_(o),w_(R)), α_(o) and                 I_(s) defined as above             -   I is the instantaneous intracavity intensity, in W/m²

$\begin{matrix} {\frac{dP}{dt} = {{{- \gamma_{c}}P} - {{\gamma_{g}\left( {v - v_{o}} \right)}\frac{\ln \left\lbrack {1 + {\frac{P}{P_{s}}{\overset{\_}{g}\left( {{v - v_{o}},w_{R}} \right)}}} \right\rbrack}{\frac{P}{P_{s}}{\overset{\_}{g}\left( {{v - v_{o}},w_{R}} \right)}}P}}} & \left( {{Eq}.\mspace{14mu} 17} \right) \end{matrix}$

-   -   -   with:         -   γ_(g)(ν−ν_(o))=cα_(o)g(ν−ν_(o),w_(R))         -   ν, ν_(o), w_(R), g_(o), g(ν−ν_(o),w_(R)), g(ν−ν_(o),w_(R)),             γ_(c) and P_(s) defined as above         -   The resonance profile g both appears in γ_(g) and in the             fraction of Eq. (17). This fraction represents the             non-linear saturation term and its value tends to 1 as             P/P_(s)→0. In such case, integration of Eq. (17) produces a             truly exponential behavior, as for conventional CRD             spectroscopy. On the other hand, this fraction tends to 0 as             P/P_(s)→∞, which corresponds to the case of strong             saturation, when the gas becomes transparent to laser             radiation.         -   In SCAR experiments, the measured decay signal extends over             several decades in amplitude, with a S/N ratio mainly             limited by the resolution of the A/D converter. Such signal             must be fitted by the power decay function resulting from             integration of Eq. (17). To this purpose, it is more             convenient to factorize intracavity power into the two terms             of Eq. (15), obtaining:

$\begin{matrix} {\frac{df}{dt} = {{- {\gamma_{g}\left( {v - v_{o}} \right)}}\frac{\ln \left\lbrack {1 + {Z_{o}{\overset{\_}{g}\left( {{v - v_{o}},w_{R}} \right)}{\exp \left( {{- \gamma_{c}}t} \right)}f}} \right\rbrack}{Z_{o}{\overset{\_}{g}\left( {{v - v_{o}},w_{R}} \right)}{\exp \left( {{- \gamma_{c}}t} \right)}}}} & \left( {{Eq}.\mspace{14mu} 18} \right) \end{matrix}$

-   -   -   with:             -   ν, ν_(o), w_(R), g_(o), g(ν−ν_(o),w_(R)), γ_(c) and                 γ_(g)(ν−ν_(o)) defined as above

$Z_{o} = \frac{P(0)}{P_{s}}$

-   -   -   -   is the saturation parameter at the beginning of the                 decay for the saturation power at the center of the                 absorbent transition.

        -   The new function ƒ takes into account deviation from a             purely exponential decay, at rate γ_(c), of the experimental             SCAR signal due to the saturated gas absorption. The much             smaller dynamic range of the ƒ function allows to determine             the gas absorption decay rate γ_(g) with a lower             uncertainty.

    -   b) inhomogeneous regime (w_(G)>>w_(L))

$\begin{matrix} {{\alpha_{g}\left( {{v - v_{o}},I} \right)} = {\alpha_{o}\frac{G\left( {{v - v_{o}},w_{G}} \right)}{\sqrt{1 + \frac{I}{I_{s}}}}}} & \left( {{Eq}.\mspace{14mu} 19} \right) \end{matrix}$

-   -   -   with:         -   ν, ν_(o), w_(G), G(ν−ν_(o),w_(G)), α_(o), I and I_(s)             defined as above

    -   Like the homogeneous case, the interaction of the molecules with         a TEM₀₀ mode of the radiation inside a cavity have to be         considered. Therefore also for the following expressions, the         validity of the above mentioned adiabatic approximation is         assumed. Instead two further cases are considered:         -   the density of the gas is still large enough to allow the             above mentioned local approximation         -   the gas density is low enough to let the molecules cross the             laser beam without any collision.         -   First case: still diffusive gas (collisions are present             during gas-laser interaction, w_(L)<w_(G))         -   The local approximation is still valid.

$\begin{matrix} {\frac{dP}{dt} = {{{- \gamma_{c}}P} - {{\gamma_{g}\left( {v - v_{o}} \right)}\frac{2}{\sqrt{1 + \frac{P}{P_{s}}} + 1}P}}} & \left( {{Eq}.\mspace{14mu} 20} \right) \end{matrix}$

-   -   -   with:             -   P_(s) and γ_(c) defined as above             -   γ_(g)(ν−ν₀)=cα_(o)G(ν−ν_(o),w_(G))         -   or in terms of Eq. 15

$\begin{matrix} {\frac{df}{dt} = {{- {\gamma_{g}\left( {v - v_{o}} \right)}}\frac{2}{\sqrt{1 + {Z_{o}{\exp \left( {{- \gamma_{c}}t} \right)}}} + 1}f}} & \left( {{Eq}.\mspace{14mu} 21} \right) \end{matrix}$

-   -   -   with:             -   Z₀, γ_(c) and γ_(g)(ν−ν_(o)) defined as above         -   Second case: non-diffusive gas (collisions are not present             during gas-laser interaction, w_(L)<<w_(G))         -   In the second case, when the pressure is very low, so that             the molecules pass through the beam without collisions, the             local approximation fails. In this case we can roughly             approximate the interaction assuming an “average” intensity             along the path, say x, of the individual molecules in the             term I/I_(S) at the denominator. The integration of the             starting equation for P becomes not analytic, however it is             well approximated by

$\begin{matrix} {\frac{dP}{dt} = {{{- \gamma_{c}}P} - {{\gamma_{g}\left( {v - v_{o}} \right)}\frac{1.4256}{\sqrt{1 + \frac{P}{2P_{s}}} + 0.4256}P}}} & \left( {{Eq}.\mspace{14mu} 22} \right) \end{matrix}$

-   -   -   with:             -   P_(s), γ_(c) and γ_(g)(ν−ν_(o)) defined as above         -   or in terms of Eq. 15:

$\begin{matrix} {\frac{df}{dt} = {{- {\gamma_{g}\left( {v - v_{o}} \right)}}\frac{1.4256}{\sqrt{1 + {0.5Z_{o}{\exp \left( {{- \gamma_{c}}t} \right)}f}} + 0.4256}f}} & \left( {{Eq}.\mspace{14mu} 23} \right) \end{matrix}$

-   -   -   with:             -   Z_(o), γ_(c) and γ_(g)(ν−ν_(o)) defined as above.

For each frequency of the electromagnetic radiation emitted by the laser, a SCAR measurement is made, i.e. the decay signal of the radiation leaking out of the cavity is measured. As a result of each SCAR experiment, for each wavelength, a saturated cavity ring-down signal for each spectral frequency, namely a S_(exp)(t; δν_(i)), is obtained. It is a result of the detection of the cavity-transmitted power due to cavity ring-down decay. The behavior of S_(exp)(t; δν_(i)) is described by the integration of a differential equation which takes into account the spectroscopic dominant processes of the targeted absorption transition depending on the gas temperature and pressure conditions, as detailed below. Imposing a specific behavior to this function which is the fit of the SCAR measurements, allows to obtain a proper determination of the gas concentration.

The behaviour over time t of the SCAR decay function S, known the selected frequency of the laser beam, the resonance frequency, the thermodynamic conditions of pressure and temperature of the target gas, depends on the following parameters:

S(t)=function of t;γ _(c),γ_(g) ,Z ₀ ,g,B,A _(d) =CP(0)

It may be factorized, for practical reasons in the fit, in U

S(t)=B+CP(0)e ^(−γ) ^(c) ^(t)ƒ(t;γ _(c),γ_(g) ,Z _(o) ,g )  (Eq. 24)

with

-   -   B is the detection background, in U     -   C is the detection conversion factor U/W     -   P(0), γ_(c), γ_(g) and g defined as above

$Z_{o} = \frac{{CP}(0)}{{CP}_{s}}$

-   -   is the saturation parameter at the begin of the decay event and         at the frequency of the absorbent transition, in U

In the following, we denote as laser frequency detuning with respect to the frequency of the absorbent transition, in Hz, as

δν_(i)=(ν_(i)−ν_(o))  (Eq. 25)

The experimental situation is a SCAR experiment in a gas mixture that can have l=0 . . . n absorptions in the spectral range measured, that is inside the cavity there is not only the target gas the concentration of which is to be measured, but also other gases with their transition lines, which are in this case a “noise” which distorts the decay function to be measured. This is the more general case, typically found in the trace gas detection experiments, where the sample gas is a mixture of gases that can absorb resonant light inside of cavity. The l=0 . . . n absorption transitions in the detected spectral range {δν_(min) . . . δν_(max)} are considered, with:

-   -   l=0 is the transition corresponding to the target gas, with the         corresponding N₁ density to be determined;     -   l=1 . . . n are the transitions of other gases, which interfere         with the target one, possibly introducing an error in the         measured absorption area.

For each laser detuning, that is for each frequency ν_(i) selected by the laser in a neighbourhood of the resonance frequency of interest ν₀ where δν_(i)=ν_(i)−ν₀ the SCAR decay can be fitted by the following curve, which is a generalization of the equation (24):

S(t;p(δν_(i)))=B(δν_(i))+A _(d)(δν_(i))e ^(−γ) ^(c) ^((δν) ^(i) ^()t)ƒ(t;γ _(c)(δν_(i)),A _(d)(δν_(i)),γ_(g)(δν_(i)),p _(l)(ν_(i)))  (Eq. 26)

with:

-   -   p(δν_(i)) is the set of parameters that define the SCAR decay         behavior:

p(δν_(i))={B(δν_(i)),A _(d)(δν_(i)),γ_(c)(δν_(i)),γ_(g)(δν_(i)),p _(l)(ν_(i)))},l=0 . . . n  (Eq. 27)

-   -   B(δν_(i)) is the detection background, in U     -   A_(d)(δν_(i)) is the amplitude of the decay signal at the begin         of the decay event, in U     -   γ_(c)(δν_(i)) and γ_(g)(δν_(i)) defined as above     -   p_(l)(ν_(i)) are n-set of parameters that defines each absorbent         transition

p _(l)(ν_(i))={(Z _(1U))_(l) ,g _(l)(ν_(i)−ν_(l) ,w _(Rl))},l=0 . . . n  (Eq. 28)

-   -   And the following factorization for Z_(o) has been used:

$\begin{matrix} {Z_{o} = {\frac{{CP}(0)}{{CP}_{s}} = {A_{d}Z_{1U}}}} & \left( {{Eq}.\mspace{14mu} 29} \right) \end{matrix}$

-   -   -   C is the detection conversion factor U/W         -   A_(d)=CP(0) is the amplitude of the decay signal at the             begin of the decay event, in U

$Z_{1U} = \frac{1}{{CP}_{s}}$

-   -   -   is the saturation parameter for an amplitude A_(d)=1 U, in U         -   ƒ the non-linear function that follow one of the rate             equations and depending of the gas conditions is described             by:             -   homogeneous regime (w_(L)≥w_(G))

$\begin{matrix} {\frac{df}{dt} = {{- {\gamma_{g}\left( {\delta \; v_{i}} \right)}}{\sum\limits_{l = 0}^{n}\frac{\ln \left\lbrack {1 + {{A_{d}\left( {\delta \; v_{i}} \right)}\left( Z_{1U} \right)_{l}{{\overset{\_}{g}}_{l}\left( {{v_{i} - v_{l}},w_{Rl}} \right)}e^{{- {\gamma_{c}{({\delta \; v_{i}})}}}t}f}} \right\rbrack}{{A_{d}\left( {\delta \; v_{i}} \right)}\left( Z_{1U} \right)_{l}{{\overset{\_}{g}}_{l}\left( {{v_{i} - v_{l}},w_{Rl}} \right)}}}}} & \left( {{Eq}.\mspace{14mu} 30} \right) \end{matrix}$

-   -   -   -   which is a generalization of Eq. 18 for l=0 . . . n                 transitions with (Z₀)_(l)=A_(d)(Z_(1U))_(l) g _(l)             -   inhomogeneous regime             -   still diffusive gas (w_(L)<w_(G))

$\begin{matrix} {\frac{df}{dt} = {{- {\gamma_{g}\left( {\delta \; v_{i}} \right)}}{\sum\limits_{l = 0}^{n}{\frac{2}{\sqrt{1 + {{A_{d}\left( {\delta \; v_{i}} \right)}\left( Z_{1U} \right)_{l}e^{{- {\gamma_{c}{({\delta \; v_{i}})}}}t}f}} + 1}f}}}} & \left( {{Eq}.\mspace{14mu} 31} \right) \end{matrix}$

-   -   -   -   which is a generalization of Eq. 21 for l=0 . . . n                 transitions with (Z₀)_(l)=A_(d)(Z_(1U))_(l)             -   non-diffusive gas (w_(L)<<w_(G))

$\begin{matrix} {\frac{df}{dt} = {{- {\gamma_{g}\left( {\delta \; v_{i}} \right)}}{\sum\limits_{l = 0}^{n}{\frac{1.4256}{\sqrt{1 + {0.5{A_{d}\left( {\delta \; v_{i}} \right)}\left( Z_{1U} \right)_{l}e^{{- {\gamma_{c}{({\delta \; v_{i}})}}}t}f}} + 0.4256}f}}}} & \left( {{Eq}.\mspace{14mu} 32} \right) \end{matrix}$

-   -   -   -   which is a generalization of Eq. 23 for l=0 . . . n                 transitions with (Z₀)_(l)=A_(d)(Z_(1U))_(l)

In the case of only the target transition is detected, Eqs. 30-32 are applied for l=0.

The above equation is in “U” unit, that is 1 U=1 physical unit, either volt or microampere, depending on the kind of front-end electronics used for the detection.

Now for the finding the function S, that is, for the finding of the values of the parameters on which the function S depends, in order to find the correct S based on the experimental data, the method of the invention is used.

The fitting procedure uses a routine that includes Eq. (26) and any of equations Eq. (30-32), depending on the physical situation, thus reproducing the decay signals to be compared with the experimental ones. This routine is called by a fitting program to estimate the values of the parameters: B, A_(d), γ_(c), γ_(g), Z_(1U) and g.

In both cases γ_(g) contains implicitly the dependence on the line-shape g, however there is an important difference in the expression of ƒ. In the homogeneous case the initial saturation parameter Z_(o)=A_(d)Z_(1U) g depends on the line-shape, while in the inhomogeneous case, with Z_(o)=A_(d)Z_(1U), it does not.

Therefore, in the inhomogeneous case, the experimental line profile γ_(g) can be obtained by fitting the each decay signal S(t) independently at a number of different laser frequencies.

Instead, unfortunately, for the homogeneous case, the line profile parameters are difficult to fit or even to impose as fixed values in the procedure. Indeed, interferences from other molecular resonances, as well as residual instrumental non-linearity, strongly affect the determination of the line profile, when the absorption is very weak, as in the case of trace gas detection. As a consequence, artifacts in the fitted line shape are observed, and further approximations must be done.

In the following, the method of the invention can be applied to any of the physical conditions, that is, both in the homogeneous and in the inhomogeneous case, however it is particularly advantageous in homogeneous case, which in particular is useful for radiocarbon dioxide detection, and their effects on concentrations retrieval.

According to a first aspect, the invention relates to a method of ring-down spectroscopy in saturated-absorption condition, for measuring a first concentration of a gas through a measurement of the spectrum of a molecular transition of said gas, the method comprising the steps of:

-   -   inserting said gas whose first concentration is to be measured         in a resonant cavity comprising two or more reflecting mirrors         arranged so as to form a closed optical path for an         electromagnetic radiation emitted by a laser source;     -   tuning the frequency of said electromagnetic radiation emitted         by said laser source so as to fix it to a value ν_(i) within a         range of frequencies [ν_(min), ν_(max)] including the resonance         frequency of said molecular transition ν₀;     -   fixing the intensity of said electromagnetic radiation in the         cavity at a value much greater than the saturation intensity         I_(s) of the molecular transition to be detected;     -   irradiating said gas by means of said electromagnetic radiation         beam emitted by said laser source having said fixed frequency         ν_(i) and intensity in said resonant cavity;     -   coupling said electromagnetic radiation to said cavity so as to         obtain a laser-cavity resonance condition;     -   changing the frequency of the electromagnetic radiation emitted         by the laser so as to switch off the laser-cavity resonance;     -   detecting an electromagnetic radiation beam in output from said         cavity after the laser-cavity resonance has been switched off;     -   recording a plurality of data representative of said output         which has the form of a decay signal;     -   considering a fitting curve S(t, δν_(i)) for the recorded decay         which depends on the following parameters:         -   B(δν_(i)) is the detection background, with             (δν_(i)=ν_(i)−ν₀)         -   A_(d)(δν_(i)) is the amplitude of the decay signal at the             beginning of the decay event,         -   γ_(c)(δν_(i)) is the cavity decay rate due to non-resonant             and non-saturable losses (empty cavity decay rate);         -   γ_(g) (δν_(i)) is contribution of the targeted molecular             transition to the decay signal;         -   g is the peak normalized line profile g(ν−ν_(o),w_(R))             centered at the molecular resonance frequency ν₀ and w_(R)             is the HWHM width of the resonance, and w_(R)=w_(L) for a             Lorentzian shape, w_(R)=w_(G) for a Gaussian shape,             w_(R)={w_(L), w_(G)} for a Voigt shape;

$Z_{o} = {\frac{{CP}(0)}{{CP}_{s}} = {A_{d}Z_{1U}}}$

-   -   -   is the saturation parameter at the beginning of the decay             event and at the frequency of the targeted molecular             transition ν₀, P(0) is the intracavity power at the             beginning of the decay signal; and

$P_{s} = {\frac{\pi \; w^{2}}{2}I_{s}}$

-   -   -   is the saturation power, where w is the spot size radius of             the laser beam, i.e. the radius for which the amplitude of             the field is 1/e times that of the axis and I_(s) the             saturation intensity;

    -   replacing Z_(1U) g(ν−ν₀,w_(R)) in the function S(t, δν_(i)) with         a constant value Z_(1ueff)=constant of a predetermined value;

    -   fitting said recorded data with a function S^(repl)(t, δν_(i))         in which Z_(1Veff)=constant replaces Z_(1U) g(ν−ν_(o),w_(R)) in         the fitting function S(t, δν_(i)).

According to this first aspect, the SCAR measurement is performed first of all selecting a laser frequency in a neighborhood of the transition frequency, that is of the absorption frequency of the target gas of interest. The absorption frequency is called as above ν₀ and the selected frequency ν_(i). The difference ν_(i)−ν₀ is called δν_(i) and it is called a “detuning” from resonance, meaning that for this selected frequency there is not a perfect laser-gas resonant condition. It is to be understood that a selected frequency of the laser radiation can almost precisely coincide with the resonant frequency, that is, it is also encompassed the case ν_(i)=ν₀.

Typically, the scan is preferably approximately from −5 to +5 HWHM Lorentzian widths or from −3 to +3 HWHM Gaussian widths, the largest of the two. For detection of radiocarbon is also preferable to extend the range beyond the interfering resonance of N₂O, then preferably the scanning should extend at least from −250 MHz to +400 MHz around the target radiocarbon resonance P(20).

The target gas inside the cavity, which is brought to a cavity-laser resonance condition, is in a saturation condition, that is, the intensity of the radiation inside the cavity is such that I>>I_(S). The value of I_(s) to be considered depending on the experimental conditions of the gas, that is, it is calculated according to one of equations (6) or (8). Pressure and temperature in the cavity are thus known parameters.

In order to measure the SCAR decay, that is, the “cavity ring down decay in saturated condition”, the laser-cavity resonance is switched off, for example switching off the laser or changing the frequency of the emitted radiation. Once switched off the laser-cavity resonance, the radiation exiting the cavity, which represents a decay, is detected and corresponding data are collected.

Indeed, as above shown describing the apparatus of the invention, from the experimental point of view, SCAR does not differ from other conventional CRD experiments. Laser light is coupled to a high-finesse Fabry-Perot cavity up to a threshold level, and then it is quickly switched off a cavity resonance. Transmitted ring-down light signal S_(exp)(t,ν) is detected and the decay rate is measured. If a molecule inside of the cavity absorbs coupled light, the loss rate due to such absorption is measured as an increase of the empty cavity loss rate. What is different in SCAR is that saturation effects of the molecular absorption induce a deviation of the ring-down signal from the perfectly exponential behavior, as expected for linear intracavity losses. Indeed, when the ring-down signal that contains saturated molecular absorption is fitted to an exponential function, the residuals differ strongly from a flat behavior. Instead, a fit to a decay function which takes into account non-linear absorption effects, gives a flat residual plot, hence explaining the physical effect of the saturated molecular transition.

Thus for each selected frequency of the laser in a neighborhood of the resonance frequency of the transition of interest in the target gas, a plurality of data is obtained. Preferably the laser emits in an infrared range.

The detected signal can be expressed in some physical unity “U”, such that i.e. 1U=1 volt or 1U=1 microampere. In the following 1U has the meaning of 1 physical unity. Therefore the signal is given by:

For each frequency of the electromagnetic radiation emitted by the laser, a SCAR measurement is made, i.e. the decay signal of the radiation leaking out of the cavity is measured. As a result of each SCAR experiment, a saturated cavity ring-down signal for each spectral frequency, namely a S_(exp), is obtained. It is a result of the detection of the cavity-transmitted power due to cavity ring-down decay. The behavior of S_(exp) is described by the integration of a differential equation which takes into account the spectroscopic dominant processes of the targeted absorption transition depending on the gas temperature and pressure conditions, as detailed below. Imposing a specific behavior to this function which is the fit of the SCAR measurements, allows to obtain a proper determination of the gas concentration.

Each detected SCAR decay data S_(exp)(t; δν_(i)) at each δν_(i) detuning is to be fitted with a function S(t; p(δν_(i))) which is in a preferred embodiment parameterized as (see again equation (26)):

S(t;p(δν_(i)))=B(δν_(i))+A _(d)(δν_(i))e ^(−γ) ^(c) ^((δν) ^(i) ^()t)ƒ(t;γ _(c)(δν_(i)),A _(d)(δν_(i)),γ_(g)(δν_(i)),Z _(1Ueff))

Which in principle depends on the following parameters p(δν_(i)), which is set of parameters that define the SCAR decay behavior:

p(δν_(i))={B(δν_(i)),A _(d)(δν_(i)),γ_(c)(δν_(i)),γ_(g)(δν_(i)),Z _(1Ueff)}

ƒ is the non-linear function that follows one of the equations (30-32) depending of the target gas experimental conditions (and imputed in the equation).

The fact that the behaviour of the function is the one of equation (26) is determined by the physics of the processes involved. The behavior of the ƒ function is to be determined by using the physics theory that takes into account the saturated absorption phenomena inside the cavity. Within this theory, the interaction between laser light and molecules for the targeted absorption transition must be analyzed in the frame of the total pressure and temperature conditions of the gas inside the cavity and of the laser power and spatial profile.

Therefore, in order to obtain the concentration of the gas of interest, the SCAR decay data are to be fitted by the function which has the form of equation (26) and at the same time satisfy any of the differential equations of Eq. (30)-(32). That is to say that the following considerations regarding the method of the invention are independent on the specific differential equation that the function ƒ has to satisfy. The differential equation to be satisfied depends on the specific physic phenomenon taking place during the experiment which depends on the experimental conditions.

However, Applicants have realized that performing a fit of eq. (26) to obtain a function S, an extremely complex numerical calculation is necessary. Substantially, the presence of so many parameters which are frequency-dependent does not allow a convergence towards a reasonable result. Therefore, Applicants have developed a method that allows real-time determination of γ_(g) for each spectral frequency.

This method tries to fix the problem of the proper frequency dependence of the saturation parameter, which can be strongly correlated with γ_(g) parameter. As a consequence, strong distortion of the wings of the spectral profile can be observed, and hence a wrong determination of the integrated absorption area. To solve this problem, the method of the invention defines a frequency independent saturation parameter, the effective saturation parameter Z_(1Ueff), that takes the place of the product Z_(1U) g(ν−ν_(o),w_(R)). This Z_(1Ueff) is to be determined in each experimental spectral situation, and it is a fixed parameter in the fitting routine for all spectral frequencies.

In other words, g(ν_(i)−ν_(o),w_(R)) is not any more a parameter of the fitting routine (more precisely, the line center and linewidth that characterizes this function are fixed).

Applicants have found that with this approximation, that is, keeping the value of the g(ν_(i)−ν_(o),w_(R)) constant to a given value initialized before the fitting routine, the fit generally converges. Although an approximation is made, the results on the concentration of gas, also at very low concentration is rather accurate, with a much higher accuracy than in standard ring-down spectroscopy.

Preferably, at the end of the fitting routine, the following parameters are obtained:

-   -   B(δν_(i)) is the detection background, with (δν_(i)=ν_(i)−ν₀),     -   A_(d)(δν_(i)) is the amplitude of the decay signal at the         beginning of the decay event,     -   γ_(c)(δν_(i)) is the cavity decay rate due to non-resonant and         non-saturable losses (empty-cavity decay rate),     -   γ_(g) (δν_(i)) is the contribution of the targeted molecular         transition to the decay signal.

In the fitting routine, all the other parameters on which the decay signal depends are of course free and from the fitting their value is determined. An important parameter is the last one which is connected to the gas concentration, as shown in equation (1).

Preferably, the method of the invention includes the step of parametrizing said S_(repl)(t, ν_(i)) as:

S(t;p(δν_(i)))=B(δν_(i))+A _(d)(δν_(i))e ^(−γ) ^(c) ^((δν) ^(i) ^()t)ƒ(t;γ _(c)(δν_(i)),A _(d)(δν_(i)),γ_(g)(δν_(i)),Z _(1Ueff))

where p(δν_(i)) is the set of parameters free during the fit

p(δν_(i))={B(δν_(i)),A _(d)(δν_(i)),γ_(c)(δν_(i)),γ_(g)(δν_(i))}

-   -   B(δν_(i)) is the detection background,     -   γ_(g)(δν_(i)) is the contribution of the targeted molecular         transition to the decay signal,     -   A_(d)(δν_(i)) is the amplitude of the decay signal,     -   γ_(c)(δν_(i)) is the cavity decay rate due to non-resonant and         non-saturable losses,     -   while the following are determined before the fitting:         -   Z_(1Ueff) is the effective saturation parameter, fixed             during the fit and equal to a constant;         -   ƒ the non-linear function that follows one of the below rate             equations depending of the gas conditions:             -   homogeneous regime (w_(L)≥w_(G))

$\frac{df}{dt} = {{- {\gamma_{g}\left( {\delta \; v_{i}} \right)}}\frac{\ln \left\lbrack {1 + {{A_{d}\left( {\delta \; v_{i}} \right)}Z_{1{Ueff}}e^{{- {\gamma_{c}{({\delta \; v_{i}})}}}t}f}} \right\rbrack}{{A_{d}\left( {\delta \; v_{i}} \right)}Z_{1{Ueff}}}}$

-   -   -   -   inhomogeneous regime                 -   still diffusive gas (w_(L)<w_(G))

$\frac{df}{dt} = {{- {\gamma_{g}\left( {\delta \; v_{i}} \right)}}\frac{2}{\sqrt{1 + {{A_{d}\left( {\delta \; v_{i}} \right)}Z_{1{Ueff}}e^{{- {\gamma_{c}{({\delta \; v_{i}})}}}t}f}} + 1}f}$

-   -   -   -   -   non-diffusive gas (w_(L)<<w_(G))

$\frac{df}{dt} = {{- {\gamma_{g}\left( {\delta \; v_{i}} \right)}}\frac{1.4256}{\sqrt{1 + {0.5{A_{d}\left( {\delta \; v_{i}} \right)}Z_{1{Ueff}}e^{{- {\gamma_{c}{({\delta \; v_{i}})}}}t}f}} + 0.4256}f}$

The parametrization of the fitting curve S has been chosen by the Applicants due to the physical constraints that determine the shape of the decay signal, and thus which should also determine the shape of the fitting curve, in this case the behavior of the ƒ function (ƒ function which is chosen among the three above listed depending on the experimental conditions). Selecting such a parametrization allows an easier convergence during the fitting procedure. The ƒ functions of Eqs. (30)-(32) have been rewritten in function of the fixed parameter Z_(1Ueff).

In FIG. 1, the flowchart of this fitting procedure is shown. This method tries to fix the problem of the proper frequency dependence of the saturation parameter (Z_(1U)), which can be strongly correlated with γ_(g) parameter. As a consequence, strong distortion of the wings of the spectral profile can be observed, and hence a wrong determination of the integrated absorption area. To solve this problem, ESPS defines a frequency independent saturation parameter, the effective saturation parameter Z_(1Ueff), equal for all l transitions, which must be determined in each experimental gas conditions following the calibration routine described below (sec. 2.1), and it is a fixed parameter in the fitting routine for all laser detunings.

Preferably, the peak normalized line profile g(ν−ν_(o),w_(R)) centered at the molecular resonance frequency ν₀ has either a Lorentzian shape, a Gaussian shape, or a Voigt shape depending on the experimental conditions.

The line profile function shapes include Lorentzian, Gaussian and Voigt functions, whose parameters are the line position, maximum height and half-width. Actual line shapes are mainly determined by Doppler, collision and proximity broadening. For each system the half-width of the shape function varies with temperature and pressure (or concentration).

The peak normalized line profile is therefore either a Lorentzian, Gaussian or Voigt function depending on the experimental conditions as better clarified below. However, for the present invention, the line profile can be any line profile, because the teaching of the Invention can be applied to any of those. In order to minimize correlation effects, this parameter is fixed in the fitting routines (more correctly, the spectral parameters of line center and linewidth that characterized the profile are fixed). The initialization value for this parameter could be determined by considering the broadening mechanisms on the targeted transition. Alternatively, an experimental determination is performed.

Preferably, said fitting is a least squares fitting. More preferably, said least square fitting uses the Levenberg-Marquardt (L-M) algorithm.

This method considers the set of parameters p^(fr)(δν_(i)) to be fitted (i.e. free in the least square fitting method, such as in the L-M algorithm) and Z_(1Ueff) fixed during fit, i.e.

p(δν_(i))={p ^(fr)(δν_(i)),Z _(1Ueff) }={B ^(fr)(δν_(i)),A _(d) ^(fr)(δν_(i)),γ_(c) ^(fr)(δν_(i)),γ_(g) ^(fr)(δν_(i)),Z _(1Ueff)}  (Eq.33)

This are the parameters on which the SCAR decay signals depend, with the exception of Z_(1Ueff) which has been set equal to a known constant n.

The L-M algorithm performs this fit by minimizing χ²(p(δν_(i))) in the t variable in the space of the free parameters p^(fr)(δν_(i)) (see FIGS. 2 and 3):

$\begin{matrix} {{\chi^{2}\left( {p\left( {\delta \; v_{i}} \right)} \right)} = {\sum\limits_{t_{k} = 0}^{t_{\max}}\frac{\left( {{S_{\exp}\left( {t_{k};{\delta \; v_{i}}} \right)} - {S\left( {t_{k};{p\left( {\delta \; v_{i}} \right)}} \right)}} \right)^{2}}{\sigma_{\exp}^{2}\left( {t_{k};{\delta \; v_{i}}} \right)}}} & \left( {{Eq}.\mspace{14mu} 34} \right) \end{matrix}$

where σ_(exp)(t; δν_(i)) is the uncertainty of the detected decay S_(exp)(t; δν_(i)).

One of the possible routines is the following (see FIG. 3):

-   -   1) Consider the detected SCAR curve {S_(exp)(t;         δν_(i)),σ_(exp)(t; δν_(i))}_(j) for a sample that contains the         target gas at the thermodynamic conditions (Pressure,         Temperature)     -   2) Initialization procedure:         -   Set the initial free parameters to:

p _(j) ^(ini)(δν_(i))={B ^(ini)(δν_(i)),A _(d) ^(ini)(δν_(i)),γ_(c) ^(ini)(δν_(i)),γ_(g) ^(ini)(δν_(i))}_(j) ^(ini)

-   -   -   Set the Z_(1Ueff) parameter to a pre-defined calibrated             value, and keep it fixed during all fit procedure.         -   Other initializations: L-M limit convergence values,             λ_(max),δχ_(max) ²,χ_(min) ² and maximum number of             interactions N_(itr)

    -   3) Minimization loop on p_(j) ^(fr)(δν_(i))

Try as initial guess of the free parameters p_(j) ^(try)(δν_(i))=p_(j) ^(ini)(δν_(i)) and uses the L-M routine to obtain a fitted set of parameters p_(j) ^(fit)(δν_(i)).

We follow the procedure depicted in the flowchart of FIGS. 2 and 3.

-   -   a) p(δν_(i))={p_(j) ^(try)(δν_(i)),Z_(1Ueff)}. Calculate the         χ²(p(δν_(i))) (Eq. 34), gradient β(p(δν_(i))) and the curvature         matrix, |α|(p(δν_(i))). Here, the temporal interval of the decay         under analysis, [0,t_(max)] is divided in t_(k) points at         intervals of Δt. For each t_(k), the routine calculates:         -   The numerical integration of the differential equation for             ƒ(t_(k); p(δν_(i))) (Eqs. 30-32 depending of the case             applies). We use a fourth-order Runge-Kutta (RK4) routine             for such numerical integration but other integration             routines can be used.         -   The expected S(t_(k); p(δν_(i))) values by using the Eq.             (26).         -   By considering S_(exp)(t_(k); δν_(i)) and σ_(exp)(t_(k);             δν_(i)) as input values, the summations of the             χ²(p(δν_(i))), β(p(δν_(i))) and |α|(p(δν_(i))) are updated.         -   The above three steps are repeated until t_(k)=t_(max).     -   b) Solve the linear Eqs. of the L-M algorithm for the present         problem, and calculate δp(δν_(i)) and χ²(p(δν_(i))+δp(δν_(i)))     -   c) Follow the L-M convergence criteria on λ and/or χ² as         described in the L-M method. If not converging come back to         step a) with an update value of λ and/or p_(j)         ^(try)(δν_(i))=p(δν_(i))+δp(δν_(i)).     -   d) At convergence, a fitted set of parameters equal to the last         guess is determined

p _(j) ^(fit)(δν_(i))=p _(j) ^(try)(δν_(i))

-   -   4) As a result, the absorption decay rate parameters is         determined: |γ_(g)(δν_(i))|_(j=1) ^(d*m) with uncertainties

|σ_(γ) _(g) (δν_(i))|_(j=1) ^(d*m)

Preferably, the method of the invention includes:

-   -   4) As a result, the absorption decay rate parameters is         determined: |γ_(g)(δν_(i))|_(j=1) ^(d*m) with uncertainties

|σ_(γ) _(g) (δν_(i))|_(j=1) ^(d*m)

Preferably, the method of the invention includes:

-   -   Changing the frequency of the electromagnetic radiation emitted         by the laser to a frequency ν_(i+d) where ν_(i+d) belongs to         [ν_(min), ν_(max)] and repeating the steps of         -   fixing the intensity of said electromagnetic radiation in             the cavity at a value much greater than the saturation             intensity I_(s) of the molecular transition to be detected;         -   irradiating said gas by means of said electromagnetic             radiation beam emitted by said laser source having said             fixed frequency ν_(i+d) and intensity in said resonant             cavity;         -   coupling said electromagnetic radiation to said cavity so as             to obtain a laser-cavity resonance condition;         -   changing the frequency of the electromagnetic radiation             emitted by the laser so as to switch off the laser-cavity             resonance;         -   detecting an electromagnetic radiation beam in output from             said cavity after the laser-cavity resonance has been             switched off;         -   recording a plurality of data representative of said output             which has the form of a decay signal;         -   fitting said recorded data with a function S^(repl)(t,             δν_(i+d)) in which Z_(1Ueff)=constant replaces Z_(1U)             g(ν−ν_(o),w_(R)) in the fitting function S(t, δν_(i+d)).

More preferably, the method includes repeating the above steps, for a frequency ν_(i+2d)=ν_(i+d) of the electromagnetic radiation emitted by the laser, as long as the frequency of the electromagnetic radiation is included in [ν_(min), ν_(max)].

The frequency of the laser electromagnetic radiation is thus tuned so that different frequencies ν_(j) can be selected, all in a neighbourhood of the resonant frequency. In this way a plurality of SCAR decay signals, one for each frequency, are obtained. Each of these curve is fitted according to the method above described, always keeping the same value of Z_(1Ueff). Therefore, using the method of the invention, a plurality of values of

|σ_(γ) _(g) (δν_(i))|_(j=1) ^(m)

Can be obtained, for example m values for m different frequencies. In case for a given frequency more than a SCAR decay signal is collected, in order to improve the accuracy, the following steps are preferably repeated, with reference to the above description:

-   -   5) Repeat steps 1) to 3) up to j=1 . . . d*m, i.e. the total         number of detected SCAR decays for the i=1 . . . m step         detunings in the δν_(min)≤δν_(i)≤δν_(max) spectral range,         repeated by d times. This is made simultaneously to the SCAR         decays acquisitions.     -   6) As a result an array of fitted absorption decay rate         parameters is determined: |γ_(g)(δν_(i))|_(j=1) ^(d*m) with         uncertainties |σ_(γ) _(g) (δν_(i))|_(j=1) ^(d*m)

Advantageously, the method includes:

-   -   obtaining for each of a plurality of m frequencies ν_(j) with         δν_(j) from a δν_(min) to a δν_(max) a value of γ_(g)(δν_(j)),         and     -   fitting said m values of γ_(g)(δν_(j)) so as to obtain a value         of the first concentration of the target gas.

As clear from equation (3), the concentration is substantially the integrated absorption area under the γ_(g)(δν_(j)) profile.

More preferably said fitting of the m values of γ_(g)(δν_(j)) includes:

-   -   Selecting as a free parameters in the fitting a parameter which         takes into account the presence of other molecular absorptions         l=1 . . . n in addition to the target molecular resonance l=0.

More preferably said fitting of the m values of γ_(g)(δν_(j)) includes:

-   -   Selecting as a free parameters in the fitting a parameter which         takes into account the presence of a polynomial background         around the resonance frequency of the target transition ν₀.

Preferably, further consideration has to be performed before making the fitting of the m values of γ_(g)(δν_(j)). In the fitting of the SCAR decay signal with the S curve, the approximation of Z_(1Ueff)=const, has been made. This approximation however leads to a line profile modification, which will be dealt later.

Further approximations are as follows. First of all, around the molecular transition of interest, called here l=0, at the given frequency ν₀, there are other transitions l=1 . . . n of other gasses which are in the same cavity together with the target gas.

A SCAR apparatus is intended as an instrument designed and built to measure concentrations of a specific gas, g contained in a mixture of gases, g+og in conditions of total pressure PT and temperature T, and that interacts with saturating absorbing light, characterized by a power P₀ at the beginning of each SCAR decay event. The calibration procedures will be valid for concentration measurements of the targeted gas, one side, to parameterize accurately such absorption in terms of the spectral parameters, (i.e absorption intensity, resonance frequency and linewidth of targeted transition), and on the saturation parameters used to detect the SCAR spectrum. All of them, except absorption intensity which depends of the gas concentration, remains invariant if the thermodynamic conditions of the gas mixture and interacting light intensity do not change. Hence, they can be calibrated for such SCAR apparatus. On the other side, it means to remove, accurately, the contribution to the total integrated absorption area due to the absorptions of other lines of other molecules present into the mixture of the target gas and the other gasses, and resonant in the frequency interval of the detected SCAR spectrum. Hence, even for such transitions, spectral and saturation parameters can be calibrated and remain invariant for such SCAR apparatus.

Further, there is another contribution which is preferably to be taken into account, which is the background contribution around the resonance frequency.

The contribution of these other factors can be written as:

For each laser detuning:

S(t;p(δν_(i)))=B(δν_(i))+A _(d)(δν_(i))e ^(−γ) ^(c) ^((δν) ^(i) ^()t)ƒ(t;γ _(c)(δν_(i)),A _(d)(δν_(i)),γ_(g)(δν_(i)),p _(l)(ν_(i)))  (Eq. 35)

with:

-   -   p(δν_(i)) is the set of parameters that define the SCAR decay         behavior:

p(δν_(i))={B(δν_(i)),A _(d)(δν_(i)),γ_(c)(δν_(i)),γ_(g)(δν_(i)),p _(l)(ν_(i))},l=0 . . . n  (Eq. 36)

-   -   B(δν_(i)) is the detection background, in U     -   A_(d)(δν_(i)) is the amplitude of the decay signal at the begin         of the decay event, in U     -   γ_(c)(δν_(i)) and γ_(g)(δν_(i)) defined as above     -   p_(l)(ν_(i)) are n set of parameters that defines each absorbent         transition

p _(l)(ν_(i))={(Z _(1U))_(l) ,g _(l)(ν_(i)−ν_(l) ,w _(Rl))},l=0 . . . n  (Eq. 37)

-   -   ƒ the non-linear function that follow one of the rate equations         and depending of the gas conditions is described by:     -   homogeneous regime (w_(L)≥w_(G))

$\begin{matrix} {\frac{df}{dt} = {{- {\gamma_{g}\left( {\delta \; v_{i}} \right)}}{\sum\limits_{l = 0}^{n}\frac{\ln\left\lbrack {1 + {{A_{d}\left( {\delta \; v_{i}} \right)}\left( Z_{1U} \right)_{l}{{\overset{\_}{g}}_{l}\left( {{v_{i} - v_{l}},w_{Rl}} \right)}e^{{- {\gamma_{c}{({\delta \; v_{i}})}}}t}f}} \right\rbrack}{{A_{d}\left( {\delta \; v_{i}} \right)}\left( Z_{1U} \right)_{l}{{\overset{\_}{g}}_{l}\left( {{v_{i} - v_{l}},w_{Rl}} \right)}}}}} & \left( {{Eq}.\mspace{14mu} 38} \right) \end{matrix}$

-   -   which is a generalization of Eq. 23 for l=0 . . . n transitions         with (Z₀)_(l)=A_(d)(Z_(1U))_(l) g _(l).     -   inhomogeneous regime     -   still diffusive gas (w_(L)<w_(G))

$\begin{matrix} {\frac{df}{dt} = {{- {\gamma_{g}\left( {\delta \; v_{i}} \right)}}{\sum\limits_{l = 0}^{n}{\frac{2}{\sqrt{1 + {{A_{d}\left( {\delta \; v_{i}} \right)}\left( Z_{1U} \right)_{l}e^{{- {\gamma_{c}{({\delta \; v_{i}})}}}t}f}} + 1}f}}}} & \left( {{Eq}.\mspace{14mu} 39} \right) \end{matrix}$

-   -   which is a generalization of Eq. 25 for l=0 . . . n transitions         with (Z₀)_(l)=A_(d)(Z_(1U))_(l)     -   non-diffusive gas (w_(L)<<w_(G))

$\begin{matrix} {\frac{df}{dt} = {{- {\gamma_{g}\left( {\delta \; v_{i}} \right)}}{\sum\limits_{l = 0}^{n}{\frac{1.4256}{\sqrt{1 + {0.5{A_{d}\left( {\delta \; v_{i}} \right)}\left( Z_{1U} \right)_{l}e^{{- {\gamma_{c}{({\delta \; v_{i}})}}}t}f}} + 0.4256}f}}}} & \left( {{Eq}.\mspace{14mu} 40} \right) \end{matrix}$

-   -   which is a generalization of Eq. 30-32 for l=0 . . . n         transitions with (Z₀)_(l)=A_(d)(Z_(1U))_(l)

In the case of only the target transition is detected, Eqs. 38-40 are applied for l=0.

The target gas concentration is, deriving from eq. (3):

$\begin{matrix} {N_{g,{l = 0}} = {{\frac{1}{c^{2}{L_{S}(T)}}{\int{{\gamma_{g,{l = 0}}\left( {v - v_{0}} \right)}{dv}}}} = {\frac{1}{c^{2}{L_{S}(T)}}{\gamma_{g,{l = 0}}(0)}{\int{{g_{S\; 0}\left( {{v - v_{0}},w_{{RS}\; 0}} \right)}{dv}}}}}} & \left( {{Eq}.\mspace{14mu} 41} \right) \end{matrix}$

with:

-   -   L_(s)(T) is the line-strength of the target molecular transition         at temperature T, in cm/molecule     -   c is the vacuum velocity of light, in cm/s     -   γ_(g,l=0) is the contribution to the cavity decay rate due to         the target molecular transition

γ_(g,l=0)(ν−ν₀)=γ_(g,l=0)(0)g _(S0)(ν−ν₀ ,w _(RS,l=0))  (Eq. 42)

-   -   with     -   γ_(g,l=0)(0) is the contribution to the cavity decay rate due to         the target molecular transition at the resonance frequency ν₀,         in s⁻¹     -   g_(S0)(ν−ν₀,w_(RS0)) is the area normalized absorption profile         in saturation conditions, in Hz

g _(S0)(ν−ν₀ ,w _(RS0))=g _(S0)(0) g _(S0)(ν−ν_(O) ,w _(RS0))  (Eq. 43)

with

-   -   g_(S0)(0) is the area normalization factor, in Hz     -   g _(S0)(ν−ν₀,w_(RS0)) is the peak normalized absorption profile         in saturation conditions, with w_(RS0) saturation-modified         linewidth.

The function g _(S0)(ν−ν₀,w_(RS0)) differs from the peak normalized absorption profile in linear absorption conditions, g ₀(ν−ν₀,w_(R0)) by a factor R, due to the fact that the approximation

g _(S0)(ν−ν₀ ,w _(RS0))=Rg ₀(ν−ν₀ ,w _(R0))  (Eq. 44)

This R factor takes into account the line profile modifications due to Z_(1Ueff) approximation at the basis of this method.

Hence, from Eq. 41, the concentration of the target gas, N_(g,l=0) is:

$\begin{matrix} {N_{g,{l = 0}} = {\frac{{\gamma_{g,{l = 0}}(0)}{g_{S\; 0}(0)}}{c^{2}{L_{S}(T)}}R}} & \left( {{Eq}.\mspace{14mu} 45} \right) \end{matrix}$

The γ₀(0)g_(S0)(0) contribution in this equation is measured by fitting the set of final fitted parameters

γ_(g)(δ v_(i))_(δ v_(min))^(δ v_(max))

(i.e. the total absorption spectrum in the explored spectral range) to a function which takes into account the expected absorption profile:

γ(ν;p _(α))=Σ_(l=0) ^(n)γ_(l)(0)g _(Sl)(0) g _(Sl)(ν−ν_(l) ,w _(RSl))+Σ_(j=0) ^(p) B _(j)(ν−ν₀)^(j)  (Eq. 46)

with

-   -   p_(α) is the array of parameters that defines the expected         absorption profile including possible background

p _(α) ={p _(g) ,p _(bk)}  (Eq. 47)

-   -   with     -   p_(g) is the array of parameters that defines the expected         absorption profile due to the l=0 . . . n saturated absorptions     -   p_(bk) is the array of parameters that defines the possible         absorption background

These are thus the parameters that are preferably considered in the fitting procedure for talking into account the other saturated absorptions and the background.

More in detail, preferably:

p _(g)=(γ_(l)(0),g _(Sl)(0),ν_(l) ,w _(RSl))_(l=0) ^(n)  (Eq. 48)

with

-   -   γ_(l)(0) are the contributions to the cavity decay rate due to         absorptions from the l=0 . . . n molecular transitions at each         resonance frequency ν_(l), in s⁻¹     -   g_(Sl)(0) are the area normalization factors of the l=0 . . . n         molecular transitions, in Hz     -   ν_(l) are the resonance frequencies of the l=0 . . . n molecular         transitions, in Hz     -   w_(RSl) are the saturation-modified linewidths of the l=0 . . .         n molecular transitions, in Hz

p _(g) ={B _(j)}_(j=0) ^(p)

with

-   -   B₁ are the j=0 . . . p background coefficients of a polynomial         background around the resonance frequency of the target         transition ν₀.

Again, the fit of γ(ν; p_(α)) to

γ_(g)(δ v_(i))_(δ v_(min))^(δ v_(max))

uses the L-M method. At this respect for the p_(α) set of parameters only some of them are free to be fitted and the other are fixed to a value determined with the calibration procedure described below:

p _(α) ={p _(α) ^(fr) ,p _(α) ^(fx)}  (Eq. 49)

with

-   -   p_(α) ^(fr)=[{γ_(l)(0),g_(Sl)(0)}_(l=0) ^(n),{B_(j)}_(j=0) ^(p)]         is the array of free parameters     -   p_(α) ^(fx)={ν_(l),w_(RSl)}_(l=0) ^(n) is the array of fixed         parameters that are known are calibrated, for example following         the procedures described below.

Preferably, said fitting of the m values of γ_(g)(δν_(j)) is a is a least squares fitting. More preferably, said least square fitting uses the Levenberg-Marquardt (L-M) algorithm.

The recipe is preferably:

-   -   1) Initialization procedure:         -   Set the initial free parameters to:

p _(α) ^(ini)=[{γ_(l)(0),g _(Sl)(0)}_(l=0) ^(n) ,{B _(j)}_(j=0) ^(p)]^(ini)  (Eq. 50)

-   -   -   Set the p_(α) ^(fx) parameters to the calibrated values             following one of the calibration recipes (see below), and             keep them fixed during all fit procedure.         -   Other initializations: L-M limit convergence values,             λ_(max),δχ_(max) ²,χ_(min) ² and maximum number of             interactions N_(itr)

    -   2) Minimization loop on p_(α) ^(fr)

    -   a) p_(α) ^(try)=p_(α) ^(ini)

    -   b) p={p_(α) ^(try),p_(α) ^(fx)}. Calculate the χ²(p):

$\begin{matrix} {{\chi^{2}(p)} = {\sum\limits_{v_{i} = v_{\min}}^{v_{\max}}\frac{\left\lbrack {{\gamma_{g}\left( {v_{i} - v_{0}} \right)} - {\gamma \left( {v_{i};p} \right)}} \right\rbrack}{\sigma_{\gamma_{g}}^{2}\left( {v_{i} - v_{0}} \right)}}} & \left( {{Eq}.\mspace{14mu} 51} \right) \end{matrix}$

-   -   -   and gradient β(p) and the curvature matrix, |α|(p).

    -   c) Solve the linear Eqs. of the L-M algorithm for the present         problem, and calculate δp and χ²(p+δp)

    -   d) Follow the L-M convergence criteria on λ and/or χ² as         described in the L-M method. If not converging, come back to         step a) with an update value of λ and/or p_(α) ^(try)=p_(α)         ^(try)+δp

    -   e) At convergence, a fitted set of parameters equal to the last         guess is determined

p _(α) ^(fit)=[{γ_(l)(0),g _(Sl)(0)}_(l=0) ^(n) ,{B _(j)}_(j=0) ^(p)]^(fit) =p _(α) ^(try)  (Eq. 52)

-   -   3) Calculation of the target gas concentration:

$N_{0} = \frac{{\gamma_{0}^{fit}(0)}{g_{S\; 0}^{fit}(0)}}{c^{2}{L_{S}(T)}}$

more preferably corrected as detailed below.

As explained above, the introduction of a Z_(1Ueff) in the method allows deviations of the linear-absorption spectral profiles g _(Sl)(ν−ν_(l),w_(RSl)) with respect to the expected ones g _(l)(ν−ν_(l),w_(Rl)) taking into account the thermodynamic conditions of the gas mixture. The purpose of this calibration routines is determine the values of the parameters ν_(l),w_(RSl) of the g _(Sl)(ν−ν_(l),w_(RSl)) profile for the l detected transitions (strictly speaking g _(Sl) and g _(l) has the same functional behavior (we mean both are Voigt or Gaussians or Lorentzians, depending of the physical situation) but the linewidth is different, w_(RSl) for g _(Sl) and w_(Rl) for g _(l)). The routine is slightly different for the target transition (l=0) and for the others (l=1 . . . n).

In this way for the array of equation (49)

p _(α) ^(fx)={ν_(l) ,w _(RSl)}_(l=0) ^(n)

all parameters are known.

The calibration routine is depicted in FIG. 5.

-   -   For the target transition (l=0)

This routine uses a gas sample mixture enriched of the target gas, that is, with a second concentration of the target gas much higher than the first concentration. The strategy is to calculate the absorption spectrum |γ_(g) ^(↑S/N)(δν_(i))|_(j=1) ^(d*m) of this enriched sample in the δν_(min)≤δν_(i)≤δν_(max) spectral range, and fit it to the expected absorption behavior γ(ν; p_(α)) (Eq. 46), but now ν₀,w_(RS0) are free parameters. From eq. (49) p_(α)={p_(α) ^(fr), p_(α) ^(fx)} with

-   -   p_(α) ^(fr)=[γ₀ ^(↑S/N)(0),g_(S0)         ^(↑S/N)(0),{γ_(l)(0),g_(Sl)(0)}_(l=1)         ^(n),ν₀,w_(RS0),{B_(j)}_(j=0) ^(p)] is the array of free         parameters     -   p_(α) ^(fx)={ν₀,w_(RS0)}_(i=1) ^(n) is the array of fixed         parameters that must be known or must be calibrated following         the procedures described below.

The recipe is the following:

-   -   1) Calculate the spectrum {γ_(g) ^(↑S/N)(δν_(i)),σ_(γ) _(g)         ^(↑S/N)(δν_(i))}_(j=1) ^(d*m) by fitting the detected SCAR decay         signals of the enriched sample: {S_(exp) ^(↑S/N)(t;         δν_(i)),σ_(exp) ^(↑S/N)(t; δν_(i))} to a S(t; p(δν_(i)))         following the method above described.     -   2) Initialization procedure:         -   Set the initial free parameters to:

p _(α) ^(ini)=[γ₀ ^(↑S/N)(0),g _(S0) ^(↑S/N)(0),{γ_(l)(0),g _(Sl)(0)}_(l=1) ^(n),ν₀ ,w _(RS0) ,{B _(j)}_(j=0) ^(p)]^(ini)  (Eq. 53)

-   -   -   Set the p_(α) ^(fx) parameters to the calibrated values             following one of the calibration recipes (see below), and             keep them fixed during all fit procedure.         -   Other initializations: L-M limit convergence values,             λ_(max),δχ_(max) ²,χ_(min) ² and maximum number of             interactions N_(itr)

    -   3) Minimization loop on pr         -   a) p_(α) ^(try)=p_(α) ^(ini)         -   b) p={p_(α) ^(try), p_(α) ^(fx)}. Calculate the χ²(p)             (Eq. 21) and gradient β(p) and the curvature matrix, |α|(p)             for this problem.         -   c) Solve the linear Eqs. of the L-M algorithm for the             present problem, and calculate δp and χ²(p+δp)         -   d) Follow the L-M convergence criteria on λ and/or χ² as             described in the L-M algorithm. If not converging, come back             to step a) with an update value of λ and/or p_(α)             ^(try)=p_(α) ^(try)+δp         -   e) At convergence, a fitted set of parameters equal to the             last guess is determined

p _(α) ^(fit)[γ₀ ^(↑S/N)(0),g _(S0) ^(↑S/N)(0),{γ_(l)(0),g _(Sl)(0)}_(l=0) ^(n),ν₀ ,w _(RS0) ,{B _(j)}_(j=0) ^(p)]^(fit) =p _(α) ^(try)  (Eq. 54)

-   -   4) ν₀=ν₀ ^(fit) and w_(RS0)=w_(RS0) ^(fit) are the calibrated         values to be used in the method of the fitting of decay signal         at very low concentrations.     -   For the interference transitions (l=1 . . . n)

Differently from the others, this routine uses a gas sample mixture totally depleted of the target gas at the same P and T conditions. In this situation there is no absorption from the target gas for this experiment. The strategy is to calculate the absorption spectrum |γ_(g) ^(↓S/N)(δν_(i))|_(j=1) ^(d*m) of this totally depleted sample in the δν_(min)≤δν_(i)=δν_(max) spectral range, and fit it to the expected absorption behavior γ(ν; p_(α)) (Eq. 46), but now with all free parameters and without the parameters of the target transition (those with l=0)

p _(α) =p _(α) ^(fr)=[{γ_(l)(0),g _(Sl)(0),ν_(l) ,w _(RSl)}_(l=1) ^(n) ,{B _(j)}_(j=0) ^(p)]  (Eq. 55)

The recipe is the following:

-   -   1) Calculate the spectrum {γ_(g) ^(↓S/N)(δν_(i)),σ_(γ) _(g)         ^(↓S/N)(δν_(i))}_(j=1) ^(d*m) by fitting the detected SCAR decay         signals of the totally depleted sample: {S_(exp) ^(↓S/N)(t;         δν_(i)),σ_(exp) ^(↓S/N)(t; δν_(i))} to a S(t; p(δν_(i)))         following the method described for low concentrations.     -   2) Initialization procedure:         -   Set the initial free parameters to:

p _(α) ^(ini)=[{γ_(l)(0),g _(Sl)(0),ν_(l) ,w _(RSl)}_(l=0) ^(n) ,{B _(j)}_(j=0) ^(p)]^(ini)  (Eq. 56)

-   -   -   Other initializations: L-M limit convergence values,             λ_(max),δχ_(max) ²,χ_(min) ² and maximum number of             interactions N_(itr)

    -   3) Minimization loop on p_(α) ^(fr)         -   a) p_(α) ^(try)=p_(α) ^(ini)         -   b) p=p_(α) ^(try). Calculate the χ²(p) (Eq. 51) and gradient             β(p) and the curvature matrix, |α|(p) for this problem.         -   c) Solve the linear Eqs. of the L-M algorithm for the             present problem, and calculate δp and χ²(p+δp)         -   d) Follow the L-M convergence criteria on λ and/or χ² as             described in the L-M method. If not converging, come back to             step a) with an update value of λ and/or p_(α) ^(try)=p_(α)             ^(try)+δp         -   e) At convergence, a fitted set of parameters equal to the             last guess is determined

p _(α) ^(fit)=[{γ_(l)(0),g _(Sl)(0),ν_(l) ,w _(RSl)}_(l=1) ^(n) ,{B _(j)}_(j=0) ^(p)]^(fit) =p _(α) ^(try)  (Eq. 57)

-   -   4) ν_(l)=ν_(l) ^(fit) and w_(RSl)=w_(RSl) ^(fit) (l=1 . . . n)         are the calibrated values to be used in the ESPS method.

Advantageously, a result of said fitting is multiplied by a correcting factor R in order to obtain then concentration of the target gas, said correcting factor taking into account the line profile modifications due to the Z_(1Ueff)=constant approximation.

The equation above is preferably corrected in the following way:

$\begin{matrix} {N_{0} = {\frac{{\gamma_{0}^{fit}(0)}{g_{S\; 0}^{fit}(0)}}{c^{2}{L_{S}(T)}}R}} & \left( {{Eq}.\mspace{14mu} 58} \right) \end{matrix}$

with R to be determined with the calibration procedures described below.

Advantageously, the method includes the steps of:

-   -   Selecting the value of Z_(1Ueff)=constant to be introduced in         the fitting by means of the following step:         -   inserting said gas at a second concentration, wherein said             second concentration is at least 10 times said first             concentration in the resonant cavity;         -   tuning the frequency of said electromagnetic radiation             emitted by said laser source so as to fix it to a value             ν_(i) within a range of frequencies [ν_(min), ν_(max)]             including said molecular transition ν₀;         -   fixing the intensity of said electromagnetic radiation in             the cavity at a value much greater than the saturation             intensity I_(s) of the molecular transition to be detected;         -   irradiating said gas by means of said electromagnetic             radiation beam emitted by said laser source having said             fixed frequency ν_(i) and intensity in said resonant cavity;         -   coupling said electromagnetic radiation to said cavity so as             to obtain a laser-cavity resonance condition;         -   changing the frequency of the electromagnetic radiation             emitted by the laser so as to switch off the laser-cavity             resonance;         -   detecting an electromagnetic radiation beam in output from             said cavity after the laser-cavity resonance has been             switched off;         -   recording a plurality of data representative of said output             which has the form of a decay signal;         -   fitting the data of the recorded decay with a curve             S_(exp)(t, δν_(i)) which depends on the following             parameters:             -   B^(high conc)(δν_(i)) is the detection background, with                 (δν_(i)=ν_(i)−ν₀)             -   A_(d) ^(high conc)(δν_(i)) is the amplitude of the decay                 signal at the beginning of the decay event, g             -   γ_(c) ^(high conc)(δν_(i)) is the cavity decay rate due                 to non-resonant and non-saturable losses (empty cavity                 decay rate);             -   γ_(g) ^(high conc)(δν_(i)) is contribution of the                 targeted molecular transition to the decay signal;             -   Z_(1Ueff) ^(high conc)=Z_(1U) g(δν_(i),w_(R)), where g                 is the peak normalized line profile g(ν−ν₀,w_(R))                 centered at the molecular resonance frequency ν₀ and                 w_(R) is the HWHM width of the resonance, and                 w_(R)=w_(L) for a Lorentzian shape, w_(R)=w_(G) for a                 Gaussian shape, w_(R)={w_(L), w_(G)} for a Voigt shape;

$Z_{o} = {\frac{{CP}(0)}{{CP}_{s}} = {A_{d}Z_{1U}}}$

-   -   -   -   is the saturation parameter at the beginning of the                 decay event and at the frequency of the targeted                 molecular transition ν₀, P(0) is the intracavity power                 at the beginning of the decay signal; and

$P_{s} = {\frac{\pi \; w^{2}}{2}I_{s}}$

-   -   -   -   is the saturation power, where w is the spot size radius                 of the laser beam, i.e. the radius for which the                 amplitude of the field is 1/e times that of the axis and                 I_(s) the saturation intensity;

    -   Selecting as Z_(1Ueff) ^(high conc)=Z_(1Ueff)=value obtained         from the above fitting.

The invention in its first aspect is based on the approximation that Z_(1U) g(δν_(i),w_(R)) is a constant. The value to be put in the fitting procedure (and which is kept fixed for the whole procedure) is determined in a calibration routine. In the calibration routine, a SCAR measurement is performed, however an enriched gas is used. In other words, the same target gas at the same resonant frequency is investigated, however the concentration of the gas is much more than the concentration of the gas at which the standard gas measurements are performed. This second concentration of calibration is at least 10 times said first concentration.

The experimental decay curves used in most parts of this calibration routine are recorded with a sample gas mixture enriched with the target gas. This enrichment must be enough to measure a spectrum γ (δν) with high Signal to Noise ratio for the target transition, even with a CRD experiment in linear-absorption conditions (Z₀<<1). The thermodynamic conditions (Pressure, Temperature) of this enriched sample mixture are the same of the sample mixture used in the trace gas detection, that is, the same conditions in which the further measurement using the fitting of the invention is to be performed at very low concentrations.

This routine, depicted in FIG. 4, uses the data of a SCAR experiment performed in the enriched gas sample mixture to determine a value of Z_(1Ueff). The strategy is to detect the SCAR decay signals of this enriched sample: {S_(exp) ^(↑S/N)(t; δν_(i)),σ_(exp) ^(↑S/N)(t; δν_(i))} to be fitted with S(t; p(δν_(i))) following the method described according to the first aspect, but with the main difference that in this calibration routine directed to obtain a value of Z_(1Ueff), all parameters of p(δν_(i)) in the fitting of the so obtained decay signal, including Z_(1Ueff), are free.

p ^(fr)(δν_(i))={B ^(fr)(δν_(i)),A _(d) ^(fr)(δν_(i)),γ_(c) ^(fr)(δν_(i)),γ_(g) ^(fr)(δν_(i)),Z _(1Ueff) ^(fr)}  (Eq. 59)

Preferably, a similar routine as in the standard SCAR signal is used.

The procedure is the following:

-   -   1) Consider the detected SCAR curve {S_(exp) ^(↑S/N)(t;         δν_(i)),σ_(exp) ^(↑S/N)(t; δν_(i))}_(j)     -   2) Initialization procedure:         -   Set the initial free parameters to:

p _(j) ^(ini)(δν_(i))={B ^(ini)(δν_(i)),A _(d) ^(ini)(δν_(i)),γ_(c) ^(ini)(δν_(i)),γ_(g) ^(ini)(δν_(i)),Z _(1Ueff) ^(ini)}_(j)  (Eq. 60)

-   -   -   Other initializations: L-M limit convergence values,             λ_(max),δχ_(max) ²,χ_(min) ² and maximum number of             interactions N_(itr)

    -   3) Minimization loop on p_(j) ^(fr)(δν_(i))         -   Try as initial guess of the free parameters p_(j)             ^(try)(δν_(i))=p_(j) ^(ini)(δν_(i)) and uses the L-M routine             to obtain a fitted set of parameters p_(j) ^(fit)(δν_(i)).         -   a) p(δν_(i))=p_(j) ^(try)(δν_(i)). Calculate the             χ²(p(δν_(i))) (Eq. 7), gradient β(p(δν_(i))) and the             curvature matrix, |α|(p(δν_(i))) as described in item 3.a)             in sec. 1A).         -   b) Solve the linear Eqs. of the L-M algorithm for the             present problem, and calculate δp(δν_(i)) and             χ²(p(δν_(i))+δp(δν_(i)))         -   c) Follow the L-M convergence criteria on λ and/or χ² as             described in the L-M method. If not converging, come back to             step a) with an update value of λ and/or p_(j)             ^(try)(δν_(i))=p(δν_(i))+δp(δν_(i))         -   d) At convergence, a fitted set of parameters equal to the             last guess is determined

p _(j) ^(fit)(δν_(i))=p _(j) ^(try)(δν_(i))  (Eq. 61)

-   -   4) Repeat steps 1 to 3 up to j=1 . . . d*m, i.e. the total         number of detected SCAR decays for the i=1 . . . m step         detunings δν_(i) repeated by d times.     -   5) For the set of final fitted parameters |p_(j)         ^(fit)(δν_(i))|_(j=1) ^(d*m), calculate |χ²(p_(j)         ^(try)(δν_(i))|_(j=1) ^(d*m) and the sum of these χ² (p_(j)         ^(fit)(δν_(i)))

$\begin{matrix} {{S\left( \chi^{2} \right)} = {\sum\limits_{j = 1}^{d*m}{\chi^{2}\left( {p_{j}^{fit}\left( {\delta \; v_{i}} \right)} \right)}}} & \left( {{Eq}.\mspace{14mu} 62} \right) \end{matrix}$

-   -   -   and             -   If this calculation is performed by the first time (i.e.                 first convergence iteration), make S_(old)(χ²)=S(χ²) and                 come back to step 1 with an updated value Z_(1Ueff)                 ^(ini)=Z_(1Ueff) ^(fit)±δZ_(1Ueff)             -   In successive iterations, it searches for a minimum of                 S(χ²) in absolute and relative terms                 (δS(χ²)=S(χ2)−S_(old)(χ²)<δS_(lim)), i.e.:                 -   If S_(old)(χ²)<S(χ²)) and/or δS(χ²)>δS_(lim), make                     S_(old)(χ²)=S(χ²) and come back to step 1 with an                     updated value Z_(1Ueff)=Z_(1Ueff)±δZ_(1Ueff)                 -   If S_(old)(χ²)≥S(χ²)) and δS(χ²)≤δS_(lim), the                     procedure converges and the calibrated value of the                     effective saturation parameter to be used in the                     first aspect of the method of the invention is:                     Z_(1Ueff)=Z_(1Ueff) ^(fit)

Preferably, the value R is calculated performing a first and a second measurement of a concentration of the same target gas at the same temperature and pressure, the first measurement in saturation absorption and at a second concentration, wherein said second concentration is at least 10 times said first concentration, and said second measurement in linear absorption at said second concentration.

The modification of the absorption spectral profiles for the method according to the first aspect as explained above allows a wrong determination of the area under the target transition by a factor R, which is preferably calibrated. The strategy is to compare this area for a spectrum of the gas sample mixture enriched by the target gas in conditions of saturated-absorption |γ_(g) ^(↑S/N)(δν_(i))|_(j=1) ^(d*m) and in conditions of linear-absorption |γ_(li) ^(↑S/N)(δν_(i))|_(j=1) ^(d*m). The condition of linear absorption are that the intensity of the radiation inside the cavity is much lower than the saturation intensity, that is I<<_(s). These spectra are recorded by performing SCAR and CRD experiments (as mentioned, SCAR measurements are substantially CRD measurements in saturation conditions) in the same enriched gas sample, but at different laser intensities, by using the same apparatus, to obtain the SCAR decay curves {γ_(g) ^(↑S/N)(δν_(i)),σ_(γ) _(g) ^(↑S/N)(δν_(i))}_(j=1) ^(d*m) when I>>I_(S) and CRD decay curves {γ_(li) ^(↑S/N)(δν_(i)),σ_(γ) _(li) ^(↑S/N)(δν_(i))}_(j=1) ^(d*m), when I<<I_(s), respectively.

The recipe is the following:

-   -   1) Calculate the γ₀ ^(↑S/N)(0),g_(S0) ^(↑S/N)(0) parameters for         a spectrum {γ_(g) ^(↑S/N)(δν_(i)),σ_(γ) _(g)         ^(↑S/N)(δν_(i))}_(j=1) ^(d*m) following the procedure described         above to calculate the saturation-modified absorption profile of         the target transition (l=0). This procedure uses the SCAR decay         signals of the enriched sample: {S_(exp) ^(↑S/N)(t;         δν_(i)),σ_(exp) ^(↑S/N)(t; δν_(i))}, detected with the SCAR         apparatus in conditions of saturated absorption Z₀>>1.     -   2) Calculate the spectrum {γ_(li) ^(↑S/N)(δν_(i)),σ_(γ) _(li)         ^(↑S/N)(δν_(i))}_(j=1) ^(d*m) for the CRD detected with the SCAR         apparatus in conditions of linear-absorption Z₀<<1 for the         enriched sample mixture.     -   a. Consider the detected CRD curve {C_(exp) ^(↑S/N)(t;         δν_(i)),σ_(exp) ^(↑S/N)(t; δν_(i))}_(j) to be fitted to a         function C(t; p(δν_(i))) which is parameterized as:

C(t;p(δν_(i)))=B(δν_(i))+A _(d)(δν_(i))e ^(−(γ) ^(c) ^(+γ) ^(li) ^()(δν) ^(i) ^()t)  (Eq. 63)

-   -   -   with:         -   p(δν_(i)) is the set of parameters that define the CRD decay             behavior:

p(δν_(i))={B(δν_(i)),A _(d)(δν_(i)),(γ_(c)+γ_(li))(δν_(i))}  (Eq. 64)

-   -   -   by using the L-M algorithm with all parameters free to be             fitted p^(fr)(δν_(i))=p(δν_(i))

    -   b. Initialization procedure:         -   Set the initial free parameters to:

p _(j) ^(ini)((δν_(i))={B(δν_(i)),A _(d)(δν_(i)),(γ_(c)+γ_(li))(δν_(i))}_(j) ^(ini)  (Eq. 65)

-   -   -   Other initializations: L-M limit convergence values,             λ_(max),δχ_(max) ²,χ_(min) ² and maximum number of             interactions N_(itr)

    -   c. Minimization loop on p_(j) ^(fr)(δν_(i))         -   Try as initial guess of the free parameters p_(j)             ^(try)(δν_(i))=p_(j) ^(ini)(δν_(i)) and uses the L-M routine             to obtain a fitted set of parameters p_(j) ^(fit)(δν_(i)).         -   c.1 p(δν_(i))=p_(j) ^(try)(δν_(i)). Calculate the             χ²(p(δν_(i)))

$\begin{matrix} {{\chi^{2}\left( {p\left( {\delta \; v_{i}} \right)} \right)} = {\sum\limits_{t_{k} = 0}^{t_{\max}}\frac{\left( {{C_{\exp}\left( {t_{k};{\delta \; v_{i}}} \right)} - {C\left( {t_{k};{p\left( {\delta \; v_{i}} \right)}} \right)}} \right)^{2}}{\sigma_{\exp}^{2}\left( {t_{k};{\delta \; v_{i}}} \right)}}} & \left( {{Eq}.\mspace{14mu} 66} \right) \end{matrix}$

-   -   -   -   gradient β(p(δν_(i))) and the curvature matrix,                 |α|(p(δν_(i))) for this problem.

        -   c.2 Solve the linear Eqs. of the L-M algorithm for the             present problem, and calculate δp(δν_(i)) and             χ²(p(δν_(i))+δp(δν_(i)))

        -   c.3 Follow the L-M convergence criteria on λ and/or χ² as             described in the L-M method. If not converging, come back to             step a) with an update value of 2 and/or p_(j)             ^(try)(δν_(i))=p(δν_(i))+δp(δν_(i)).

        -   c.4 At convergence, a fitted set of parameters equal to the             last guess is determined

p _(j) ^(fit)(δν_(i))=p _(j) ^(try)(δν_(i))  (Eq. 67)

-   -   d. Repeat steps 2) to 4) up to j=1 . . . d*m, i.e. the total         number of detected CRD decays for the i=1 . . . m step detunings         in the δν_(min)≤δν_(i)≤δν_(max) spectral range, repeated by d         times. This is made simultaneously to the CRD decays         acquisitions.     -   e. As a result an array of fitted cavity decay rate parameters         is determined: |(γ_(c)+γ_(li))(δν_(i))|_(j=1) ^(d*m) with         uncertainties |σ_(γ) _(li) (δν_(i))|_(j=1) ^(d*m)     -   3) Calculate the (γ_(li))₀ ^(↑S/N)(0),g₀ ^(↑S/N)(0) parameters         by fitting this CRD spectrum |(γ_(c)+γ_(li))((δν_(i)),σ_(γ)         _(li) (δν_(i))|_(j=1) ^(d*m) to the expected linear-absorption

γ(ν;p _(li))=γ_(c)+(γ_(li))₀ ^(↑S/N)(0)g ₀ ^(↑S/N)(0) g ₀(ν−ν₀ ,w _(R0))+Σ_(l=1) ^(n)(γ_(li))_(l)(0)g _(l)(0) g _(l)(ν−ν_(l) ,w _(Rl))+Σ_(j=1) ^(p) B _(j)(ν−ν₀)^(j)  (Eq. 68)

-   -   with         -   p_(li) is the array of parameters that defines the expected             absorption profile including possible background

p _(li)=[γ_(c),(γ_(li)))₀ ^(↑S/N)(0),g ₀ ^(↑S/N)(0),ν₀ ,w _(R0),(γ_(li))_(l) ,g _(l)(0),ν_(l) ,w _(Rl))_(l=1) ^(n) ,{B _(j)}_(j=1) ^(p)]  (Eq. 69)

-   -   -   with         -   γ_(c) is the cavity decay rate due to losses different of             linear-absorption by the gas mixture, in s⁻¹         -   (γ_(li))₀ ^(↑S/N) is the contribution to the cavity decay             rate due to absorption of the target transition at resonance             frequency ν₀, in s⁻¹         -   g₀ ^(↑S/N)(0) is the area normalization factor of the             absorption of the target transition, in Hz         -   (γ_(li))(0) are the contributions to the cavity decay rate             due to absorptions of the l=1 . . . n interfering             transitions at each resonance frequency ν_(l), in s⁻¹         -   g_(l)(0) are the area normalization factors of the l=1 . . .             n interfering transitions, in Hz         -   ν_(l) are the resonance frequencies of the l=0 . . . n             transitions, in Hz         -   w_(Rl) are the linear-absorption linewidths of the l=0 . . .             n molecular transitions, in Hz         -   B_(j) are the j=1 . . . p background coefficients of a             polynomial background around the resonance frequency of the             target transition ν₀

Again, this fit uses the L-M method. At this respect for the p_(li) set of parameters only some of them are free to be fitted and the other are fixed to a value:

p _(li) =[p _(li) ^(fr) ,p _(li) ^(fx)]  (Eq. 70)

with

-   -   p_(li) ^(fr)=[γ_(c),(γ_(li))₀ ^(↑S/N)(0),g₀         ^(↑S/N)(0),ν₀,w_(R0),{(γ_(li))_(l),g_(l)(0)}_(l=1)         ^(n),{B_(j)}_(j=1) ^(p)] is the array of free parameters     -   p_(li) ^(fx)=[{ν_(l),w_(RSl)}_(l=1) ^(n)] is the array of fixed         parameters

The recipe is:

-   -   a. Initialization procedure:         -   Set the initial free parameters to:

p _(li) ^(ini)=[γ_(c)(γ_(li))₀ ^(↑S/N)(0),g ₀ ^(↑S/N)(0),ν₀ ,w _(R0),{(γ_(li))_(l) ,g _(l)(0)}_(l=1) ^(n) ,{B _(j)}_(j=1) ^(p)]^(ini)  (Eq. 71)

-   -   -   Set the pix parameters to the values calculated from known             data reported by molecular databases         -   Other initializations: L-M limit convergence values,             λ_(max),δχ_(max) ²,χ_(min) ² and maximum number of             interactions N_(itr)

    -   b. Minimization loop on p_(α) ^(fr)         -   b.1. p_(li) ^(try)=p_(li) ^(ini)         -   b.2. p={p_(li) ^(try),p_(li) ^(fx)}. Calculate the χ²(p)

$\begin{matrix} {{\chi^{2}(p)} = {\sum\limits_{v_{i} = v_{\min}}^{v_{\max}}\frac{\left\lbrack {{\left( {\gamma_{c} + \gamma_{li}} \right)\left( {v_{i} - v_{0}} \right)} - {\gamma \left( {v_{i};p} \right)}} \right\rbrack}{\sigma_{\gamma_{l\; i}}^{2}\left( {v_{i} - v_{0}} \right)}}} & \left( {{Eq}.\mspace{14mu} 72} \right) \end{matrix}$

-   -   -   -   and gradient β(p) and the curvature matrix, |α|(p) for                 this problem.

        -   b.3. Solve the linear Eqs. of the L-M algorithm for the             present problem, and calculate δp and χ²(p+δp)

        -   b.4. Follow the L-M convergence criteria on λ and/or χ² as             described in the L-M method. If not converging, come back to             step a) with an update value of λ and/or p_(α) ^(try)=p_(α)             ^(try)+δp

        -   b.5. At convergence, a fitted set of parameters equal to the             last guess is determined

p _(li) ^(fit)=[γ_(c),(γ_(li))₀ ^(↑S/N)(0),g ₀ ^(↑S/N)(0),ν₀ ,w _(R0),{(γ_(li))_(l) ,g _(l)(0)}_(l=1) ^(n) ,{B _(j)}_(j=0) ^(p)]^(fit) =p _(α) ^(try)  (Eq. 53)

-   -   4) Calculate the R-factor to be used in the calculation of         concentrations by using

$N_{0} = {\frac{{\gamma_{0}^{fit}(0)}{g_{S\; 0}^{fit}(0)}}{c^{2}{L_{S}(T)}}R\text{:}}$

$\begin{matrix} {R = \frac{{\gamma_{0}^{{\uparrow S}/N}(0)}{g_{S\; 0}^{{\uparrow S}/N}(0)}}{\left( \gamma_{li} \right)_{0}^{{\uparrow S}/N}(0){g_{0}^{{\uparrow S}/N}(0)}}} & \left( {{Eq}.\mspace{14mu} 74} \right) \end{matrix}$

with

-   -   γ₀ ^(↑S/N)(0),g_(S0) ^(↑S/N)(0) from p_(α) ^(fit) (i.e. Eq. 54)     -   (γ_(li))₀ ^(↑S/N)(0), g₀ ^(↑S/N)(0) from p_(li) ^(fit) (i.e. Eq.         72)

According to a second aspect, the invention relates to a method of ring-down spectroscopy in saturated-absorption condition, for measuring a first concentration of a target gas through a measurement of the spectrum of a molecular transition of said target gas, the target gas being in a mixture together with other gasses the method comprising the steps of:

-   -   Repeating for a number m of fixed frequencies vi spaced each         other by a of a frequency step in a range of frequencies         including the frequency of the molecular transition of the         target gas, the gases in the mixture having l=0, . . . , n         absorptions in the measured spectral range, where l=0 is the         target absorption, the following steps:         -   Repeating d times at the same frequency the following steps:             -   inserting said gas whose first concentration is to be                 measured in a resonant cavity comprising two or more                 reflecting mirrors arranged so as to form a closed                 optical path for an electromagnetic radiation emitted by                 a laser source;             -   tuning the frequency of said electromagnetic radiation                 emitted by said laser source so as to fix it to a value                 ν_(i) within a range of frequencies [ν_(min), ν_(max)]                 including said molecular transition ν₀;             -   fixing the intensity of said electromagnetic radiation                 in the cavity at a value much greater than the                 saturation intensity I_(s) of the molecular transition                 to be detected;             -   irradiating said gas by means of said electromagnetic                 radiation beam emitted by said laser source having said                 fixed frequency ν_(i) and intensity in said resonant                 cavity;             -   coupling said electromagnetic radiation to said cavity                 so as to obtain a laser-cavity resonance condition;             -   changing the frequency of the electromagnetic radiation                 emitted by the laser so as to switch off the                 laser-cavity resonance;             -   detecting an electromagnetic radiation beam in output                 from said cavity after the laser-cavity resonance has                 been switched off;             -   recording a plurality of data representative of said                 output, obtaining a decay signal for the fixed                 frequency;     -   collecting the d*m SCAR decay signals obtained;     -   Fitting at the same time the d*m SCAR decay signals with d*m         fitting curves, considering a fitting curve S(t; p(ν))=[S(t;         p(ν_(i))_(j))]_(j=1) ^(d*m) for the recorded decay signals which         depends on the following parameters p(ν)=[p(ν_(i))_(j)]_(j=1)         ^(d*m):

p(ν)=[p(ν)_(j)]_(j=1) ^(d*m) =[B _(j) ,A _(d) _(j) ,γ_(c) _(j) ,[{γ_(l)(δν_(il)),Z _(l1U) g _(l)(δν_(il) ,w _(Rl))}_(j)]_(l=0) ^(n)]_(j=1) ^(d*m) ,i=1 . . . m

-   -   [B_(j)(δν_(il))]_(j=1) ^(d*m) are the detection backgrounds of         the j=1 . . . d*m recorded SCAR decays; and δν_(il)=ν_(i)−ν_(l)         is the detuning of the i (=1 . . . m) scanned frequency of the         laser radiation with respect to the resonance frequency of the l         resonant transition,     -   [A_(d) _(j) (δν_(il))]_(j=1) ^(d*m) are the amplitudes of the         j=1 . . . d*m recorded SCAR decays,     -   [γ_(c) ^(j)(δν_(il))]_(j=1) ^(d*m) are the cavity decay rates of         the j=1 . . . d*m recorded SCAR decays due to non-resonant and         non-saturated gas absorption losses,     -   [[γ_(l)(δν_(il))_(j)]_(l=0) ^(n)]_(j=1) ^(d*m) are the arrays of         the absorption decay rates of the l=0 . . . n saturated         transitions at the δν_(il) detuning for the j=1 . . . d*m         recorded SCAR decays,     -   [[{Z_(1Ul) g _(l)(ν_(i)−ν_(l),w_(Rl))}_(j)]_(l=0) ^(n)]_(j=1)         ^(d*m) are the arrays of the l=0 . . . n saturation parameters         of the j=1 . . . d*m recorded SCAR decays; g         _(l)(ν_(i)−ν_(l),w_(Rl)) is the peak normalized line profile         centered at the molecular resonance frequency ν_(l) and w_(Rl)         is the HWHM width of the l-resonance;

$Z_{o} = {\frac{{CP}(0)}{{CP}_{s}} = {A_{d}Z_{1U}}}$

-   -   is the saturation parameter at the beginning of the decay event         and at the frequency of the targeted molecular transition ν₀,         P(0) is the intracavity power at the beginning of the decay         signal; and

$P_{s} = {\frac{\pi \; w^{2}}{2}I_{s}}$

-   -   is the saturation power, where w is the spot size radius of the         laser beam, i.e. the radius for which the amplitude of the field         is 1/e times that of the axis and I_(s) the saturation         intensity;     -   Dividing the above mentioned parameters in a first and a second         group;     -   fitting said recorded data with a function S(t; p(ν))=[S(t;         p(ν_(i))_(j))]_(j=1) ^(d*m) by         -   keeping the first parameter group fixed and equal to a set             of pre-determined values and performing a first fitting only             considering the second group as free parameters in a first             fitting step; and         -   keeping the second parameter group fixed to a value given in             the first fitting step and performing a second fitting             considering only the second group as free parameters in a             second fitting step.

Differently from the first aspect, the second aspect of the invention follows a global approach to determine the absorption area of the target gas from all detected SCAR decays

[{S_(exp)(t; v_(i)), σ_(exp)(t; v_(i))}_(j)]_(j = 1)^(d * m).

It means that these input data must be known a priori, and hence when the complete frequency scan of the experiment is performed. In this sense, the second aspect of the invention gives a concentration measurement of the target gas a posteriori (i.e. not in real time with the data acquisition). On the contrary, it gives under certain circumstances more precise results by taking into account better the spectral behavior, and hence avoiding some correlations.

The acquisition of the various SCAR decay signals takes place as described already with reference to the first aspect of the invention. The details therefore are not herein repeated. Everything regarding the acquisition is substantially the same as before, the only difference is that the acquisition ends when a plurality of SCAR decays for m different frequencies is taken. Preferably for a single frequency more than a single SCAR decay signal is recorded, however this is not necessary in the present invention.

The second aspect of the invention includes a step of fitting, all together, this large number of input data

{S_(exp)(t; v), σ_(exp)(t; v)} = [{S_(exp)(t; v_(i)), σ_(exp)(t; v_(i))}_(j)]_(j = 1)^(d * m)

(i.e. I=1 . . . m frequency steps repeated d times) with an equally large number of fit functions S(t; p(ν))=[S(t; p(ν_(i))_(j))]_(j=1) ^(d*m) with each element p(ν_(i))_(j) of the total parameter array p(ν) having seven parameters to be fitted in the most simple case, where interference lines are not present in the spectral range ν₁≤ν_(i)≤ν_(m). It means that the fit routine should manage, at least, d*m*7 parameters This situation is even worse if interfering absorptions are included, that is further absorption lines with l=1 . . . n, as already discussed with reference to the first aspect of the invention, incrementing consequently the number of parameters to be fitted.

To simplify the problem and to better take into account the effect of the line profile g in the saturation parameter Z₀ and its correlation with γ_(g), the second aspect of the method of the invention considers a set of global parameters p_(gb)(ν) equal for all [S(t; p(ν_(i))_(j))]_(j=1) ^(d*m) functions and j=1 . . . d*m sets of local parameters [p_(lc) _(j) ]_(j=1) ^(d*m) to take into account small frequency-dependent variations around a mean value of some global parameters:

p(ν)=[p _(gb)(ν),[p _(lc) _(j) ]_(j=1) ^(d*m)]  (Eq.75)

The experimental situation is a SCAR experiment, that is, in saturation conditions, in a gas mixture that can have l=0 . . . n absorptions in the measured spectral range. l=0 is the resonant transition of the target gas and l=1 . . . n are the resonant interfering transitions of other gases present in the sample.

The total decay function S(t; p(ν)) is parameterized in the following way:

                                        (Eq.  76) ${{S\left( {t;{p(v)}} \right)} = {\left\lbrack {S\left( {t;{p\left( v_{i} \right)}_{j}} \right)} \right\rbrack_{j = 1}^{d*m} = \left\lbrack {B_{j} + {A_{d_{j}}e^{{- \gamma_{c_{j}}}t}{f_{j}\left( {{t;A_{d_{j}}},\gamma_{c_{j}},\left\lbrack \left\{ {{\gamma_{l}\left( {\delta \; v_{il}} \right)},{Z_{l\; 1U}{{\overset{\_}{g}}_{l}\left( {{\delta \; v_{il}},w_{Rl}} \right)}}} \right\}_{j} \right\rbrack_{l = 0}^{n}} \right)}}} \right\rbrack_{j = 1}^{d*m}}},\mspace{79mu} {i = {1\mspace{14mu} \ldots \mspace{14mu} m}}$

with

-   -   δν_(il)=ν_(i)−ν_(l) is the detuning of the i (=1 . . . n)         scanned frequency of the laser radiation with respect to the         resonance frequency of the l resonant transition, in Hz     -   p(ν) is the set of parameters that define the SCAR decay         behavior in the explored spectral range:

p(ν)=[p(ν_(i))_(j)]_(j=1) ^(d*m) =[B _(j) ,A _(d) _(j) ,γ_(c) _(j) ,[{γ_(l)(δν_(il)),Z _(1U) g _(l)(δν_(il) ,w _(Rl))}_(j)]_(l=0) ^(n)]_(j=1) ^(d*m) ,i=1 . . . m   (Eq.77)

-   -   -   [B_(j)]_(j=1) ^(d*m) are the detection backgrounds of the             j=1 . . . d*m recorded SCAR decays; in U         -   [A_(d) _(j) ]_(j=1) ^(d*m) are the amplitudes of the j=1 . .             . d*m recorded SCAR decays, in U         -   [γ_(c) _(j) ]_(j=1) ^(d*m) are the cavity decay rates of the             j=1 . . . d*m recorded SCAR decays due to non-resonant and             non-saturated gas absorption losses, in s⁻¹         -   [[γ_(l)(δν_(il))_(j)]_(l=0) ^(n)]_(j=1) ^(d*m) are the             arrays of the absorption decay rates of the l=0 . . . n             saturated transitions at the δν_(il) detuning for the j=1 .             . . d*m recorded SCAR decays, in s⁻¹         -   [[{Z_(1Ul) g _(l)(ν_(i)−ν_(l),w_(Rl))}_(j)]_(l=0)             ^(n)]_(j=1) ^(d*m) are the arrays of the l=0 . . . n             saturation parameters at 1U of the j=1 . . . d*m recorded             SCAR decays, in U⁻¹

    -   ƒ=[ƒ_(j)]_(j=1) ^(d*m) are the non-linear functions that follows         one of the below rate equations depending of the gas conditions:         -   homogeneous regime (w_(L)≥w_(G))

$\begin{matrix} \left\lbrack {\frac{{df}_{j}}{dt} = {- {\sum\limits_{l = 0}^{n}{{\gamma_{l}\left( {\delta \; v_{il}} \right)}_{j}\frac{\ln \left\lbrack {1 + {A_{d_{j}}Z_{l\; 1U}{{\overset{\_}{g}}_{l}\left( {{\delta \; v_{il}},w_{Rl}} \right)}e^{{- \gamma_{c_{j}}}t_{f_{j}}}}} \right\rbrack}{A_{d_{j}}Z_{l\; 1U}{{\overset{\_}{g}}_{l}\left( {{\delta \; v_{il}},w_{Rl}} \right)}e^{{- \gamma_{c_{j}}}t}}}}}} \right\rbrack_{j = 1}^{d*m} & \left( {{Eq}.\mspace{14mu} 78} \right) \end{matrix}$

-   -   -   inhomogeneous regime             -   still diffusive gas (w_(L)<w_(G))

$\begin{matrix} \left\lbrack {\frac{{df}_{j}}{dt} = {- {\sum\limits_{l = 0}^{n}{{\gamma_{l}\left( {\delta \; v_{il}} \right)}_{j}\frac{2}{\sqrt{1 + {A_{d_{j}}Z_{l\; 1U}e^{{- \gamma_{c_{j}}}t}f_{j}}} + 1}f_{j}}}}} \right\rbrack_{j = 1}^{d*m} & \left( {{Eq}.\mspace{14mu} 79} \right) \end{matrix}$

-   -   -   -   non-diffusive gas (w_(L)<<w_(G))

$\begin{matrix} \left\lbrack {\frac{{df}_{j}}{dt} = {- {\sum\limits_{l = 0}^{n}{{\gamma_{l}\left( {\delta \; v_{il}} \right)}_{j}\frac{1.4256}{\sqrt{1 + {0.5A_{d_{j}}Z_{l\; 1U}e^{{- \gamma_{c_{j}}}t}f_{j}} + 0.4256}}f_{j}}}}} \right\rbrack_{j = 1}^{d*m} & \left( {{Eq}.\mspace{14mu} 80} \right) \end{matrix}$

It is to be noted that more than one of those parameters is strongly correlated for all SCAR decay experimental curves. Further other parameters such as [B_(j)]_(j=1) ^(d*m),

[A_(d_(j))]_(j = 1)^(d * m), [γ_(c_(j))]_(j = 1)^(d * m)

depends very slight on the frequency ν_(i) and they are quasi-equal for all j curves. Therefore, for the B_(j), A_(d) _(j) and γ_(c) _(j) parameters, which have small variations from decay to decay, it can be written

B _(j) =B+δB _(j)  (Eq. 81)

-   -   B is the mean value of the background of the j=0 . . . d*m         recorded SCAR decays, in U     -   {δB_(j)}_(j=1) ^(d*m) are small variations around B for each j=0         . . . d*m recorded SCAR decay, in U

A _(d) _(j) =Ā _(d) +δA _(d) _(j)   (Eq. 82)

-   -   Ā_(d) is the mean value of the amplitude of the j=0 . . . d*m         recorded SCAR decays, in U

{δ A_(d_(j))}_(j = 1)^(d * m)

-   -   are small variations around Ā_(d) for each j=0 . . . d*m         recorded SCAR decay, in U

γ_(c) _(j) =γ _(c)+δγ_(c) _(j)   (Eq. 83)

-   -   γ _(c) is the mean value of the cavity decay rate of the j=0 . .         . d*m recorded SCAR decays due to non-absorbent gas and         non-saturated absorption losses, in s⁻¹

{δ γ_(c_(j))}_(j = 1)^(d * m)

-   -   are small variations around γ _(c) for each j=0 . . . d*m         recorded SCAR decay, in s⁻¹

The parameters are thus divided in a first group of parameters, the first group preferably comprising:

$\begin{matrix} {\left\lbrack p_{{lc}_{j}} \right\rbrack_{j = 1}^{d*m} = \left\lbrack {{\delta \; B_{j}},{\delta \; A_{d_{j}}},{\delta \; \gamma_{c_{j}}}} \right\rbrack_{j = 1}^{d*m}} & \left( {{Eq}.\mspace{14mu} 84} \right) \end{matrix}$

And the other in a second group. Due to the difference in variation of these parameters, that is, due to the fact that the parameters of the first group have a “local” variation and the parameter of the second group have a “global” variation, it is assumed that the procedure of the invention can be performed in two separate fits.

In the first fit the second group of parameters is kept fixed. Only the local parameters, or first group, are real “free” parameters in the fitting procedure. At the end of the fitting procedure a series of values for the local parameters is obtained, in a much simplified way because not all parameters are considered in a single fitting procedure. Further, using these values, i.e. the values in output of the fitting procedure, in the second fitting step, the global parameters now becomes the real “free parameters” in the second fitting, where the values of the local parameters are kept fixed to the value at the outcome of the first fitting procedure.

Preferably, the method includes the step of calculating the first concentration of the target gas N₀ by:

$\begin{matrix} {N_{0} = \frac{{\gamma_{0}^{fit}(0)}{g_{0}^{fit}(0)}}{c^{2}{L_{S}(T)}}} & \; \end{matrix}$

-   -   Where c is the speed of light in vacuum and     -   L_(s)(T) is the line-strength of the absorbent molecular         transition at the temperature T.

The concentration of the target gas is obtained by eq. (3) and the values obtained from the fitting procedure (a global and a local value—the fitting values are indicated with an apex “fit”) without any further correction procedure.

Preferably, the method of the invention includes

-   -   Introducing two further parameters to the second group of         parameters, said parameters being:         -   d is the contribution to the cavity decay rate due to             residual detection non-linearities, in s⁻¹         -   Z_(d) is the “equivalent”-saturation parameter at 1U due to             residual detection non-linearities, in U⁻¹     -   So that the second group of parameters is equal to:

p _(gb)(ν)=[ B,Ā _(d),γ _(c) ,Z _(d) ,d,[Z _(1Ul),γ_(l)(0),g _(l)(0), g _(l)(δν_(il) ,w _(Rl))]_(l=0) ^(n)]

Where γ_(l)(δν_(il))_(j)=γ_(l)(0)g _(l)(0) g _(l)(δν_(il) ,w _(Rl))  (eq. 85)

with

-   -   {γ_(l)(0)}_(l=0) ^(n)={cα_(l)(0)}_(l=0) ^(n) are the absorption         decay rates for each transition l=0 . . . n at the resonance         frequency ν_(i)=ν_(l), in s⁻¹,     -   {g_(l)(0)}_(l=0) ^(n) are the area normalization factors for         each transition l=0 . . . n at the resonance frequency         ν_(i)=ν_(l),     -   {g _(l)(ν_(i)−ν_(l),w_(Rl))}_(l=0) ^(n) are the peak normalized         line profiles for each transition l=0 . . . n, centered at the         resonance frequency ν_(i) and linewidth w_(Rl).

d and Z_(d) are, respectively, the area and saturation parameter of a possible spectral background. For this distortion background, it has been assumed a non-linear dependence equal to that of the saturated-absorption that characterized the gas absorption but with a constant saturation parameter for all probed frequencies. It is assumed that it is a background due to not-cancelled residual-random non-linearities. Adding these parameters in the fit allows to avoid any further correction to the fitting results.

Further, preferably the method includes adding a further additional parameter [C_(l) ^(γ)]_(l=0) ^(n), which is the amplitude factors of a function G_(j)(ν_(i); [C_(l) ^(γ)]_(l=0) ^(n)), proportional to the spectrum profile [g _(l)(δν_(il),w_(Rl))]_(l=0) ^(n)

G _(j)(ν_(i) ;[C _(l) ^(γ)]_(l=0) ^(n))=Σ_(l=0) ^(n) C _(l) ^(γ) g _(l)(δν_(il) ,w _(Rl))  (Eq.85a)

-   -   is used to update the local parameters

[δ γ_(c_(j))]_(j = 1)^(d * m).

The [C_(l) ^(γ)]_(l=0) ^(n) are amplitude factors to take into account the correlation between γ_(c) and γ_(g) for each spectral frequency ν_(l), both for the targeted absorption and for the l absorptions of other gases, respectively. They have not a physical meaning and they were introduced for improve calculation strategy. By adding this parameter, the cross correlation between γ_(c) and γ_(g), is taken into account in the local fits. Indeed, the global/local approach in the parametrization of the fit function imply that γ_(g)(ν_(i)) behavior is taken into account by the global parameters, whereas γ_(c)(ν_(i)) behavior is taken into account by the local parameters. Because the global parameters are kept fixed during fit of local parameters, and viceversa, this cross correlation is slowly propagated in the proper way in each of the correlated parameters, doing the GLFS routine very slow to converge. The introduction of the auxiliar global parameters allows us to transfer in a proper and quicker fashion this correlation, accelerating the convergence process.

Preferably, the method of the invention preferably includes parametrizing the non-linear function ƒ as:

-   -   homogeneous regime (w_(L)≥w_(G))

$\begin{matrix} \begin{bmatrix} {\frac{{df}_{j}}{dt} = {{- {\sum\limits_{l = 0}^{n}{{\gamma_{l}(0)}{g_{l}(0)}\frac{\ln \begin{bmatrix} {1 + \left( {{\overset{\_}{A}}_{d} + {\delta \; A_{d_{j}}}} \right)} \\ {Z_{l\; 1U}{{\overset{\_}{g}}_{l}\left( {{\delta \; v_{il}},w_{Rl}} \right)}e^{- {({{\overset{\_}{\gamma}}_{c} + {\delta \; \gamma_{c_{j}}}})}^{t}}f_{j}} \end{bmatrix}}{\left( {{\overset{\_}{A}}_{d} + {\delta \; A_{d_{j}}}} \right)Z_{l\; 1U}e^{- {({{\overset{\_}{\gamma}}_{c} + {\delta \; \gamma_{c_{j}}}})}^{t}}}}}} -}} \\ {d\frac{\ln \left\lbrack {1 + {\left( {{\overset{\_}{A}}_{d} + {\delta \; A_{d_{j}}}} \right)Z_{d}e^{- {({{\overset{\_}{\gamma}}_{c} + {\delta \; \gamma_{c_{j}}}})}^{t}}f_{j}}} \right\rbrack}{\left( {{\overset{\_}{A}}_{d} + {\delta \; A_{d_{j}}}} \right)Z_{d}e^{- {({{\overset{\_}{\gamma}}_{c} + {\delta \; \gamma_{c_{j}}}})}^{t}}}} \end{bmatrix}_{j = 1}^{d*m} & \left( {{eq}.\mspace{14mu} 86} \right) \end{matrix}$

and g _(l)(ν_(il),w_(Rl)) is a Voigt function

-   -   inhomogeneous regime     -   still diffusive gas (w_(L)<w_(G))

$\begin{matrix} \left\lbrack {\frac{{df}_{j}}{dt} = {\begin{bmatrix} {{- {\sum\limits_{l = 0}^{n}\frac{2\; {\gamma_{l}(0)}{g_{l}(0)}{{\overset{\_}{g}}_{l}\left( {{\delta \; v_{il}},w_{Rl}} \right)}}{\sqrt{1 + {\left( {{\overset{\_}{A}}_{d} + {\delta \; A_{d_{j}}}} \right)Z_{l\; 1U}e^{- {({{\overset{\_}{\gamma}}_{c} + {\delta \; \gamma_{c_{j}}}})}^{t}}f_{j}} + 1}}}} -} \\ \frac{2d}{\sqrt{1 + {\left( {{\overset{\_}{A}}_{d} + {\delta \; A_{d_{j}}}} \right)Z_{d}e^{- {({{\overset{\_}{\gamma}}_{c} + {\delta \; \gamma_{c_{j}}}})}^{t}}f_{j}} + 1}} \end{bmatrix}f_{j}}} \right\rbrack_{j = 1}^{d*m} & \left( {{Eq}.\mspace{14mu} 87} \right) \end{matrix}$

and (δν_(il),w_(Rl)) is a Gaussian function

-   -   non-diffusive gas (w_(L)<<w_(G))

$\begin{matrix} {\begin{bmatrix} {\frac{{df}_{j}}{dt} =} \\ {\begin{bmatrix} {{- {\sum\limits_{l = 0}^{n}\frac{1.4256\; {\gamma_{l}(0)}{g_{l}(0)}{{\overset{\_}{g}}_{l}\left( {{\delta \; v_{il}},w_{R\; l}} \right)}}{\sqrt{1 + {0.5\left( {{\overset{\_}{A}}_{d} + {\delta \; A_{d_{j}}}} \right)Z_{l\; 1U}e^{{- {({{\overset{\_}{\gamma}}_{c} + {\delta \; \gamma_{c_{j}}}})}^{t}}f_{j}}}} +}}} -} \\ 0.4256 \\ \frac{1.4256d}{\sqrt{1 + {0.5\left( {{\overset{\_}{A}}_{d} + {\delta \; A_{d_{j}}}} \right)Z_{d}e^{{- {({\gamma_{c} + {\delta \; \gamma_{c_{j}}}})}^{t}}f_{j}}}} + 0.4256} \end{bmatrix}f_{j}} \end{bmatrix}_{j = 1}}^{d*m} & \left( {{Eq}.\mspace{14mu} 88} \right) \end{matrix}$

-   -   and g _(l)(δν_(il),w_(Rl)) is a Gaussian function.

Preferably, said first fitting step is a least squares fitting. More preferably, said least square fitting uses the Levenberg-Marquardt (L-M) algorithm.

The method according to the second aspect of the invention follows a sequential strategy to calculate the final parameter values by fitting alternatively the local parameters while keeping fixed the global ones, and then by fitting the global parameters while keeping fixed the local ones. This sequence is repeated up to final convergence of all procedure as explained below.

Each element S(t; p(ν_(i))_(j)) of the total decay function S(t; p(ν)) is calculated from the following parametrization by using Eqs. 81-85:

S(t;p(ν))=[ B+δB _(j)+(Ā _(d) ±δA _(d) _(j) )e ^(−γ) ^(c) ^((δν) ^(i) ^()t)·ƒ_(j)(t;Ā _(d),γ _(c) ,Z _(d) ,d,δA _(d) _(j) ,δγ_(c) _(j) ,[C _(l) ^(γ) ,Z _(1Ul),γ_(l)(o),g ₁(0), g _(l)(δν_(il) ,w _(Rl))]_(l=0) ^(n))]  (Eq. 89)

with the set of parameters p(ν) reduced to a set of global parameters, p_(gb)(ν) and a set of arrays of local parameters, [p_(lc)]_(j=1) ^(d*m) to be determined by fitting the recorded decays:

p _(gb)(ν)=[ B,Ā _(d),γ _(c) ,Z _(d) ,d,[C _(l) ^(γ) ,Z _(1Ul),γ_(l)(0),g _(l)(0), g _(l)(δν_(il) ,w _(Rl))]₁₌₀ ^(n)]  (Eq.90)

[p _(lc) _(j) ]_(j=1) ^(d*m) =[δB _(j) ,δA _(d) _(j) ,δγ_(c) _(j) ]_(j=1) ^(d*m)  (Eq.91)

after the ƒ function is evaluated by the numerical integration of the rate equation that applies to the case under analysis (Eqs. 86-88). A fourth-order Runge-Kutta (RK4) routine for such numerical integration has be used, but other integration routines can be used as well.

The first fitting step where a first group of parameters is free and the second group is fixed is depicted in the flowchart diagram of FIG. 6 and it is described in the following:

-   -   1) Consider all the detected SCAR curves [{S_(exp)(t;         ν_(i)),σ_(exp)(t; ν_(i))}_(j)]_(j=1) ^(d*m) in the frequency         interval ν₁<ν_(i)<ν_(m), for a sample that contains the target         gas at the thermodynamic conditions (P, T)     -   2) Initialization procedure:         -   Local parameters initialization

[p _(lc) _(j) ^(ini)]_(j=1) ^(d*m) =[δB _(j) ^(ini) ,δA _(d) _(j) ^(ini),δγ_(c) _(j) ^(ini)]_(j=1) ^(d*m)=[0,0,0]  (Eq. 92)

-   -   -   Global parameters initialization: for p_(gb)(ν) two subset             of parameters are considered:         -   (a) The subset of free parameters p_(gb) ^(fr)

p _(gb) ^(fr) =[B,Ā _(d),γ _(c) ,d,[C _(l) ^(γ),γ_(l)(0),g _(l)(0)]_(l=0) ^(n) ]=p _(gb) ^(ini)  (Eq. 93)

-   -   -   -   which is initialized to an initial expected set of                 values p_(gb) ^(ini)

        -   (b) The subset of fixed parameters p_(gb) ^(cal)

p _(gb) ^(cal)(ν)=[Z _(d) ,[Z _(1Ul) ,g _(l)(δν_(il) ,w _(Rl))]_(l=0) ^(n)]^(cal)  (Eq.94)

-   -   -   -   which is initialized to the calibrated set of values                 p_(gb) ^(cal) determined following the procedure                 described below. They remaining unchanged across all fit                 procedure.

        -   Other initializations: L-M limit convergence values,             λ_(max),δχ_(max) ²,χ_(min) ² and maximum number of             interactions N_(itr)

    -   3) Minimization loop on local parameters p_(lc) _(j)         -   Keeping fixed p_(gb) ^(fr) during this step, a L-M             minimization loop is allowed on pt, by fitting each             experimental SCAR decay, {S_(exp)(t; ν_(i)),σ_(exp)(t;             ν_(i))}_(j) to the expected S(t; p(ν_(i))_(j)) function,             independently of the other detected decays. In this way, the             space of free parameters is reduced to 3 for each j fit,             easily to manage by L-M algorithms.         -   a) Keep fixed all global parameters to p_(gb)             ^(fx)(ν)=[p_(gb) ^(fr),p_(gb) ^(cal)(ν)]         -   b) Consider a trial set of local parameters p_(lc) _(j)             ^(try)=p_(lc) _(j) ^(ini) with δγ_(c) _(j) ^(try) calculated             by

δγ_(c) _(j) ^(try)=δγ_(c) _(j) ^(ini) +G _(j)(ν_(i) ;[C _(l) ^(γ)]_(l=0) ^(n))  (Eq. 95)

-   -   -   -   with G_(j)(ν_(i); [C_(l) ^(γ)]_(l=0) ^(n)) given by (Eq.                 85a). After the calculation of Eq. 95, the parameters                 [C_(l) ^(γ)]_(l=0) ^(n) are reset to zero value.

        -   c) Perform a L-M minimization routine in the set of             parameters p(ν_(i))_(j)=[p_(gb) ^(fx)(ν_(i)),p_(lc) _(j)             ^(try)] for the free local parameters p_(lc) _(j) ^(try),             i.e.             -   c1) Calculate the element S(t; p(ν_(i))_(j)), the χ²                 (p_(lc) _(j) ^(try))

$\begin{matrix} {{\chi^{2}\left( p_{{lc}_{j}}^{try} \right)} = {\sum\limits_{t_{k} = 0}^{t_{\max}}\frac{\left( {{S_{\exp}\left( {t_{k};v_{i}} \right)} - {s\left( {t_{k};{p\left( v_{i} \right)}_{j}} \right)}} \right)^{2}}{\sigma_{\exp}^{2}\left( {t_{k};v_{i}} \right)}}} & \left( {{Eq}.\mspace{14mu} 96} \right) \end{matrix}$

-   -   -   -   -   the gradient β(p_(lc) _(j) ^(try)), and the                     curvature matrix, |α|(p_(lc) _(j) ^(try))

            -   c2) Solve the linear Eqs. of the L-M algorithm for the                 present problem, and calculate δp_(lc) _(j) ,

χ²(p _(lc) _(j) ^(try) +δp _(lc) _(j) ) and δχ²=χ²(p _(lc) _(j) ^(try) +δp _(lc) _(j) )−χ²(p _(lc) _(j) ^(try))

-   -   -   -   -   If δχ²>0 then λ=k*λ and p_(lc) _(j) ^(try)=p_(lc)                     _(j) ^(try)                 -   If δχ²<0 then λ=λ/k and p_(lc) _(j) ^(try)=p_(lc)                     _(j) ^(try)+δp_(lc) _(j)

            -   c3) Follow the L-M convergence criteria on λ and/or χ²                 and/or δχ², i.e.                 -   If λ<λ_(max) and χ²>χ_(min) ² and δχ²>δχ_(max) ²,                     the solution does not converge, and come back to                     step c1) with an update value of λ and/or p_(lc)                     _(j) ^(try)=p_(lc) _(j) ^(try)+δp_(lc) _(j)                 -   If λ>λ_(max) and/or χ²<χ_(min) ² and/or δχ²>δχ_(max)                     ², the solution converges

            -   c4) At convergence, a fitted set of local parameters                 p_(lc) _(j) ^(fit) of the equal to the last guess p_(lc)                 _(j) ^(try) is determined. λ is reset to the initial                 value.

        -   d) Repeat steps b) to c) up to j=1 . . . d*m, i.e. the total             number of detected SCAR decays for the i=1 . . . m frequency             steps in the ν₁≤ν_(i)≤ν_(m) spectral range, repeated by d             times.

Therefore in this first step only the local parameters were allowed to vary, that is the minimization of χ² is made only on the local parameters, while the others are kept fixed.

Preferably, said second fitting step is a least squares fitting. More preferably, said least square fitting uses the Levenberg-Marquardt (L-M) algorithm.

This second fitting procedure is depicted in FIG. 7. In this second fit, the local parameters are set to the values found in the first fitting step, while the global parameters are free to vary in the second fitting procedure.

-   -   1) Minimization loop on free global parameters p_(gb) ^(fr)     -   a) Calculate the global weighted mean deviations of local         parameters P _(lc)=[δB,δA _(d),δγ _(c)]

$\begin{matrix} {\overset{\_}{\delta \; B} = \frac{\sum\limits_{j = 1}^{d*m}{{\sigma^{- 2}\left( {\delta \; B} \right)}_{j}^{fit}\delta \; B_{j}^{fit}}}{\sum\limits_{j = 1}^{d*m}{\sigma^{- 2}\left( {\delta \; B} \right)}_{j}^{fit}}} & \left( {{Eq}.\mspace{14mu} 97} \right) \\ {{\overset{\_}{\delta \; A}}_{d} = \frac{\sum\limits_{j = 1}^{d*m}{{\sigma^{- 2}\left( {\delta \; A_{d}} \right)}_{j}^{fit}\delta \; A_{d_{j}^{fit}}}}{\sum\limits_{j = 1}^{d*m}{\sigma^{- 2}\left( {\delta \; A_{d}} \right)}_{j}^{fit}}} & \; \\ {{\overset{\_}{\delta \; \gamma}}_{c} = \frac{\sum\limits_{j = 1}^{d*m}{{\sigma^{- 2}\left( {\delta \; \gamma_{c}} \right)}_{j}^{fit}\delta \; \gamma_{c_{j}^{fit}}}}{\sum\limits_{j = 1}^{d*m}{\sigma^{- 2}\left( {\delta \; \gamma_{c}} \right)}_{j}^{fit}}} & \; \end{matrix}$

-   -   -   with σ(p_(lc) ^(fit))=[σ(δB)_(j) ^(fit),σ(δA_(d))_(j)             ^(fit),σ(δγ_(c))_(j) ^(fit)]_(j=1) ^(d*m) the calculated             uncertainties of the fitted local parameters as described in             step 3) of the first fitting procedure.

    -   b) Update the fitted local parameters by:

δB _(j) ^(fit) =δB _(j) ^(fit)−δB

δA _(d) _(j) ^(fit) =δA _(d) _(j) ^(fit)−δA _(d)

δγ_(c) _(j) ^(fit)=δγ_(c) _(j) ^(fit)−δγ _(c)  (Eq. 98)

-   -   -   and keep them fixed during the rest of the step 4)

    -   c) Considers a trial of free global parameters p_(gb)         ^(fr)=p_(gb) ^(try)==p_(gb) ^(ini)

    -   d) Perform a L-M minimization routine in the set of parameters

p(v) = [p_(gb)^(try), p_(gb)^(cal)(v), [p_(lc_(j))^(fit)]_(j = 1)^(d * m)]

-   -   for the free global parameters p_(gb) ^(try), i.e.         -   d1) Calculate the elements [S(t; p(ν_(i))_(j))]_(j=1)             ^(d*m), the χ² (p_(gb) ^(try))

$\begin{matrix} {{\chi^{2}\left( p_{gb}^{try} \right)} = {\sum\limits_{j = 1}^{d*m}{\sum\limits_{t_{k} = 0}^{t_{\max}}\frac{\left( {{S_{\exp}\left( {t_{k};v_{i}} \right)} - {S\left( {t_{k};{p\left( v_{i} \right)}_{j}} \right)}} \right)^{2}}{\sigma_{\exp}^{2}\left( {t_{k};v_{i}} \right)}}}} & \left( {{Eq}.\mspace{14mu} 99} \right) \end{matrix}$

-   -   -   -   the gradient β(p_(gb) ^(try)), and the curvature matrix,                 |α|(p_(gb) ^(try))

        -   d2) Solve the linear Eqs. of the L-M algorithm for the             present problem, and calculate δp_(gb),

χ²(p _(gb) ^(try) +δp _(gb)) and δχ²=χ²(p _(gb) ^(try) +δp _(gb))−χ²(p _(gb) ^(try))

-   -   -   -   If δχ²>0 then λ=k*λ and p_(gb) ^(try)=p_(gb) ^(try)             -   If δχ²<0 then λ=λ/k and p_(gb) ^(try)=p_(gb)                 ^(try)+δp_(gb)

        -   d3) Follow the L-M convergence criteria on λ and/or χ²             and/or δχ², i.e.             -   If λ<λ_(max) and χ²>χ_(min) ² and δχ²>δχ_(max) ², the                 solution does not converge, and come back to step 3) of                 this recipe with an updated value of λ and/or p_(gb)                 ^(try)=p_(gb) ^(try)+δp_(gb), and [p_(lc) _(j)                 ^(try)=p_(lc) _(j) ^(fit)]_(j=1) ^(d*m)             -   If λ<λ_(max) and χ²<χ_(min) ² and/or δχ²<δχ_(max) ², try                 a new minimization loop from the step 3) up to 4)d2)                 with an updated value of λ and/or p_(gb) ^(try)=p_(gb)                 ^(try)+δp_(gb), and [p_(lc) _(j) ^(try)=p_(lc) _(j)                 ^(try)]_(j=1) ^(d*m), to verify if this is the best fit                 of all set of fitted parameters. After this loop go to                 step 4)e).

    -   e) If χ²<χ_(min) ² and/or δχ²<δχ_(max) ² and/or λ>λ_(max) the         final convergence of the fit is reached and a set of the fitted         parameters

[p_(gb)^(fit), [p_(lc_(j))^(fit)]_(j = 1)^(d * m)]

-   -   of the equal to the last guess

[p_(gb)^(try), [p_(lc_(j))^(try)]_(j = 1)^(d * m)]

-   -   is determined.     -   2) Calculate the concentration of the target gas by

$\begin{matrix} {N_{0} = \frac{{\gamma_{0}^{fit}(0)}{g_{0}^{fit}(0)}}{c^{2}{L_{S}(T)}}} & \left( {{Eq}.\mspace{14mu} 100} \right) \end{matrix}$

-   -   -   Using the results of the first and second fitting step.

Advantageously, the method of the invention includes the step of:

-   -   Dividing the second group of parameter in a first sub-group of         free parameters and a second sub-group of calibration parameters         which are set equal to calibration values which remain constant         during the first and second fitting steps.

As shown in the explanation of an embodiment of the first fitting step above, in particular in (Eq. 94), not all parameters denominated “global parameters” (that is, the second group of parameters) are freely varying during the second step of the fitting method. Some of them are equal to a value which is calculated and fixed before the fitting routine.

More preferably, the sub-group of calibration parameters includes:

p _(gb) ^(cal)(ν)=[Z _(d) ,[Z _(1U1) ,g _(l)(δν_(il) ,w _(Rl))]_(l=0) ^(n)]^(cal)

-   -   Where {₁(ν_(i)−ν_(l),w_(Rl))}_(l=0) ^(n) are the peak normalized         line profiles for each transition l=0 . . . n, centered at the         resonance frequency ν_(l) and linewidth w_(Rl), and     -   Z_(d) is the “equivalent”-saturation parameter at 1U due to         residual detection non-linearities, in U⁻¹.

Calibration routines of the fixed set of parameters, p_(gb) ^(cal) (see Eq. 94), used in the second aspect of the method of the invention for the analysis of the SCAR signals are preferably implemented. This calibration is preferably performed for each SCAR apparatus for the same thermodynamic conditions (Pressure, Temperature) of the sample mixture used in the trace gas detection to which the method according to the second aspect will be applied.

A first calibration routine is relative to the calibration of

[Z _(1U0) ,g ₀(δν_(i0) ,w _(R0))]^(cal)

The first calibration concerns the saturation parameter at 1U, Z_(1U0), and absorption lineshape g ₀(δν_(i0),w_(R0)) of the transition of the target gas (l=0) that absorb cavity-resonant light in the recorded spectral range.

The experimental decay curves used in this calibration routine are recorded with a sample gas mixture enriched with the target gas to be calibrated, that is in a condition in which the concentration of the target gas in the cavity is equal to a second concentration much higher than the first concentration. This enrichment must be enough to measure a spectrum γ₀ (ν) with high S/N for the target transition. The thermodynamic conditions (P,T) of these enriched sample mixture are the same of the sample mixture used in the trace gas detection (that is, when the first concentration is calculated). This approach to consider a high S/N spectrum for the target transition allows on one side to determine the calibrated parameters for this transition with high precision and accuracy and on the other side to consider negligible the contribution of the other k=1, . . . , n transitions and detected residual non-linearities (i.e d=0).

The recipe, schematically depicted in FIG. 8, is the following:

-   -   1) Consider all the SCAR curves (j=1, . . . , d*m) detected in         the spectral range ν₁≤ν_(i)≤ν_(m) for the enriched sample for         the target gas {S_(exp) ^(↑S/N)(t, ν),σ_(exp) ^(↑S/N)(t;         ν)}=[{S_(exp) ^(↑S/N)(t; ν_(i)),σ_(exp) ^(↑S/N)(t;         ν_(i))}_(j)]_(j=1) ^(d*m)     -   2) Fit this data by using the method of the invention according         to the second aspect described in the section above with the set         of parameters p(ν)=p^(↑S/N)(ν):

p ^(↑S/N)(ν)=[p _(gb) ^(↑S/N)(ν),[p _(lc) _(j) ]_(j=1) ^(d*m)]  (Eq.101)

-   -   -   with

p _(gb) ^(↑S/N)(ν)=[ B,Ā _(d),γ _(c) ,C ₀ ^(γ) ,Z _(1U0),γ₀ ^(↑S/N)(0),g ₀ ^(↑S/N)(0), g ₀(ν_(i0) ,w _(R0))]  (Eq.102)

[p _(lc) _(j) ]_(j=1) ^(d*m) =[δB _(j) ,δA _(d),δγ_(c) _(j) ]_(j=1) ^(d*m)

-   -   -   i.e. keeping free all spectral parameters that define the             target transition and fixed to zero the spectral parameters             of the other lines and d=0.

    -   3) At the convergence of the above step, the calibrated values         p_(gb) ₀ ^(cal)=[Z_(1U0),g ₀(δν_(i0),w_(R0))]^(cal) equal to the         fitted values [Z_(1U0),g ₀(δν_(i0),w_(R0))]^(fit) are         determined.

A second calibration routine is relative to

[Z _(d) ,[Z _(1Ul) ,g _(l)(δν_(il) ,w _(Rl))]_(l=1) ^(n)]^(cal)

In this case the routine concerns the calibration of the saturation parameters and line-shapes of all other absorbing transitions different from the target one (l=1 . . . n), as well as the parameters that define the residual non-linearity, Z_(d).

Differently from the first calibration routine, this routine uses a gas sample mixture totally depleted of the target gas at the same Pressure and Temperature conditions of those used in the trace gas experiment. In this situation there is no absorption from the target gas for the detected SCAR signals,

{S_(exp)^( ↓ S/N)(t; v), σ_(exp)^( ↓ S/N)(t; v)} = [{S_(exp)^( ↓ S/N)(t; v_(i)), σ_(exp)^( ↓ S/N)(t; v_(i))}_(j)]_(j = 1)^(d * m),

and the parameters related with the target gas (l=0) must be considered negligible.

The recipe, schematically depicted in FIG. 9, is the following:

-   -   1) Consider all the SCAR curves (j=1, . . . , d*m) detected in         the spectral range ν₁≤ν_(i)≤ν_(m) for the totally depleted         sample for the target gas,

{S _(exp) ^(↓S/N)(t;ν),σ_(exp) ^(↓S/N)(t;ν)}=[{S _(exp) ^(↓S/N)(t;ν _(i)),σ_(exp) ^(↓S/N)(t;ν _(i))}_(j)]_(j=1) ^(d*m)

-   -   2) Fit this data by using the method of the invention according         to the second aspect described in the section above with the set         of parameters p(ν)=p^(↓S/N)(ν):

p ^(↓S/N)(ν)=[p _(gb) ^(↓S/N)(ν),[p _(lc) _(j) ]_(j=1) ^(d*m)]  (Eq.103)

with

p _(gb) ^(↓S/N)(ν)=[ B,Ā _(d),γ _(c) ,d,Z _(c) ,[C _(l) ^(γ) ,Z _(1Ul),γ_(l) ^(↓S/N)(0),g _(l) ^(↓S/N)(0), g _(l)(δν_(il) ,w _(Rl))]_(l=1) ^(n)]  (Eq.104)

[p _(lc) _(j) ]_(j=1) ^(d*m) =[δB _(j) ,δA _(d) _(j) ,δγ_(c) _(j) ]_(j=1) ^(d*m)

-   -   3) At the convergence of the above step, the calibrated values         [p_(gb) _(l) ^(cal)]_(l=1) ^(cal)=[Z_(d),[Z_(1Ul),g         _(l)((δν_(il),w_(Rl))]_(l=1) ^(n)]^(cal) equal to the fitted         values [Z_(d),[Z_(1Ul),g _(l)(δν_(il),w_(Rl))]_(l=1) ^(n)]^(fit)         are determined.

BRIEF DESCRIPTION OF THE DRAWINGS

The invention will be now described in a non-limiting manner with reference to the appended drawing, in which:

FIGS. 1-5 are different method steps of the method of the invention according to its first aspect;

FIGS. 6-9 are different method steps of the method of the invention according to its second aspect;

FIG. 10 is a schematic diagram of an apparatus to perform the method of the invention according to the first or to the second aspect;

FIG. 11 is a first measurement performed and fitted according to the first aspect of the method of the invention;

FIG. 12 is a fitting of a measurement similar to the one in FIG. 11, but with a larger frequency scan, in CRD conditions performed according to a prior art method;

FIG. 13 is a fitting of the same measurement of the one of FIG. 12 in SCAR conditions but according to a method not of the invention;

FIG. 14 is the fitting of the same measurement of the one of FIG. 12 with a fitting according to the invention;

FIG. 15 is a second measurement performed and fitted according to the first aspect of the method of the invention;

FIG. 16 is a fitting of the same measurement of the one in FIG. 11 in CRD conditions performed according to a prior art method;

FIG. 17 is a fitting of the same measurement of the one of FIG. 11 in SCAR conditions but according a method not of the invention;

FIG. 18 is the fitting of FIG. 11 with a fitting according to the invention;

FIG. 19 is a third measurement performed on a highly enriched sample and fitted according to the second aspect of the method of the invention; and

FIG. 20 is a fitting of the same measurement of the one in FIG. 19 in CRD conditions performed according to a prior art method.

DETAILED DESCRIPTION OF AN EMBODIMENT

With initial reference to FIG. 10, reference numeral 1 indicates an apparatus for SCAR spectroscopy including a laser source 2, for example a continuous wave (CW) coherent laser source generated by a frequency difference tunable over a predetermined range.

Preferably, the radiation emitted by laser source 2 has a wavelength in the mid-infrared, however other wavelengths may be used. The mid-infrared has the advantage of having the strongest molecular absorption.

The type of laser source 2 used in the present invention is for example described in the article written by Galli et al., Opt. Lett. 35, 3616 (2010). Other types of laser sources may be used, provided that the intensity of radiation I inside the cavity is much greater than the intensity of saturation Is of the molecular transition to be detected, i.e. I>>I_(s).

For example, in a range of wavelengths of 4-5 μm, for the transitions of CO₂, which have an Einstein coefficient A of about 200 s⁻¹, at a pressure of about 12 mbar, Voigt enlargement condition, the power emitted by the laser must be greater than 20 mW, preferably greater than 100 mW.

Device 1 further includes a resonant cavity 3, for example, a cavity having a length of 1 m, provided at opposite ends thereof with two reflecting mirrors 4 a and 4 b. Preferably, the reflectivity of the mirrors is greater than 99.9%, even more preferably it is greater than 99.99%.

The gas of which the concentration has to be measured is introduced into cavity 3, for example through a duct 5 which connects cavity 3 to a suitable container, such as a cylinder 6.

Apparatus 1 further includes a photodetector 7 suitably arranged for detecting the radiation beam outgoing from cavity 3 as well as a diffuser element 8 interposed between cavity 3 and photodetector 7.

The diffuser element 8 is adapted to diffuse the laser beam exiting cavity 3 before it impinges on photodetector 7.

Experiment 1

The SCAR apparatus is composed of an IR laser emitting around 4.5 μm, where ¹⁴CO₂ absorbs, a high-finesse Fabry-Perot cavity and an IR-detector followed by amplifying and digitalizing electronics.

The IR radiation is provided by an Optical-Frequency-Comb (OFC) assisted Difference-Frequency-Generation (DFG) continuous-wave (CW) coherent source (“Ti:sapphire laser intracavity difference-frequency generation of 30 mW cw radiation around 4.5 μm”, I. Galli et al, Opt. Lett. 35, 3616 (2010)). The DFG process occurs inside the cavity of a Ti:Sapphire laser operating around 850 nm (pump laser), single-mode controlled by an injected extended-cavity diode-laser (ECDL). A Nd:YAG laser at 1064 nm, amplified up to 10 W using an Yb-doped fiber amplifier provides the DFG signal laser. It is mixed with the intracavity Ti:Sapphire radiation through a periodically-poled lithium niobate non-linear crystal. The frequency of the ECDL is phase-locked to the Nd:YAG frequency by direct digital synthesis, using the OFC to cover the frequency gap (about 70 THz) between the two CW lasers (“Ultra-stable, widely tunable and absolutely linked mid-IR coherent source”, I. Galli et al., Opt. Express 17, 9582 (2009)). In this way, the linewidth of the IR generated radiation is given by a fraction of the narrow Nd:YAG linewidth (about 5 kHz in 1 ms integration time), thus allowing a highly efficient coupling of the IR radiation to the high-finesse Fabry-Perot cavity. Moreover, the frequency chain used to cover the frequency difference between both signal and pump lasers of the DFG process, includes a microwave synthesizer. In this way, frequency of the generated IR radiation is scanned and tuned in a synthesized way by changing the microwave frequency. In addition, the Nd:YAG frequency is stabilized against the nearest tooth of the OFC. As a consequence, the IR frequency is absolutely traceable against the primary frequency standard with a precision of 6×10⁻¹³ in 1 s and an accuracy of 2×10⁻¹². Moreover, the intracavity DFG boosts the generated IR power up to 30 mW around 4.5 μm wavelength, which provides the required power to saturate the ¹⁴CO₂ transitions.

The measurement cell is a cylindrical vacuum chamber 1.2 m long and 10 cm in internal diameter. It is enclosed in a polystyrene box, which can be filled with dry-ice pellets for cooling the chamber down to 195 K. The chamber houses a Fabry-Perot optical resonator resting inside of it on 4 cantilevered legs, which dampen vibrations (for frequencies >20 Hz) in all 3 spatial directions. The mechanical frame of the Fabry-Perot resonator is made of 3 Invar bars connected by 2 circular flanges, leaving 8 L internal volume available for the gas. At both ends of the frame, properly machined flanges house the mirrors, with high reflectivity dielectric coatings at 4.5 μm wavelength and a 6-m of radius of curvature. The total losses for each mirror (transmission plus absorption/scattering) amount to ^(˜)270 parts per million (ppm) and the achieved optical finesse is higher than 11,000. The mirror mounting flanges have screws for coarse alignment and PZT for fine adjustment of both the alignment and the cavity length. The mirror spacing is 1 m and the corresponding free spectral range is 150 MHz.

A N₂-cooled InSb detector is used to detect the radiation transmitted by the cavity. It is followed by a transimpedance amplifier (Z=32000 Ohm) with a final bandwidth of about 1 MHz. A 18-bit digitizing oscilloscope with a sampling rate of 10 Ms/s is used to digital convert the analogical detected signal for further process and analysis.

For the SCAR experiment, the IR radiation is efficiently coupled (about 86% as expected) in the high-finesse Fabry-Perot cavity by using a couple of lenses to achieve the TEM₀₀ mode propagation inside of the cavity. In this conditions, taking into account the mirror losses and an available power of about 20 mW before the cavity input, we estimate an inside cavity power in resonance of about 40 W.

Then, the frequency of the IR coherent source is changed to be in resonance with the cavity. OFC-assisted frequency stability of the IR radiation and the high mechanical stability of the cavity length allow to maintain a long term resonance condition without need of any active lock of the cavity length to the IR frequency of the laser.

A double-pass acusto-optic modulator placed at the output of Nd:YAG laser is used to quickly switch off resonance the input IR light when a threshold coupling level (about 3V at the output of the amplified detector) is reached, and thus the transmission cavity decay is detected by a N2-cooled InSb detector during 100 μs. The output tension of the amplified detector is digitalized, and thus recording the single SCAR-decay-event. In order to improve S/N ratio of the detected decay, compatible with digitalization resolution, 128 consecutive SCAR-decay-events are averaged, giving the SCAR-decay-signals to be analyzed by the fit routines, as described in this invention.

For the SCAR spectroscopy of ¹⁴C¹⁶O₂, the (00⁰1-00⁰0) P(20) rovibrational transition around 4.5 μm was targeted. It is a quasi-isolated absorption line of this molecule with respect to other absorptions from other CO₂ isotopologue. Nevertheless, the (05⁵1-05⁵0) P(19)e line of ¹³C¹⁶O₂, blue-shifted by ^(˜)230 MHz with respect to the frequency of the target transition, is an interference line that must be taken into account in the analysis to get an accurate concentration measurement of radiocarbon dioxide.

The Fabry-Perot cell is filled with a CO₂ sample at total pressure P=12 mbar and temperature T=195 K (corresponding to ^(˜)0.15 L volume at standard thermodynamic conditions). At this temperature, the intensity of the nearby ¹³C¹⁶O₂ interference line is decreased by more than three orders of magnitude, minimizing its interference effect. The CO₂ sample is a natural mix of almost all carbon dioxide isotopologues including ¹⁴C¹⁶O₂, at the present natural abundance (about 1.2×10⁻¹²). Other gases must be avoided and in particular ¹⁴N₂ ¹⁶O, which has absorption line almost at the same frequency of the target transition. N₂O concentration below 0.3 ppb in the measured CO₂ sample is required to produce a negligible interference with the target line.

In these conditions of temperature and pressure, the deformation from a pure exponential produced by radiocarbon dioxide absorption along decay signals is of the order of 1 μV out of 3 V. This set very stringent limits on the residual non-linearity, which can be born. At this respect particular attention was putted to get a linear response of the detection system for the more than 6 decades as require radiocarbon-dioxide detection at ppt level or better. An optical diffuser was placed between the cavity and the detector, in order to uniform the way to illuminate the detector area. In addition, a numerical calibration of residual non-linearities is allowed in the SCAR decays obtained in vacuum conditions. In this case, the decay must be a perfect exponential, and deviations are taken into account by adding a Fourier function to the expected exponential behavior. This correction is then applied to all recorded SCAR decays.

On the other side, the measurement is possible thanks to improve the noise present in the signal to be captured, which, in good approximation, has zero, or at least constant, average. In addition to the 128 consecutive SCAR decays averaged by the digitation electronics as described above, ten of these acquisitions for each laser frequency are again averaged for a total of 1280 consecutive decay events. The decay resulting of the total average, which is stored in memory for further analysis as described in this invention, increases the resolution of digitization by approximately 35 times.

To obtain the SCAR decay signals for ¹⁴CO₂ concentration measurements, the following steps have been performed:

-   -   1) The cavity-cell is filled with the CO₂ sample at the         thermodynamic conditions described above. The IR laser is         operated at the required power for saturation of the target         transition (see above) and coupled to the cavity as described         above.     -   2) The frequency of the IR coherent source is set to a value 66         227 000 MHz, which is 400 MHz apart of the known frequency of         the target P(20) transition of radiocarbon dioxide. This is         achieved by setting the proper OFC repetition rate and taking         into account the lock frequency chain described above.     -   3) The Fabry-Perot is bringing in resonance with the IR         radiation by changing the cavity length by means of the tension         applied to the PZT transducers. The transmitted intensity is         detected with the N₂-cooled InSb detector.     -   4) When the transimpedance tension at the output of the         amplified detector reach the 3 V threshold condition, the laser         frequency is quickly and automatically switched-off the cavity         resonance by means of the AOM as described above.     -   5) The SCAR decay is detected by 100 μs, digitized and averaged         with consecutive decay signals. The averaged signal is         temporally stored in the digitizing electronics.     -   6) The laser is quickly and automatically switched-on the cavity         resonance by means of the AOM as described above, waiting to         threshold condition be again reached.     -   7) The steps 4 to 6 are repeated by 128 consecutive decay         events. This step is repeated by ten times, where the result of         each 128-average decay is averaged with the consecutive one. At         the end, a SCAR decay signal for the present frequency of the         spectrum (i.e. the laser frequency), resulting of the average of         1280 consecutive events is stored for the further analysis         described in this invention.     -   8) The laser frequency is changed-up by a frequency step of 10         MHz, and the steps from 3 to 7 are repeated up to the laser         frequency is scanned-up by 740 MHz. (i.e. 73 frequency steps in         the upward direction).     -   9) The laser frequency is changed-down by a frequency step of 10         MHz, and the steps from 3 to 7 are repeated up to the laser         frequency is scanned-down by 740 MHz. (i.e. 73 frequency steps         in the downward direction).     -   10) As a final result, 146 SCAR decay signals, two for each         scanned IR frequency, are stored to be analyzed by the methods,         subject of this invention.

The results of this experiment and the fit according to the first aspect of the invention can be seen in FIG. 11. A spectrum of the ¹⁴C¹⁶O₂ P(20) transition (black trace) recorded at the present natural abundance (named “2010”) is shown. In this figure, two more spectra are shown, corresponding to a ¹⁴C-enriched sample (green trace) and a ¹⁴C-depleted sample (blue trace). ¹⁴C-enriched spectra were used to accurately measure the center frequency of the P(20) transition, as well as the most favorable thermodynamic conditions of the gas sample for ¹⁴C detection in subnatural abundance.

More quantitatively, the fit of the 2010 CO₂ spectrum to the expected Voigt profile yields an area of 1:50(8) ms⁻¹ MHz. Assuming for the P(20) the calculated line strength, S=3.10(15)*10⁻¹⁸ cm, a ¹⁴C¹⁶O₂ natural abundance concentration of 1.24(10) ppt is measured.

For the calibration it has been used the target gas at a second concentration of: 58 times the first concentration. Further, in the calibration the target gas was at a temperature: 195K, and pressure 11.6 mbar. The parameter obtained from the calibration are:

line center: 66227382.3 MHz FWHM Gaussian width: 123.54 MHz FWHM Lorentzian width: 74.32 MHz Z_(1Ueff)=8.2 V⁻¹ Ratio R of the spectral area: (effective area)=(area measured by SCAR)×0.925.

FIG. 12 shows a similar measurement, but with a larger frequency scan, performed with a standard CRD routine, that is, not in saturation condition, and a fitting of the result. It is clear that the variations (drift and oscillation) of the curve for a standard CRD measurement is much higher. The FIG. 12 shows two different measurements which are partially overlapped, however they do not match.

FIG. 13 shows SCAR measurements of the same sample in which, instead of using the method according to the first aspect, the line profile has been selected a priori. The tails of the curve as shown are not correct and they move upwards.

FIG. 14 shows the SCAR fitting method according to the first aspect of the invention, which has a much more correct behaviour.

Experiment 2

The experimental set up is the same as in experiment 1.

Sub-Doppler measurements on low-pressure ¹⁷O¹²C¹⁶O at natural abundance (7.5*10⁻⁴) have been performed. FIG. 15 shows the resolved hyperfine structure of the (00⁰1-00⁰0) R(0) transition of ¹⁷O¹²C¹⁶O, due to interaction between the ¹⁷O electric quadrupole (which is non-null for a nuclear spin I=5/2) and the electric field gradient at the nucleus position. This spectrum was recorded in about 3 h with 11 forward-backward frequency scans in 20-kHz steps.

In Table I the fit results for the line centers are reported.

TABLE I Measured absolute frequencies ν_(ΔF) of the hyperfine triplet with 1−σ uncertainties and relative intensities. ΔF ν_(ΔF) (kHz) rel. intens. −1 70 174 358 594.9 (6.2) 0.207 (11) 0 70 174 357 409.8 (3.0) 0.345 (9)  +1 70 174 358 229.8 (3.6) 0.448 (13)

The fitted FWHM of each Lorentzian dip is Γ=217.1(5:5) kHz.

The measured line-center frequency ν_(c)=70 174 358 037.3 (3:9) kHz is consistent with the value 70 174 368 MHz reported by HITRAN, within the declared 3-30 MHz uncertainty. Thus, we improved the frequency accuracy by more than 3 orders of magnitude.

FIG. 16 shows for comparison the same measurement of the one in FIG. 11 according to CRD, that is in a linear regime, and with its fitting.

FIG. 17 shows SCAR fitting of the same measurement of the one of FIG. 11, in which, instead of using the method according to the first aspect, the line profile has been imposed a priori. The tails of the curve as shown are not correct and they move upwards.

FIG. 18 shows the SCAR fitting method according to the first aspect of the invention, which has a much more correct behaviour.

Experiment 3

In this experiment, the fitting is performed according to the second aspect of the invention.

The apparatus used is the one described according to the first experiment.

SCAR spectroscopy in a ¹⁴C¹⁶O²-enriched CO₂ sample as in the first experiment. The laser frequency was scanned across the P(20) line, and the SCAR decay signals where fitted by using the above described fit function.

The sample was much more enriched than in the first experiment (about 6400 times with respect to natural abundance) in order to minimize the effects of the S/N ratio in the fitted parameters uncertainty and to calibrate the SCAR procedure. Moreover, to check the linearity of the detection, this very high radiocarbon abundance has been changed by reducing it with successive dilutions. Another set of measurements was done with gas at natural abundance, obtained by fermentation of the same cane sugar. All these new measurements were done at temperature T=195 K and pressure P=12 mbar.

In this procedure, at the highest abundance, the line center, the Lorentzian and Gaussian width, the spectral area cα_(o) of the line profile and the saturation parameter Z_(1U) were fitted as “global” variables, while the amplitude, the background and the cavity decay constant of each single decay signals were fitted as “local” variables.

The values obtained for Z_(1U), line widths and line center at the highest abundance where then used as fixed parameter for the fit of subsequent measured signals at lower abundances.

The decay signals were detected at various laser frequencies, covering a total span of about 650 MHz around the P(20) line frequency, at 10 MHz steps. Due to technical reasons of our experimental setup, continuous scans larger than 400 MHz are not allowed, and such 650 MHz span is the result of the partially overlapping of two consecutive 400 MHz scans at the blue and red sides of the line center. At this frequency span, the interfering line (05⁵1-05⁵0) P(19) of ¹³C¹⁶O₂ on the “blue” side of the P(20) line must be considered. Its interfering effects at 195 K may be neglected for highly-enriched samples, nevertheless they have been considered.

The interfering line parameters ν′_(o), w′_(L), w′_(G) and Z′_(1U) were fixed to a value determined with the following procedure: first preliminary values for the global parameters have been calculated by using the enriched-sample SCAR spectrum without considering the interfering line. Then, a SCAR spectrum recorded with a CO₂ sample with radiocarbon dioxide at natural abundance has been fitted by fixing all “global” parameters except the P(20) line area, d, and all interfering line parameters. Indeed, at this ¹⁴C¹⁶O₂ abundance and the values of T=195K and P=12 mbar, the area of the interference line is almost equal to the P(20) one, and hence its effects are better taken into account. Therefore, with these first estimations of the parameters of the interfering line, the global parameters of the enriched spectrum are again calculated. This procedure is repeated several times, till final convergence.

In FIG. 19, one of the three SCAR spectra recorded for the about 6400 times enriched sample is shown. The fitted global parameters for the three spectra were all self-consistent within the errors. The averaged values are shown in Table 2.

These results have been checked by comparing them with those measured for the spectra recorded with the same CO₂ enriched sample at the same thermodynamic conditions, but near the linear-absorption regime.

The laser power was not as low as to decrease too much the S/N ratio. In order to fit the decay signals, we used the same global procedure. However, because, at low power, the two decay parameters γ_(c) and γ_(g) are quite correlated, the signal is essentially obtained by reducing the SCAR procedure to the standard CRD procedure, where the overall γ=γ_(c)+γ_(g) is considered, and by subtracting an averaged value for the cavity decay rate out of resonance. In this way, the non-linearity due to a residual saturation was corrected and an experimental linear absorption line profile was generated. A successive Voigt fit also takes into account, in this case, optical fringes that modulate the γ_(c) contribution.

FIG. 20 shows the average of three profiles obtained with successive measurements and the corresponding fit. The residuals do not allow to shed light on the real shape of the profile, alternative to the Voigt profile. Indeed, the residuals obtained in the left and right scans have different values in the common frequency range: likely they are dominated by an incomplete neutralization of the fringes. The main discrepancy from a Voigt profile is an asymmetric residual in the wings. Also in this regime the values of the fitted Voigt profile obtained from the three independent measurements are consistent within the errors.

The averaged values are shown in Table 3.

These values do not coincide with those obtained at high power. For the linewidths, the discrepancies can be ascribed to the local approximation of the laser field-molecules interaction done in the theory. Indeed, the residuals for the saturated case (FIG. 19) clearly show a systematic small discrepancy from a Voigt profile. Anyway, such discrepancies are small and, for the calibration purpose, its effect will be included in the ratio between the areas estimated in the two regimes. The ratio between the area obtained in the linear regime and the area obtained at high power with the SCAR procedure is 1.0284(8).

Again, this is the correction factor to get the “true” radiocarbon dioxide concentration from the SCAR measurements in the cell at our experimental conditions with this new analysis procedure. Its discrepancy can be considered as a measure of the goodness of the present approach.

TABLE 2 Fitted parameters of the SCAR spectrum of the P(20) transition of ¹⁴C¹⁶O₂ in a 6375 times enriched CO₂ sample (p = 12 mbar, T = 195 K) Line center, ν_(o) - 66227000 MHz 382.1143(47) MHz Spectral Area 8968.44(32) ms⁻¹MHz Lorentzian HWHM, w_(L) 49.927(25) MHz Gaussian HWHM, w_(G) 41.6905(42) MHz Saturation parameter Z_(1V) 9.030(24) V⁻¹

TABLE 3 Fitted parameters of the linear-CRD spectrum of the P(20) transition of ¹⁴C¹⁶O₂ in a 6375 times enriched CO₂ sample (p = 12 mbar, T = 195 K). Line center, ν_(o)- 66227000 MHz 382.34 (7)  MHz Spectral Area  9224 (7) ms⁻¹MHz Lorentzian HWHM, w_(L) 53.30 (9) MHz Gaussian HWHM, w_(G)  39.22 (13) MHz 

1. A method of ring-down spectroscopy in saturated-absorption condition, for measuring a first concentration of a gas through a measurement of the spectrum of a molecular transition of said gas, the method comprising the steps of: inserting said gas whose first concentration is to be measured in a resonant cavity comprising two or more reflecting mirrors arranged so as to form a closed optical path for an electromagnetic radiation emitted by a laser source; tuning the frequency of said electromagnetic radiation emitted by said laser source so as to fix it to a value ν_(i) within a range of frequencies [ν_(min), ν_(max)] including the resonance frequency of said molecular transition ν₀; fixing the intensity of said electromagnetic radiation in the cavity at a value much greater than the saturation intensity I_(s) of the molecular transition to be detected; irradiating said gas by means of said electromagnetic radiation beam emitted by said laser source having said fixed frequency ν_(i) and intensity in said resonant cavity; coupling said electromagnetic radiation to said cavity so as to obtain a laser-cavity resonance condition; changing the frequency of the electromagnetic radiation emitted by the laser so as to switch off the laser-cavity resonance; detecting an electromagnetic radiation beam in output from said cavity after the laser-cavity resonance has been switched off; recording a plurality of data representative of said output, obtaining a decay signal for the fixed frequency; considering a fitting curve S(t, δν_(i)) for the recorded decay signal which depends on the following parameters: B(δν_(i)) is the detection background, with (δν_(i)=ν_(i)−ν₀); A_(d)(δν_(i)) is the amplitude of the decay signal at the beginning of the decay event; γ_(c)(δν_(i)) is the cavity decay rate due to non-resonant and non-saturable losses (empty cavity decay rate); γ_(g)(δν_(i)) is contribution of the targeted molecular transition to the decay signal; g is the peak normalized line profile g(ν−ν₀,w_(R)) centered at the molecular resonance frequency ν₀ and w_(R) is the HWHM width of the resonance, and w_(R)=w_(L) for a Lorentzian shape, w_(R)=w_(G) for a Gaussian shape, w_(R)={w_(L), w_(G)} for a Voigt shape; $Z_{o} = {\frac{{CP}(0)}{{CP}_{s}} = {A_{d}Z_{1U}}}$ is the saturation parameter at the beginning of the decay event and at the frequency of the targeted molecular transition ν₀, P(0) is the intracavity power at the beginning of the decay signal; and $P_{s} = {\frac{\pi \; w^{2}}{2}I_{s}}$ is the saturation power, where w is the spot size radius of the laser beam, i.e. the radius for which the amplitude of the field is 1/e times that of the axis and I_(s) the saturation intensity; replacing Z_(1U) g(ν−ν₀,w_(R)) in the function S(t, δν_(i)) with a constant value Z_(1ueff)=constant of a predetermined value; fitting said recorded data with a function S^(repl)(t, δν_(i)) in which Z_(1Ueff)=constant replaces Z_(1U) g(ν−ν_(o),w_(R)) in the fitting function S(t, δν_(i)).
 2. The method according to claim 1, including the step of: obtaining from said fit the following parameters: B(δν_(i)) is the detection background, with (δν_(i)=ν_(i)−ν₀); A_(d)(δν_(i)) is the amplitude of the decay signal at the beginning of the decay event; γ_(c)(δν_(i)) is the cavity decay rate due to non-resonant and non-saturable losses (empty cavity decay rate); γ_(g)(δν_(i)) is contribution of the targeted molecular transition to the decay signal.
 3. The method according to claim 1, including the step of parametrizing said S^(repl)(t, ν_(i)) as: S ^(repl)(t;p(δν_(i)))=B(δν_(i))+A _(d)(δν_(i))e ^(−γ) ^(c) ^((δν) ^(i) ^()t)ƒ(t;γ _(c)(δν_(i)),A _(d)(δν_(i)),γ_(g)(δν_(i)),Z _(1Ueff)) where p(δν_(i)) is the set of parameters free during the fit p(δν_(i))={B(δν_(i)),A _(d)(Sδν _(i)),γ_(c)(δν_(i)),γ_(g)(δν_(i))} B(δν_(i)) is the detection background, γ_(g)(δν_(i)) is contribution of the targeted molecular transition to the decay signal, A_(d)(δν_(i)) is the amplitude of the decay signal, γ_(c)(δν_(i)) is the cavity decay rate due to non-resonant and non-saturable losses, while the following are determined before the fitting: Z_(1Ueff) is the effective saturation parameter, fixed during the fit and equal to a constant; ƒ is the non-linear function that follows one of the below rate equations depending of the gas conditions: homogeneous regime (w_(L)≥w_(G)) $\frac{df}{dt} = {{- {\gamma_{g}\left( {\delta \; v_{i}} \right)}}\frac{\ln \left\lbrack {1 + {{A_{d}\left( {\delta \; v_{i}} \right)}Z_{1\; {Ueff}}e^{{- {\gamma_{c}{({\delta \; v_{i}})}}}t}f}} \right\rbrack}{{A_{d}\left( {\delta \; v_{i}} \right)}Z_{1{Ueff}}}}$ inhomogeneous regime still diffusive gas (w_(L)<w_(G)) $\frac{df}{dt} = {{- {\gamma_{g}\left( {\delta \; v_{i}} \right)}}\frac{2}{\sqrt{1 + {{A_{d}\left( {\delta \; v_{i}} \right)}Z_{1{Ueff}}e^{{- {\gamma_{c}{({\delta \; v_{i}})}}}t}f}} + 1}f}$ non-diffusive gas (w_(L)<<w_(G)) $\frac{df}{dt} = {{- {\gamma_{g}\left( {\delta \; v_{i}} \right)}}\frac{1.4256}{\sqrt{1 + {0.5{A_{d}\left( {\delta \; v_{i}} \right)}Z_{1{Ueff}}e^{{- {\gamma_{c}{({\delta \; v_{i}})}}}t}f}} + 0.4256}f}$
 4. The method according to claim 1, wherein the peak normalized line profile g(ν−ν₀,w_(R)) centered at the resonance frequency ν₀ of the molecular transition has either a Lorentzian shape, a Gaussian shape, or a Voigt shape depending on the experimental conditions.
 5. The method according to claim 1, wherein said fitting is a least squares fitting.
 6. The method according to claim 5, wherein said least square fitting uses the Levenberg-Marquardt (L-M) algorithm.
 7. The method according to claim 1, including: changing the frequency of the electromagnetic radiation emitted by the laser to a frequency ν_(i+d) where ν_(i+d) belongs to [ν_(min), ν_(max)] and repeating the steps of: fixing the intensity of said electromagnetic radiation in the cavity at a value much greater than the saturation intensity I_(s) of the molecular transition to be detected; irradiating said gas by means of said electromagnetic radiation beam emitted by said laser source having said fixed frequency ν_(i+d) and intensity in said resonant cavity; coupling said electromagnetic radiation to said cavity so as to obtain a laser-cavity resonance condition; changing the frequency of the electromagnetic radiation emitted by the laser so as to switch off the laser-cavity resonance; detecting an electromagnetic radiation beam in output from said cavity after the laser-cavity resonance has been switched off; recording a plurality of data representative of said output which has the form of a decay signal; fitting said recorded data with a function S^(repl)(t, δν_(i+d)) in which Z_(1Ueff)=constant replaces Z_(1U) g(ν−ν₀,w_(R)) in the fitting function s(t, δν_(i+d)).
 8. The method according to claim 7, including repeating the steps of claim 7, for a frequency ν_(i+2d)=ν_(i+d) of the electromagnetic radiation emitted by the laser, as long as the frequency of the electromagnetic radiation is included in [ν_(min), ν_(max)].
 9. The method according to claim 6, including: obtaining for each of a plurality of m frequencies ν_(j) with δν_(i) from a δν_(min) to a δν_(max) a value of γ_(g)(δν_(j)), and fitting said m values of γ_(g)(δν_(j)) so as to obtain a value of the first concentration of the target gas.
 10. The method according to claim 9, wherein said fitting includes: selecting as a free parameters in the fitting a parameter which takes into account the presence of other molecular absorptions l=1 . . . n in addition to the target molecular resonance l=0.
 11. The method according to claim 9, wherein said fitting includes: selecting as a free parameters in the fitting a parameter which takes into account the presence of a polynomial background around the resonance frequency of the target transition ν₀.
 12. The method according to claim 9, wherein said fitting is at least squares fitting.
 13. The method according to claim 12, wherein said least square fitting uses the Levenberg-Marquardt (L-M) algorithm.
 14. The method according to claim 9, wherein a result of said fitting is multiplied by a correcting factor R in order to obtain the concentration of the target gas, said correcting factor taking into account the line profile modifications due to the Z_(1Ueff)=constant approximation.
 15. The method according to claim 1, including the step of: selecting the value of Z_(1Ueff)=constant to be introduced in the fitting by means of the following step: inserting said gas at a second concentration, wherein said second concentration is at least 10 times said first concentration in the resonant cavity; tuning the frequency of said electromagnetic radiation emitted by said laser source so as to fix it to a value ν_(i) within a range of frequencies [ν_(min), ν_(max)] including said molecular transition ν₀; fixing the intensity of said electromagnetic radiation in the cavity at a value much greater than the saturation intensity I_(s) of the molecular transition to be detected; irradiating said gas by means of said electromagnetic radiation beam emitted by said laser source having said fixed frequency ν_(i) and intensity in said resonant cavity; coupling said electromagnetic radiation to said cavity so as to obtain a laser-cavity resonance condition; changing the frequency of the electromagnetic radiation emitted by the laser so as to switch off the laser-cavity resonance; detecting an electromagnetic radiation beam in output from said cavity after the laser-cavity resonance has been switched off; recording a plurality of data representative of said output which has the form of a decay signal; fitting the data of the recorded decay with a curve S(t, δν_(i)) which depends on the following parameters: B^(high conc)(δν_(i)) is the detection background, with (δν_(i)=ν_(i)−ν₀) A_(d) ^(high conc)(δν_(i)) is the amplitude of the decay signal at the beginning of the decay event, g γ_(c) ^(high conc)(δν_(i)) is the cavity decay rate due to non-resonant and non-saturable losses (empty cavity decay rate); γ_(g) ^(high conc)(δν_(i)) is contribution of the targeted molecular transition to the decay signal; Z_(1Ueff) ^(high conc)=Z_(1U) g(δν_(i),w_(R)), where g is the peak normalized line profile g(ν−ν₀,w_(R)) centered at the molecular resonance frequency ν₀ and w_(R) is the HWHM width of the resonance, and w_(R)=w_(L) for a Lorentzian shape, w_(R)=w_(G) for a Gaussian shape, w_(R)={w_(L), w_(G)} for a Voigt shape; $Z_{o} = {\frac{{CP}(0)}{{CP}_{s}} = {A_{d}Z_{1U}}}$ is the saturation parameter at the beginning of the decay event and at the frequency of the targeted molecular transition ν₀, P(0) is the intracavity power at the beginning of the decay signal; and $P_{s} = {\frac{\pi \; w^{2}}{2}I_{s}}$ is the saturation power, where w is the spot size radius of the laser beam, i.e. the radius for which the amplitude of the field is 1/e times that of the axis and I_(s) the saturation intensity; selecting as Z_(1Ueff) ^(high conc)=Z_(1Ueff)=value obtained from the above fitting.
 16. The method according to claim 14, wherein the value R is calculated performing a first and a second measurement of a concentration of the same target gas at the same temperature and pressure, the first measurement in saturation absorption and at a second concentration, wherein said second concentration is at least 10 times said first concentration, and said second measurement in linear absorption at said second concentration.
 17. The method according to claim 1, wherein said electromagnetic radiation emitted by said laser is in the Infrared range.
 18. A method of ring-down spectroscopy in saturated-absorption condition, for measuring a first concentration of a target gas through a measurement of the spectrum of a molecular transition of said target gas, the target gas being in a mixture together with other gasses the method comprising the steps of: repeating for a number m of fixed frequencies vi spaced each other by a of a frequency step in a range of frequencies including the frequency of the molecular transition of the target gas, the gases in the mixture having l=0, . . . , n absorptions in the measured spectral range, where l=0 is the target absorption, the following steps: repeating d times at the same frequency the following steps: inserting said gas whose first concentration is to be measured in a resonant cavity comprising two or more reflecting mirrors arranged so as to form a closed optical path for an electromagnetic radiation emitted by a laser source; tuning the frequency of said electromagnetic radiation emitted by said laser source so as to fix it to a value ν_(i) within a range of frequencies [ν_(min), ν_(max)] including said molecular transition ν₀; fixing the intensity of said electromagnetic radiation in the cavity at a value much greater than the saturation intensity I_(s) of the molecular transition to be detected; irradiating said gas by means of said electromagnetic radiation beam emitted by said laser source having said fixed frequency ν_(i) and intensity in said resonant cavity; coupling said electromagnetic radiation to said cavity so as to obtain a laser-cavity resonance condition; changing the frequency of the electromagnetic radiation emitted by the laser so as to switch off the laser-cavity resonance; detecting an electromagnetic radiation beam in output from said cavity after the laser-cavity resonance has been switched off; and recording a plurality of data representative of said output, obtaining a decay signal for the fixed frequency; collecting the d*m SCAR decay signals obtained; fitting at the same time the d*m decay signals with d*m fitting curves considering a fitting curve S(t; p(ν))=[S(t; p(ν_(i))_(j))]_(j=1) ^(d*m) for the recorded decay signals which depends on the following parameters p(ν)=[p(ν_(i))_(j)]_(j=1) ^(d*m): ${{p(v)} = {\left\lbrack {p\left( v_{i} \right)}_{j} \right\rbrack_{j = 1}^{d*m} = \left\lbrack {B_{j},A_{d_{j}},\gamma_{c_{j}},\left\lbrack \left\{ {{\gamma_{l}\left( {\delta \; v_{il}} \right)},{Z_{l\; 1U}{{\overset{\_}{g}}_{l}\left( {{\delta \; v_{il}},w_{Rl}} \right)}}} \right\}_{j} \right\rbrack_{l = 0}^{n}} \right\rbrack_{j = 1}^{d*m}}},\mspace{79mu} {i = {1\mspace{14mu} \ldots \mspace{14mu} m}}$ [B_(j)(δν_(il))]_(j=1) ^(d*m) are the detection backgrounds of the j=1 . . . d*m recorded SCAR decay signals; and δν_(il)=ν_(i)−ν_(l) is the detuning of the i (=1 . . . m) scanned frequency of the laser radiation with respect to the resonance frequency of the l resonant transition, [A_(d) _(j) (δν_(il))]_(j=1) ^(d*m) are the amplitudes of the j=1 . . . d*m recorded SCAR decay signals, [γ_(c)(δν_(il))]_(j=1) ^(d*m) are the cavity decay rates of the j=1 . . . d*m recorded SCAR decay signals due to non-resonant and non-saturated gas absorption losses, [[γ_(l)(δν_(il))_(j)]_(l=0) ^(n)]_(j=1) ^(d*m) are the arrays of the absorption decay rates of the l=0 . . . n saturated transitions at the δν_(il) detuning for the j=1 . . . d*m recorded SCAR decays, $\left\lbrack \left\lbrack \left\{ {Z_{1U\; l}{{\overset{\_}{g}}_{l}\left( {{v_{i} - v_{l}},w_{Rl}} \right)}} \right\}_{j} \right\rbrack_{l = 0}^{n} \right\rbrack_{j = 1}^{d*m}$ are the arrays of the l=0 . . . n saturation parameters of the j=1 . . . d*m recorded SCAR decay signals; g _(l)(ν_(i)−ν_(l),w_(Rl)) is the peak normalized line profile centered at the molecular resonance frequency ν_(i) and w_(Rl) is the HWHM width of the l-resonance; $Z_{o} = {\frac{{CP}(0)}{{CP}_{s}} = {A_{d}Z_{1U}}}$ is the saturation parameter at the beginning of the decay event and at the frequency of the targeted molecular transition ν₀, P(0) is the intracavity power at the beginning of the decay signal; and $P_{s} = {\frac{\pi \; w^{2}}{2}I_{s}}$ is the saturation power, where w is the spot size radius of the laser beam, i.e. the radius for which the amplitude of the field is 1/e times that of the axis and I_(s) the saturation intensity; dividing the above mentioned parameters in a first and a second group; and fitting said recorded data with a function S(t; p(ν))=[S(t; p(ν_(i))_(j))]_(j=1) ^(d*m) by keeping the first parameter group fixed and equal to a set of pre-determined values and performing a first fitting only considering the second group as free parameters in a first fitting step; and keeping the second parameter group fixed to a value given in the first fitting step and performing a second fitting considering only the second group as free parameters in a second fitting step.
 19. The method according to claim 18, including the step of calculating the first concentration of the target gas N₀ by: $N_{0} = \frac{{\gamma_{0}^{fit}(0)}{g_{0}^{fit}(0)}}{c^{2}{L_{S}(T)}}$ where c is the speed of light in vacuum and L_(s)(T) is the line-strength of the absorbent molecular transition at the temperature T.
 20. The method according to claim 18, including the step of parametrizing said S(t; p(ν))=[S(t; p(ν_(i))_(j))]_(j=1) ^(d*m) as: ${{S\left( {t;{p(v)}} \right)} = {\left\lbrack {S\left( {t;{p\left( v_{i} \right)}_{j}} \right)} \right\rbrack_{j = 1}^{d*m} = \left\lbrack {B_{j} + {A_{d_{j}}e^{{- \gamma_{c_{j}}}t}{f_{j}\left( {{t;A_{d_{j}}},\gamma_{c_{j}},\left\lbrack \left\{ {{\gamma_{l}\left( {\delta \; v_{il}} \right)},{Z_{l\; 1U}{{\overset{\_}{g}}_{l}\left( {{\delta \; v_{il}},w_{Rl}} \right)}}} \right\}_{j} \right\rbrack_{l = 0}^{n}} \right)}}} \right\rbrack_{j = 1}^{d*m}}},\mspace{79mu} {i = {1\mspace{14mu} \ldots \mspace{14mu} m}}$ with ƒ=[ƒ_(j)]_(j=1) ^(d*m) are the non-linear functions that follows one of the below rate equations depending of the gas conditions: homogeneous regime (w_(L)≥w_(G)) $\left\lbrack {\frac{{df}_{j}}{dt} = {- {\sum\limits_{l = 0}^{n}{{\gamma_{l}\left( {\delta \; v_{il}} \right)}_{j}\frac{\ln \left\lbrack {1 + {A_{d_{j}}Z_{l\; 1U}{{\overset{\_}{g}}_{l}\left( {{\delta \; v_{il}},w_{Rl}} \right)}e^{{- \gamma_{c_{j}}}t}f_{j}}} \right\rbrack}{A_{d_{j}}Z_{l\; 1U}{{\overset{\_}{g}}_{l}\left( {{\delta \; v_{il}},w_{Rl}} \right)}e^{{- \gamma_{c_{j}}}t}}}}}} \right\rbrack_{j = 1}^{d*m}$ inhomogeneous regime still diffusive gas (w_(L)<w_(G)) $\left\lbrack {\frac{{df}_{j}}{dt} = {- {\sum\limits_{l = 0}^{n}{{\gamma_{l}\left( {\delta \; v_{il}} \right)}_{j}\frac{2}{\sqrt{1 + {A_{d_{j}}Z_{l\; 1U}e^{{- \gamma_{c_{j}}}t}f_{j}}} + 1}f_{j}}}}} \right\rbrack_{j = 1}^{d*m}$ non-diffusive gas (w_(L)<<w_(G)) $\left\lbrack {\frac{{df}_{j}}{dt} = {- {\sum\limits_{l = 0}^{n}{{\gamma_{l}\left( {\delta \; v_{il}} \right)}_{j}\frac{1.4256}{\sqrt{1 + {0.5A_{d_{j}}Z_{l\; 1U}e^{{- \gamma_{c_{j}}}t}f_{j}}} + 0.4256}f_{j}}}}} \right\rbrack_{j = 1}^{d*m}$
 21. The method according to claim 18, including: approximating the parameters B_(j),A_(d) _(j) and γ_(c) _(j) as having small variations from decay to decay: B _(j) =B+δB _(j) S is the mean value of the background of the j= . . . d*m recorded SCAR decays, in U {δB_(j)}_(j=1) ^(d*m) are small variations around B for each j=0 . . . d*m recorded SCAR decay, in U A _(d) _(j) =Ā _(d) +δA _(d) _(j) Ā_(d) is the mean value of the amplitude of the j=0 . . . d*m recorded SCAR decays, in U {δ A_(d_(j))}_(j = 1)^(d * m) are small variations around Ā_(d) for each j=0 . . . d*m recorded SCAR decay, in U γ_(c) _(j) =γ_(c)+δγ_(c) _(j) γ _(c) is the mean value of the cavity decay rate of the j=0 . . . d*m recorded SCAR decays due to non-absorbent gas and non-saturated absorption losses, in s⁻¹ {δγ_(c) _(j) }_(j=1) ^(d*m) are small variations around γ _(c) for each j=0 . . . d*m recorded SCAR decay, in s⁻¹ selecting [p_(lc_(j))]_(j = 1)^(d * m) = [δ B_(j), δ A_(d_(j)), δ γ_(c_(j))]_(j = 1)^(d * m) as the first group of parameters.
 22. The method according to claim 21, including: introducing two further parameters to the second group of parameters, said parameters being: d is the contribution to the cavity decay rate due to residual detection non-linearities, in s⁻¹ Z_(d) is the “equivalent”-saturation parameter at 1U due to residual detection non-linearities, in U⁻¹ so that the second group of parameters is equal to: p _(gb)(ν)=[ B,Ā _(d),γ _(c) ,Z _(d) ,d,[Z _(1Ul),γ_(l)(0),g _(l)(0), g _(l)(δν_(il) ,w _(Rl))]_(l=0) ^(n)] Where γ_(l)(δν_(il))_(j)=γ_(l)(0)g _(l)(0) g _(l)(δν_(il) ,w _(Rl)) with {γ_(l)(0)}_(l=0) ^(n)={cα_(l)(0)}_(l=0) ^(n) are the absorption decay rates for each transition l=0 . . . n at the resonance frequency ν_(i)=ν_(l), in s⁻¹, {g_(l)(0)}_(l=0) ^(n) are the area normalization factors for each transition l=0 . . . n at the resonance frequency νi=ν_(l), {g _(l)(ν_(i)−ν_(l),w_(Rl))}_(l=0) ^(n) are the peak normalized line profiles for each transition l=0 . . . n, centered at the resonance frequency ν₁ and linewidth w_(Rl).
 23. The method according to claim 22, including a further second parameter [C_(l) ^(γ)]_(l=0) ^(n), which is the amplitude factor of a function G_(j)(ν_(i); [C_(l) ^(γ)]_(l=0) ^(n)), proportional to the spectrum profile [g _(l)(δν_(il),w_(Rl))]_(l=0) ^(n) G _(j)(ν_(i) ;[C _(l) ^(γ)]_(l=0) ^(n))=Z _(l=0) ^(n) C _(l) ^(γ) g _(l)(δν_(il) ,w _(Rl)), which is used to update the local parameters [δ γ_(c_(j))]_(j = 1)^(d * m).
 24. The method according to claim 22, including the step of parametrizing the non-linear function f as: homogeneous regime (w_(L)≥w_(G)) $\begin{bmatrix} {\frac{{df}_{j}}{dt} = {{- {\sum\limits_{l = 0}^{n}{{\gamma_{l}(0)}{g_{l}(0)}\frac{\ln \begin{bmatrix} {1 + \left( {{\overset{\_}{A}}_{d} + {\delta \; A_{d_{j}}}} \right)} \\ {Z_{l\; 1U}{{\overset{\_}{g}}_{l}\left( {{\delta \; v_{il}},w_{Rl}} \right)}e^{- {({{\overset{\_}{\gamma}}_{c} + {\delta \; \gamma_{c_{j}}}})}^{t}}f_{j}} \end{bmatrix}}{\left( {{\overset{\_}{A}}_{d} + {\delta \; A_{d_{j}}}} \right)Z_{l\; 1U}e^{- {({{\overset{\_}{\gamma}}_{c} + {\delta \; \gamma_{c_{j}}}})}^{t}}}}}} -}} \\ {d\frac{\ln \left\lbrack {1 + {\left( {{\overset{\_}{A}}_{d} + {\delta \; A_{d_{j}}}} \right)Z_{d}e^{- {({{\overset{\_}{\gamma}}_{c} + {\delta \; \gamma_{c_{j}}}})}^{t}}f_{j}}} \right\rbrack}{\left( {{\overset{\_}{A}}_{d} + {\delta \; A_{d_{j}}}} \right)Z_{d}e^{- {({{\overset{\_}{\gamma}}_{c} + {\delta \; \gamma_{c_{j}}}})}^{t}}}} \end{bmatrix}_{j = 1}^{d*m}$ and g _(l)(δν_(il),w_(Rl)) is a Voigt function inhomogeneous regime still diffusive gas (w_(L)<w_(G)) $\left\lbrack {\frac{{df}_{j}}{dt} = {\begin{bmatrix} {{- {\sum\limits_{l = 0}^{n}\frac{2\; {\gamma_{l}(0)}{g_{l}(0)}{{\overset{\_}{g}}_{l}\left( {{\delta \; v_{il}},w_{Rl}} \right)}}{\sqrt{1 + {\left( {{\overset{\_}{A}}_{d} + {\delta \; A_{d_{j}}}} \right)Z_{l\; 1U}e^{- {({{\overset{\_}{\gamma}}_{c} + {\delta \; \gamma_{c_{j}}}})}^{t}}f_{j}}} + 1}}} -} \\ \frac{2d}{\sqrt{1 + {\left( {{\overset{\_}{A}}_{d} + {\delta \; A_{d_{j}}}} \right)Z_{d}e^{- {({{\overset{\_}{\gamma}}_{c} + {\delta \; \gamma_{c_{j}}}})}^{t}}f_{j}}} + 1} \end{bmatrix}f_{j}}} \right\rbrack_{j = 1}^{d*m}$ and g _(l)(δν_(il),w_(Rl)) is a Gaussian function non-diffusive gas (w_(L)<<w_(G)) $\left\lbrack {\frac{{df}_{j}}{dt} = {\begin{bmatrix} {{- {\sum\limits_{l = 0}^{n}\frac{1.4256\; {\gamma_{l}(0)}{g_{l}(0)}{{\overset{\_}{g}}_{l}\left( {{\delta \; v_{il}},w_{Rl}} \right)}}{\sqrt{1 + {0.5\left( {{\overset{\_}{A}}_{d} + {\delta \; A_{d_{j}}}} \right)Z_{l\; 1U}e^{- {({{\overset{\_}{\gamma}}_{c} + {\delta \; \gamma_{c_{j}}}})}^{t}}f_{j}}} + 0.4256}}} -} \\ \frac{1.4256d}{\sqrt{1 + {0.5\left( {{\overset{\_}{A}}_{d} + {\delta \; A_{d_{j}}}} \right)Z_{d}e^{- {({{\overset{\_}{\gamma}}_{c} + {\delta \; \gamma_{c_{j}}}})}^{t}}f_{j}}} + 0.4256} \end{bmatrix}f_{j}}} \right\rbrack_{j = 1}^{d*m}$ and g _(l)(δν_(il),w_(Rl)) is a Gaussian function.
 25. The method according to claim 18, wherein said first fitting step is a least squares fitting.
 26. The method according to claim 25, wherein said least square fitting uses the Levenberg-Marquardt (L-M) algorithm.
 27. The method according to claim 18, wherein said second fitting step is a least squares fitting.
 28. The method according to claim 27, wherein said least square fitting uses the Levenberg-Marquardt (L-M) algorithm.
 29. The method according to claim 18, including the step of: dividing the second group of parameter in a first sub-group of free parameters and a second sub-group of calibration parameters which are set equal to calibration values which remain constant during the first and second fitting steps.
 30. The method according to claim 29, wherein the sub-group of calibration parameters includes: p _(gb) ^(cal)(ν)=[Z _(d) ,[Z _(1U1) ,g _(l)(δν_(il) ,w _(Rl))]_(l=0) ^(n)]^(cal) where {g _(l)(ν_(i)−ν_(l),w_(Rl))}_(l=0) ^(n) are the peak normalized line profiles for each transition l=0 . . . n, centered at the resonance frequency ν_(l) and linewidth w_(Rl), and Z_(d) is the “equivalent”-saturation parameter at 1U due to residual detection non-linearities, in U⁻¹. 